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PREFACE 


The three-dimensional rigid-lid model is Intended to be used for 
hydrothermal predictions of closed basins subjected to a heated discharge 
together with various other inflows and outffows. This volume has been 
written in order to assist any prospective user in applying the model to 
specific sites. Derivation of the governing equations and various other 
details have been omitted. The programs are fairly general and only 
one subroutine and a data file has to be rewritten for specific cases. 

This work was sponsored by the National Aeronautics and Space 
Administration (NASA-KSC) and the Environmental Protection Agency 
(EPA-RTP) . 


NOTICE 

The information in this document has been funded wholly or in 
part by the United States Environmental Protection Agency under 
assistance agreement 78-DX-0166 to the University of Miami, Coral 
Gables, Florida. It has been subject to the Agency's peer and 
administrative review and it has been approved for publication as 
an EPA document. Mention of trade names or commercial products does 
not constitute endorsement or recommendation for use. 



ABSTRACT 


The three-dimensional rigid-lid model was developed by the thermal 
pollution group at the University of Miami and verified for accuracy at 
various sites. The model results have been found to be fairly accurate 
in all the verification runs. The model is intended to be used as a 
predictive tool in futur a sites and this manual has been written to enable 
any user to be able to apply it without difficulty. 
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SECTION 1 
INTRODUCTION 


The need for mathematical modeling In predicting and monitoring 
thermal pollution was discussed In previous reports by Vezlroglu et al, 
(1973# 1974). Predictive studies of ecosystems can only be made by 
mathematical models. A prior knowledge of the effects of disturbances 
is essential for environmental impact studies. Thus# the mathematical 
model is a crucial tool In decisions Involving power plant siting# land 
development# etc. 

The University of Miami team undertook development of a methodology 
using remote sensing and numerical modeling to study thermal pollution. 
The use of remotely-sensed data in modeling has been discussed by Sen- 
gupta et al, (1974). The remote sensing effort has been discussed in 
detail In previous publications. This volume has been written so as to 
enable a user to apply the mathematical model to new sites for predictive 
purposes. 

The hydrodynamics and thermodynamics of an ecosystem are con- 
trolled by geometry, meteorological conditions and physical characteristics 
of the water such as density, salinity and turbidity. In this model the 
effects of salinity and turbidity have been neglected. Hence, the govern- 
ing equations are composed of the three-dimensional Navier-Stokes equa- 
tions and the energy equation. Various assumptions can be made for 
different situations leading to simplification or elimination of equations. 

The main simplifying assumption in this case is the rigid-lid assumption. 
This means that surface height fluctuations are not simulated Ly this 
model, and this is a reasonable assumption for most applications (e.g.. 
Lakes) . 

The rigid-lid model has the following capabilities: 

1. It predicts the wind-driven circulation. 

2. It predicts the circulation caused by inflows and outflows to the 
domain. 

3. It predicts the thermal effects in the domain. 

4. It combines the aforementioned processes. 

The calibration procedure consists of comparing ground-truth cor- 
rected airborne radiometer data with surface isotherms predicted by the 
model. 


SECTION 2 


RECOMMENDATIONS 


Various numerical models have been developed to study the effects 
of heated discharge and meteorological conditions on bodies of water. 

Most of these models are one or two dimensional. These models have a 
high computational speed but only give horizontally or vertically averaged 
values of temperatures. 

Three-dimensional models, however, have a much finer resolution 
but they consume larger computer time. The three-dimensional rigid- 
lid model can be used to obtain detailed temperature and velocity distri- 
butions in a domain where surface gravity waves are small compared to 
the depth of the domain. This model, as compared to free-surface mo- 
dels, runs faster since surface gravity waves are eliminated by the 
rigid-lid assumption. 

A proper method of using this model would be to run a one-dimen- 
sional model initially to obtain a rough picture of the temperatures and 
then using this model to obtain a better resolution, the 1-D results being 
used as ambient conditions. 

The following improvements have been suggested for the model. 

1. Since all natural flows are turbulent, proper turbulent closures are 
needed to make the model meaningful. At present, the simplest 
possible closures, namely constant eddy viscosities and eddy diffu- 
sivities, have been used. However, better results may be obtained 
by using a higher order closure. 

2. At present, the model uses uniform horizontal grids and stretched 
vertical grids. Nonuniform horizontal grids couid be introduced for 
better resolution near the boundaries. 

3. The program has been written to be run as a batch- job on the com- 
puter. It could be made interactive so as to enable the user to run 
it on a terminal. However, this would require some modifications in 
order to reduce the storage space. 
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SECTION 3 


PROCRAM DESCRIPTION AND FLOW CHART 


DESCRIPTION OF PROGRAM ALGORITHM 

The governing equations for a body of water which are derived from 
the basic laws of conservation of mass, momentum and energy are shown 
in Table 1. These equations incorporate a vertically-stretched coordinate 
system so as to make the model general enough to handle any kind of 
bottom topography. The problem is set up as an initial value problem. 
The initial values of the water velocities and temperatures are specified 
and the model is run so as to give the values of the above quantities !n 
subsequent time periods using an explicit scheme. The sequence of the 
calculations are as follows: 

1. The initial values of the velocities and temperatures are read into 
the program, the region of interest within the basin being classified 
into interior, corner or boundary points. (Subroutines used are 
READ 3K, INITIA, INITIT, HEIGHT.) 

2. The data, which includes the boundary conditions such as the 
various meteorological parameters like surface wind speed, air 
temperature, humidity and solar radiation are read into the program 
using subroutine READ 2. 

3. Depending on the site chosen, the various discharges (volume flow 
rate, velocities and temperatures) in and out of the basin are read 

*, into the model. These are incorporated in the subroutine INLET 1. 

4. The momentum, continuity and energy equations are now solved to 

determine the velocities and temperatures in the subsequent time 
steps. The predictive equation for pressure {viz., the Poisson 
equation) is solved iteratively to determine the pressures at various 
points of the domain. (Note: Because of the rigid-lid assumption, 

the surface or lid pressure Is no longer atmospheric.) 

THE PROGRAM FLOW CHART IS SHOWN IN FIGUFc 1 

The various subroutines used are as well as a brief description of 
their functions are shown in Tables 2 and 3. 

Symbols Used in Governing Equations 

(Quantities with bar are dimensional) 
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Table 1. Governing Equations 


Continuity Equation: 
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Hydrostatic Equation: 

Ty = E u (,+P)h 

Energy Equation: 
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SECTION 4 


LIST OF PROGRAM SYMBOLS USED IN MAIN PROGRAM 
DESCRIPTION OF MAIN VARIABLES 

A. A - constant In equation of state, p = A + BT + CT 2 
AREF - reference eddy viscosity 

AA “ value of 'V' at plume Inlet 

ABR - t /Rossby number 

AH - 1 /Reynolds number 

Al - coefficient in front of pressure term 

AKT - (K s )(H ref )/(B z ) 

AP - coefficient in front of pressure term 
ARBP - arbitrary pressure 

AV ' 7'ib where e= E 

E 

A3 - normalized vertical eddy coefficient of viscosity 
ANGLE - wind direction angle 

B. B - constant in equation of state, p = A + BT + CT 2 

BB - value of 'V' at plume inlet (at 1=10) 

BZ - pC p B y 

BV - normalized vertical eddy diffusivity, normalized with respect 
to reference eddy diffusivity 

C. C - constant ir equation of state, p = A + BT + CT 2 

CC - value of y (constant) 

CW - temperature gradient at vertical boundaries 
CB - temperature gradient at the bottom 
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D - U at previous time step 


D1TZ 

3T 
" TT 


DPX - 

3P 

dX 


DPY - 

3P 

3y 



3P 


DPSX 

s 

ax. 



3 P_ 


DPSY 

9 

3y 


DT - 

time increment 

DX - 

increment in 

x-direction 

DY - 

increment in 

y-dlrection 

DZ - 

increment in 

Z-direction 

D1HUX 

3(hu) 

3 X 



D1HVY - 3 (hv) 

9 y 

DiHUUX - 

d X 

d,huvy - iaaud 

D 1 HV V Y - ■- A" 1 

d y 

D1UY - 

d y 

D1VX - |£ 

D2UX - 0 
02VX - §£- 

divwx - 

D1UZ - ||j. 
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D2UZ 


3 2 u 

nr 


D1VZ - 


D2VZ - 


a 2 v 

TXT 


D1A3Z 


DLZ - 


[DX 2 ) (DY 2 ) 


E - V at previous time step 

EPS - convergence criterion 

EUL - Enler number 

EX - residual error in pressure iteration 

FH - forcing function in pressure equation 

FW - factor in wind stress calculation formula 

G - dummy variable for V (for future time step) 

H - dummy variable for U (for future time step) 


HI - 


HREF 


HX - 


HY - 


nondlmensional depth * |4 
- reference depth 


IN - maximum number of grid points in x-direction 

IWN - maximum number of half-grid points in x-direction„ IWN 
IN - 1 


i - index of x-axis, main grid 

ITN - index for number of iterations 

IW - index for x-axis, half grid 

IRUN - index for number of runs 
= 0, first run 

= 1, from second time onwards 
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JSGNX, ISGNY - determine signs of TAUX and TAUX respectively 

J. J - index for y-axis, main grid 
JW - index for y-axis, half grid 

JWN - maximum number of half-grid points In y-dlrection 
JWN - JN - 1 

JN - maximum number of main grid points in y-direction 

K. K - index for Z-axis 

KSTORE - specified usage of tape for storing results 

KN - maximum number of main grid points in Z-directlon 

KISS - surface heat transfer coefficient (nondimensional) 

L - maximum length of the domain 

LN - number of time steps to be computed 

LLN - total number of time steps /LN 

M. MAR - number to describe general location of a point in the main 
grid 

MRH - number of describe general location of a point in the half 
grid 

MAXIT - maximum number of iterations 

O. OMEGA - relaxation factor 

P. P - nondimensional pressure 

PN - New pressure, nondimensional 

PiNTH - dummy variable for pressure (future time step) 

R. R - dimensional density at main grid points 
RE - Reynolds number 

RB - Rossby number 

RINTX - density integrated with respect to x 

RINTY - density integrated with respect to y 
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res 








RO - nondimensional density at main grid points 
ROW - nondimensional density at half grid points 
RREF - reference density (gm/cc) 

RW - dimensional density at half grid points (gm/cc) 

RADN - solar radiation (w/m 2 ) 

T. T - nondimensional temperature at main grid points 
TO - initial temperature (dimensional) (°C) 

TAMB - ambient temperature (dimensional) (°C) 

TAIR - air temperature (dimensional) (°C) 

TAI - coefficient in front of convective terms In the energy equa- 
tion, = 1. 

1 

TAH - ur~ where P„ = R x P^ 

P e e e r 

TAV - p where E 

e 

TE - equilibrium temperature (dimensional) (°C) 

TTOT - total time elapsed 
TAUX - 3u/3y (nondimensional) 

TAUY - 3v/3y (nondimensional) 

TEM - dimensional temperature at main grid points 
TEMW - dimensional temperature at half-grid points 
TREF - reference temperature 

TW - nondimensional temperature at half-grid points 
TLL - temperature at the discharge point (nondimensional) 

TSU - water surface temperature (nondimensional) 

TDEW - dewpoint temperature (dimensional) 

U. U - velocity in x-direction (nondimensional) 

7. V - velocity in y-direction (nondimensional) 
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VVIS - vertical eddy viscosity (nondlmenslonal) 

W. W - velocity In Z-dlrection (nondlmenslonal) 

WH - W at half-grid points 

WHLDT - time derivative of WH at Hd (l.e., |*g(WH)/Z = 0) 

X. XI NT - integral of x terms on the right-hand side of Poisson's 

equation 

X - horizontal coordinate across discharge 

Y. YINT - integral of y terms on the right-hand side of Poisson's 

equation 

Y - horizontal coordinate across discharge 

Z. Z - vertical coordinate 
MARKER MATRICES 

The following number convention is used for the MAR = matrix sys 
tern, which classifies points (or nodes) on the main grid system = 

(Refer to Figure). 

MAR = 0, points outside the region of interest. 

MAR = 1, point on the far y-boundary. 

MAP. = 2, point on the near y-boundary. 

MAR = 3, point on the near x-boundary. 

MAR - 4, point on the far x-boundary. 

MAR = 5, outside corner on near x-boundary and far y-boundary. 

MAR » 6, inside corner on far x-boundary and far y-boundary. 

MAR a 7, outside corner on near x-boundary and near y-boundary. 
MAR = 8, inside corner on near x-boundary and near y-boundary. 

MAR = 9, outside corner on far x-boundary and near y-boundary. 

MAR s 10/ outside corner on far x-boundary and far y-boundary. 

MAR = 11, points in the interior of the region of interest. 
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trlx 

MRH 

MRH 

MRH 

MRH 

MRH 

MRH 

MRH 


The following number convention is used to describe the MRH (ma~ 
for the half-grid system). 

= 1, corner at far x-boundary and far y-boundary. 

- 2, points on near y-boundary. 

= 3, points on near x-boundary. 

= 4, corner at near x and near y-boundaries. 

= 6, far corner on x-axis. 
a 7, corner at far x and y-boundaries. 

= 9, interior grid points. 


12 




MOTH 


SECTION 5 

PREPARATION OF RUNS 


This section presents the steps to be followed In order to run the 

mode! for a particular location. 

1. The boundaries are chosen depending on the particular situation , the 
general Idea being to Include all inflows and outflows. If a heated 
discharge enters the body of water the region of Interest must be 
chosen so as to include this since It Is a major factor In determining 
the size and spread of the resulting plume. 

2. The grid size Is chosen depending on the resolution required. The 
user should remember that the choice of the grid size directly deter- 
mines the maximum allowable time step since this Is directly re- 
lated by the various stability criteria. (See choice of time step in 
Section 6.) 

3. Specify number of full-grid points IN, JN, KN and number of half- 
grid points IWN, JWN. Since the actual domain may be smaller than 
the total rectangular region, INxJNxKN, the marker matrices MAR 
and MRH are used to specify the domain so that points outside the 
domain of interest skip the subsequent calculations. 

4. IRUN is specified (= 0 for the first run, * t for subsequent runs). 
KSTORE is specified to indicate whether any tape has been assigned 
to store results of the run. 

KSTORE » 0 if no tape has been assigned. 

~ 1 if tape has been assigned. 

LLN is specified to denote the number of hours of simulation 1,o be 
carried out. 

5. The depths at various places within the domain are specified using 
subroutine HEIGHT. The various Inflows and outflows to the domain 
are specified using INLET t. (For details please refer to Biscayne 
Bay run, Sengupta et ai. (1975).) 

6. The various data like solar radiation, wind speed, wind direction and 
dewpoint temperature are specified in a data file which is made by 
the main program. 

For further details see the next section. 
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SECTION 6 
INPUT DATA 


The data that is required for the execution of the main program Is 
listed in Table 3 in the order It appears. Note, the data input symbols 
have already been defined in Section 4. Moreover, the following remarks 
should be observed. 

* Free format Is used for all data Input. 

* Distinction must be made for integer and real number. 

* The order of the cards must be followed. 


SECTION 7 


PLOTTING PROGRAMS 


The plotting programs for the 3-D rigid-lid model are distinct from 
the main program and subroutines used to run it. The user has an 
option of either using a tape (Unit 8) during running the main program 
TMAINN to store the results or )ust run it without storing the results. 

For making subsequent continuation runs of TMAINN all that is required 
Is the result of the last hour In the previous run. For plotting, however, 
one needs the results of all the hours for which results are to be plotted. 
These results are used as Input data to run the various plotting programs. 

DESCRIPTION OF PLOT PROGRAMS 

The following are the main plotting programs. 

PLOT - plots surface isotherms- 

PLUV - plots u, v components of the velocities (I.e./'K* sections). 

PLUW - plots u, w components of the ve?«,i titles (i.e. , 'j' sections). 

PLVW - plots v, w components of the velocities (i.e., 'I' section/?'. 

SUBROUTINES 

The various plot programs and subroutines are shown in Table 4. 

Other subroutines seen in these programs (e.g., ARROHD, FLINE, 
etc.) are standard FORTRAN subroutines used for plotting, using a 
CALCOMP x,y plotter, and are hence ommitted in the above listing. 
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APPENDIX A 


EXAMPLE CASE 


INTRODUCTION 

The area of interest Is Lake Keowee in South Carolina, which was 
formed from 1968 through 1971 by damming the Little and Keowee rivers. 
The lake is located about 40 km wist of Greenville and constitutes Duke 
Power Company's Keowee-Toxaway complex. 

Lake Keowee has two arms connected by a canal (maximum depth 
30.5 m) . There are three power plants on the lake, namely, the Oconee 
Nuclear Station, Keowee hydro station and Jocassee-pumped storage 
station. The Oconee Nuclear Station is a three unit steam-electric sta- 
tion with an installed capacity of generating 2580 MW. The Oconee 
Nuclear Station draws in condenser -ceoiing water from the lower arm 
of Lalke Keowee and discharges the heated effluent to the upper arm of 
the lake. The intake structure for the condenser-cooling water allows 
water from 20 to 27 m depth (full pond) to pass through. The discharge 
structure has an opening from 9 to 12 meters below the water surface 
(full pond) through which the CCW returns directly to the upper branch 
of the lake. 

Lake Jocassee is located north of Lake Keowee and is used as a reser 
voir f^r Jocassee- pumped storage station. Lake Keowee also serves as 
the (oiver pond for this station. The Jocassee station has reversible 
turbines with a maximum generating flow (into Lake Keowee) of about 
820 m 3 /sec and a maximum pumping flow (out of Lake Keowee into Lake 
Jocassee) of about 775 m 3 /sec, the net flow into Lake Keowee from Jocas- 
see being about 15.5 m 3 /sec. 

Lake Keowee has a full pond elevation of 243.8 m above MSL. At 
full pond it has a volume of approximately 1.18 x 10* m 3 , an area of 
74 km 2 , a mean depth of 15. 8 m and a shoreline of about 480 km. The 
outflow from Lake Keowe* is through Keowee hydro station and may vary 
from approximately 1.4 m- /sec (leakage) to 560 m 3 /sec. Maximum allow- 
able draw-down of the lake is 7.6 m. 

A map of the area of interest is shown in Figure 3. 

PROBLEM STATEMENT 

The object’ve of the present work is to find the three-dimensional 
temperature and velocity distributions in the region where the effects 
of the thermal discharge are noticeable. The effects of Jocassee-pumped 


storage station, Keowee hydro station as wall as the meteorological condi- 
tions have been incorporated. 

The region of Interest Is chosen 'to include the effects of the Oconee 
Nuclear Station discharge, the outflow through Keowee dam and the impact 
of the J ocas see- pumped storage station on the velocity and temperature 
distributions In Lake Keowee. The depth of the domain is cut off at 16 
meters, since this is the level at which the thermocline occurs. Hence, 
for running the model, a constant depth region is considered. The plan 
view of the domain is shown in Figure 4. (Note: For variable depth 

refer to Biscay ne Bay simulation studies by the University of Miami 
thermal pollution group.) In this figure, AB is an open boundary which 
takes care of the flow from or to the Jocassee-pumped storage station. 

*C' shows the position of the flow In the canal connecting the two arms 
of the lake. 'D' is the discharge point for the Oconee Nuclear Station 
and 'E* is the outflow from Keowee hydro station. 

The inclusion of the above results in a domain 2895.6 m x 2438.4 m 
in the horizontal plane. The horizontal grid size (in x and y directions) 
is 152. 4 m x 152. 4 m, giving a total of 20 x 17 (= 340) nodes in the 
horizontal plane, out of which 293 lie in the region of interest. The 
16 m constant depth region of interest is divided into 4 equal slices of 

4 m each, giving at total of 5 nodes in the vertical (Z) /direction. 

Hence, there are 293 x 5 nodes (grid points) in the region of interest. 
This region is specified using the MAR and MRH marker matrices (Figure 

5 and Figure 6). 

Boundary Conditions 

On the Jocassee effect boundary, the flowL .velocity (varying with 
time) is specified. Open-boundary condition (|— - = 0) is specified for 
the temperature. y 

The same is done for the Keowee hydro boundary. The only differ- 
ence is that the values specified are at three points in the vertical plane 
(i.e. , at K = 1, 2 and 3) since this region covers the discharge area. 

For the Oconee Nuclear Station, the discharge velocity as well as 
the discharge temperature is specified at the discharge point. 

Open-boundary conditions are specified for the temperature and 
velocity at the canal. This, however, leads to a possible violation of 
mass balance in the region of interest. This mass unbalance will 
actually show up as a variation in the water level in the lake which is 
beyond the capability of the rigid-lid model. . 

At all solid boundaries as well as the artificial bottom (since the 
bottom Is cut off at 16 m) perfect insulation (temperature gradient = 0) 
and zero velocity conditions are assumed. 

At the surface, the vertical component of the velocity is specified 
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as zero (rigid-lid constraint) . Surface wind shear stress and heat 
transfer coefficient are specified. 

Initial Conditions 

The Initial values of the water velocities are assumed to be zero. 
The initial temperature of the lake is assumed to be equal to the am- 
bient water temperature (determined by running a one-dimensional 
model) and is taken to be uniform throughout the domain. 

CALCULATION OF PARAMETERS AND INPUT DATA 

Reference Quantities 

Reference length = L = maximum length of the domain = 2895.6 m. 
Reference horizontal eddy viscosity A re ^ = 0.002 L 4/3 

= 38311. 48 cm 2 sec. 

For better agreement with data the value chosen Is 60,000 cm 2 /sec. 

Reference depth = H = 16 m. 

Reference vertical A y = 0.002 x (H) 4/3. 

Eddy viscosity = 37.43 cm 2 /sec. 

Reference velocity = V re ^. = 30 cm /sec. 

Reference temperature =* T re f = 10.0°C. 

Reference time = L/V re f = 9652 sec. 

Calculation of Inflows and Outflows Into the Domain (Used in INLET1) 

Oconee Nuclear Station Discharge Velocity — 

The discharge is considered to take place through a point at a 
depth of 12 m (k = 3). The discharge velocity is calculated as follows: 


12 m 


The total discharge into the basin is equal to: 
(100 £2 x V x 152.4 x 12) a Q 


h 




152.4 m 
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where Q = average discharge irr m 3 /sec 


.\V » 


8144.1 

~W~ m/MC 


= 7.42207 cm /sec 


The average values of Q over 24 hrs Is taken since the variation 
is negligible. 


Nondimensional discharge velocity = 


ref 


s l0 


a 0. 24740 


Keowee Hydro Discharge Velocity — 

The outflow through the Keowee hydro station is through a channel 
152.4 m x 12 m. 


The volume flowrate Q = (152. 4 x 12 x V) m 3 /sec 
where V a discharge velocity (m/sec) 

. * .V = [Q/(152. 4 x 12)] m/sec = ( a ^YZx lO O ' cm/se€ 

Q is specified as a function of time in INLET 1. 

The procedure for nondimensional! zation is similar. 

Jocassee Flow Velocity — 

The entire flow to or from the J ocas see-pumped storage station is 
assumed to take place through the entire upper boundary (AB in Figure 
4). The flow through this area (shown below) is assumed to be uniform 
and is assumed to take place simultaneously with the outflow through the 
Jocassee station. 


16 m 


r 

V a q/[ (16x1 3x152. 4) x 100] cm /sec. 

Q a flow through Jocassee (m 3 /sec). 

Q is positive when Jocassee is generating (i.e., the flow is into 
the region of interest) and negative when pumping (i.e./ flow out of 
region of interest). 
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SAMPLE INPUT 


The following are the inputs to TMAINN contained in the data file 
IPUT (which includes values calculated earlier) . 


Input 

# 

No. of Data 
In Card 

Symbol 


Value 

1 

3 

(RUN 

s 

0 



KSTORE 

3 

1 



LLN 

3 

3 

2 

2 

VVIS 

3 

37.43/60,000 = 0.00062384 



ABR 

3 

0.78 

3 

4 

Al 

3 

1.0 



AH 


60,000 _ . 01228172 




So x 2 6 9 S x loo 2 



AV 

= 

(209|iiij AH _ 402.08304 



AP 

3 

1.0 

4 

4 

EPS 

3 

0.001 



MAXIT 

3 

60 



OMEGA 

3 

1.8 



ARBP 

3 

1.0 

5 

3 

DX 

3 

152. 4/2895.6 = 0.05263 



DY 

= 

0.05263 



DZ 

3 

4/16 = 0.25 

6 

3 

TAI 

= 

1.0 



TAH 

3 

AH = 0.01228172 



TAV 

= 

AV = 402.08304 

7 

3 

A 

3 

1.000428 



B 

3 

-0.000019 
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Input 

# 

No. of Data 
In Card 

Symbol 


Value 



C 

3 

-0.0000046 

8 

1 

TO 

3 

10.0 

9 

3 

EUL 

= 

98 9 iLijlMM - 1742.222 



CW 

SI 

0.0 



CB 

ss 

0.0 

10 

2 

AA 

=■' 

0. 24740 



CC 

3 

16/16 = 1.0 

11 

1 

TLL 

3 

31.7 - 10 

— ns — 

12 

1 

TAU 

3 

0.0152 cm 2 /sec 

13 

1 

DT 


Criterion (convective) 




3 

. Ax 152.4 x 100 

4t * U " 30 





» 504 secs > 504 secs 





Hence, convective criterion 
dominates; choose AT = 300 secs 





DT = a = 0.03108164 





Note: choose best time step by 





trial and error 

14 

1 

CTTOT 

= 

Vef /360 ° = 2 * 68ir,tl 

15 

1 

ISOTOP 

55 

0 

16 

6 

WS 





TSU 





TD'HW 


See Table 5 



RADN 





ISGNX 





I5GNY 



17 

1 

ANGILF. 


See Table 5 
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LAKE KEOWEE APPLICATION-EXECUTION DECK 

The following execution dock Is for use In the UN I VAC 1100 computer 
et the University of Miami. These may have to be modified if a different 
computer is used. 

(ALL PROGRAMS AND SUBPROGRAMS COMPILED AND STORED IN FILE) 
First Run 

1. 0 ASG, AX FILE. 

(THE FILE IS ASSIGNED FOR THE RUN) 

2. 9 ASG,T 8, 16N, TAPENAME. 

(A TAPE FILE NAMES '8' IS BEING ASSIGNED. THE TAPE IS 
9-TRACK, AND THE REEL NUMBER IS 'TAPENAME') 

3. @ PRT,S FILE. TMAINN 

(THE MAIN PROGRAM IS PRINTED) 

4. 9 PACK FILE. 

(THE FILE IS PACKED) 

5. 9 PREP FILE. 

(ENTRY POINT TABLE IS PREPARED) 

6. 9 MAP,S 

7. IN FILE. TMAINN 

8. LIB FILE. 

9. END 

10. 9 XQT 

11. 0 

(VALUE FOR JR UN, FIRST RUN: IRUN=0) 

12. 24 

(NUMBER OF HOURS REQUIRED, MINIMUM=1 HOUR, MAX-24) 
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iw-iiem^iipeiw 


13. 0 

(0 IF MAGNETIC TAPE IS REQUIRED TO STORE RESULT* IF 
NOT* ANY NUMBER) 

14. 0 ADD FILE. INPUT 

(INPUT DATA FILE FOR THE PARTICULAR RUN) 

15. 0 FIN 

EXECUTION DECK FOR PLOT PROGRAMS 

1. 0 ASG.AX FILE. 

2. 9 ASG,T 8., 16N, TAPENAME. 

3. & ASG*T 11., 16N, PLOTTAPE. 

(A MAGNETIC TAPE FILE NAMED '11' IS BEING ASSIGNED. THE 
TAPE IS 7-TRACK AND THE REEL NUMBER IS 'PLOTTAPE'. THE 
PLOTS ARE STORED ON THIS TAPE) 

4. 9 PRT*S FILE. PLOTTER 

(THE PLOT PROGRAM IS PRINTED) 

5. 9 PACK FILE. 

6. 0 PREP FILE. 

7. 0 MAP,S 

8. IN FILE. PLOTTER 

9. LIB FILE. 

10. END 

11. 0 XQT 

12. 0 ADD FILE. INPUT 

13. 0 FIN 
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Tibia 2. Subroutine Required In Main Program TMAINN 


No. 

Name 

Description 

Remarks. . 

t 

DVISV 

Computes D1VY, D2VY, D1VX 
and D2VX. 

Called by subroutine INTE. 
Schemes used similar to 
DVISU. 

2 

DVISU 

Computes D1UX, D 2 UX, and 
D1UY. 

Called by INTE. f^, 
are computed at Interior, 
boundary or corner pts 
by scherrie similar tc the 
one used in DINERU. 

3 

DVVY 

Computes D1HVVY. 

i 

Called by INTE. ~ (hvv) 
is computed for interior, 
boundary or corner by a 
scheme similar to the one 
used in DINERU. 

4 

DUVY 

Computes D1HUVY. 

Called by INTE. ~ (huv) 
is computed for interior, 
boundary and corner pts 
by a scheme similar to the 
one used in DINERU. 

5 

DINERU 

Computes D1HUUX and 
D1HUVX. 

Called by INTE. The re- 
sults are used in Poisson 
equation for pressure. 

6 

TPRINK 

l 

Prints temperatures at a grid 
point. 

Called by TMAINN. 

7 

PRUV 

Prints the values of U and V 
at all main grid points. 

Called by TMAINN. 

8 

PRiTEX 

Prints the No. of iterations 
(ITN) and final residual error 
in solving the Poisson equation 

Called by TMAINN. 

9 

TPRIN1 j 

Prints the input para meters. 

Called by TMAINN. 

10 

STORE 2 

Stores values of input para- 
meters and physical quantities 
on tape #8 

Called by TMAINN. 

IT 

1 

RWR 

Computes real vertical veloci- 
ties from modified vertical 
velocities used in equations at 
integral grid points. 

Called by TMAINN. 
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Table 2. Subroutines Required in Main Program TMAINN (Continued) 


No. 

Name 

...1 1 l 1 in n I' f nn •' 'IT 1 1 — i 1 ™ 1 ■ 1 " l 

Description 

Remarks 

12 

RWRH 

Computes real vertical veloci- 
ties at half-grid points. 

Called by TMAINN. 

13 

DENSTY 

Uses the equation of state and 
computes density field from the 
temperature field. 

Called by TMAINN. 

14 

TEQB 

Allows for vertical mixing at a 
particular grid point. Program 
Is called by TMAINN. 

If the temp at the grid 
pt just above it is (ess 
and the difference is more 
than a specified maximum, 
the two temperatures are 
averaged. 

15 

OLDT 

Sets the values of temperature 
field at time step 'n' equal to 
the temperature field at (n+1) 
after all computations for time 
step f n s are completed. 


16 

TEMB2 

Computes temperatures at the 
boundary points in the domain 
of interest. 

Called by TMAINN. 

17 

TEMI4 

Computes temperatures at the 
interior points of the domain of 
interest. 

Called by TMAINN. 

18 

RWH 

Computes vertical velocities at 
half-grid points. 

Called by TMAINN. 

19 

OLDUV 

Sets the values of D and E 
equal to U and V respectively 
in order to retain values of U 
and V at one time step lag. 

Called by TMAINN. 

28 

UVTOP 

Computes U and V at the top 
using wind stress boundary 
conditions. 

Called by TMAINN. Com- 
putations are made for 
MAR = 11 only (internal 
grid points). 

21 

UVT 

Computes U and V for variable 
density at successive time 
steps. 

Called by TMAINN. 
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Table 2, Subroutines Required In Mein Program TMAINN (Continued) 


No. 

Name 

Description 

Remarks 

22 

PRE1L 

Computes pressure for far 
field from Poisson's Equation 
at half-grid points. 

Called by TMAINN. 

23 

FORCE 

Computes R.H.S. of Poisson's 
Equation at half-grid points. 

Called by TMAINN. 

24 

DPSXY 

Computes DPSX and DPSY. 

Called by TMAINN. 

25 

ROINTY 

Computes Y In the Poisson's 
Equation. p 

Called by TMAINN. 

26 

ROINTX 

Computes X in the Poisson's 
Equation. p 

Called by TMAINN. 

27 

CORINT 

Adds integral of Coriolis' 
component XI NT and YINT. 

Called by TMAINN. 

28 

INTE 

Computes XI NT, YINT, DPSX, 
and DPSY. 

Called by INTE. 

29 

WHATIJ 

Computes the values of W at 
I, i from the values of WH at 
IW, JW. 

1 Called by TMAINN. 

30 

WHTOP 

Sets the value of WH equal to 
zero at the surface. 

Called by TMAINN. 

31 

ERROR 

. 

• 

Calculates "Hlrt and Harlow" 
correction term at half-grid 
points and at the surface 
(WHLDT). 

Called by TMAINN. 

32 

READ 2 

Reads in input parameters and 
physical quantities stored on 
tape #7. 

Corresponds to store 2 
Called in by TMAINN. 

33 

INLET 1 

I 

Puts in velocities u and v 
pheme discharge, etc. into 
the model. 

Called by TMAINN. 

34 

HEIGHT 

Inputs depths of the basin 
into the model. 

This subroutine is for i 
constant depths model. 
Called by TMAINN, 
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35 INITIT Sets initial temperature field. Sets the temperature field 

equal to ref temp s t all 
grid points. 

Called by TMAINN. 

36 INITIA Initializes values of U, V, WH, Called by TMAINN. 

W, D, E and PINTH. 

37 READ3K Classifies region of interest Called by TMAINN. 

Into interior, corner and 
boundary points using matrix 
MAR. 
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IPUT 


Data files containing values of 
Input data for the respective 
days. 







Read data cards 


First run? (I RUN = 0) 


READ 2 - Reads tape for 
data of previous run 


INLET 1 - Inputs V, U, 
TD to start run 


READ3K - Classifies the re- 
gion into interior, corner or 
boundary points. 

I NIT I A - Initializes U, V, W, 
WH, D, E and P. 

INITIT - Initializes T and p. 
HEIGHT - Inputs depth. 
INLET1 - Inputs V, U, TD 
TD to start run. 


ERROR 

. 

Computes WHLDT 

WHTOP 

- 

Sets WH equal to zero 

WHATIJ 

- 

Computes W 

INTE 

- 

Computes XI NT, YINT, DPSX and DPSY 

ROINTX 

- 

Computes Xp 

ROINTY 

- 

Computes Yp 

DPSXY 

- 

Computes DPSX, DPSY 

FORCE 

- 

Computes FH 

PRE1L 

- 

Computes PN 

UVTOP 


Computes U and V at the top 

OLDUV 

- 

Sets D = U and E = V 

RWH 

- 

Computes vertical velocities at half-grid points 

WHATIJ 

- 

Computes W at grid points 

TEMI4 

- 

Computes temperatures at interior points 

TEMB2 

- 

Computes temperatures at boundary points 

OLDT 

- 

Updates temperatures 

TEQB 

- 

Allows for vertical mixing 

DENSTY 

- 

Computes density field 

INLET! 


Inputs V, U and TD every time step 


A-,new- hour? 


, Yes 


RWRH - Computes real vertical velocities 

RWR - Computes real vertical velocities 

STORE 2 - Stores computed results on magnetic tape (unit 8) 

TPRIN1 - Prints input parameters 

PRITEX - Prints number of iterations and error 


Figure 1. Flow chart (main program) 
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Figure 1 (Continued) . Flow chart (main program) 
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Table 3. Input Data to TMAINN 


Input 

# 

No. of Data 
In Card 

Symbol 

Definition /Value 

1 

3 

IRUN 

s 

0 for first run 



LLN 

= 

No of hours of simulation 



KSTORE 

= 

0 If no tape is assigned 

1 if tape is assigned 

2 

2 

VVIS 

= 

Nondimensional vertical eddy 
viscosity 

3 

4 

ABR 

Al 

s 

1 /Rossby No. = ir-~ 
u ref 

Coefficient in front of inertia 
term = 1.0 



AH 

= 

1 /Reynolds No. = 

Ref eddy hoz viscosity 

U . . L ~ 

ref 


4 

AV 

s 

(1 /e 2 Re) (e = H/L) 



AP 

= 

Coefficient in fron of pressure 
term =1.0 

4 

4 

EPS 

= 

Convergence factor = 0.001 



MAXIT 

= 

Maximum number of iterations 
for Poisson Equation 



OMEGA 

= 

Relaxation factor = 1.8 


| 

ARBP 

= 

Arbitrary pressure = 1.0 

5 

3 

DX 

= 

Horizontal grid spacing (x dir.) 



DY 

= 

Horizontal grid spacing (y dir.) 
Ay/L 



DZ 

= 

Vertical grid spacing (z dir.) 
Az/H 

6 

3 

TAI 

= 

Coefficient of convective terms 
in energy equation = 1.0 



TAH 

= 

Horizontal eddy diffusivity 
AH (usually) 





Tibia 3. Input Diti to TMAINN (Continued) 


Input 

» 

No. of Date 
In Cird 

Symbol 

Definition/Value 



TAV 

= 

Vertical eddy diffusivity 




= 

AV (usually) 

7 

3 

A 

55 

1.000428 These are coefficients 




= 

-0.000019 in the equation of 





-0.0000046 state for water where 
p = A + BT + CT 2 (gm/cc) 

8 

1 

TO 

- 

Reference temperature (°C) 

9 

3 

EUL 

= 

' rer 



CW 


Temperature gradient at vertical 
boundary 



CB 

B 

Temperature gradient at the 
bottom 

10 

2 

AA 

s 

Nondimensionai discharge velocity 




B 

(discharge velocity) AJ re f 



CC 

= 

No dimensional depth = h/H^ 

11 

1 

TLL 

B 

Nondimensionai discharge tem- 
perature = (T d - T ) /T 

12 

1 

TAU 

= 

Surface shear stress (from Wilson 
Curve) (Refer to Figure 7) 

13 

1 

DT 

- 

Nondimensionai time step 




= 

4T(L/U ref ) 

14 

1 

CTTOT 

- 

Converts nondimensionai time to 
hours 

15 

1 

ISTOP 

= 

Number of hours of previous run 

16 

6 

WS 

= 

Wind spoed (m/sec) 



TSU 

- 

Air temperature (°C) 



TDEW 

= 

Dewpoint temperature (°C) 



RADN 

b 

Incident solar radiation (w/m 2 ) 
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Table 3. Input Data to TMA1NN (Continued) 


Input 

* 

No. of Data 
In Card 

Symbol 

Definition /Value 



ISGNX 

= +1 If x component of W Is 
negative 

= -1 If x component of W is 
positive 

= +1 if y component of W is 
negative 

= -1 If y component of W is 
positive 

17 

1 

ANGLE 

= Direction of W (degrees) with 
respect to the s x axis 
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Tabic 4. Plotting Programs 


No. 

Name 

Program Description 

Remarks 

t | 

PLOT 

Plots surface Isotherms 


2 

PLUV 

Plots velocities. K section 


3 

PLUW 

Plots velocities, j section 


4 

PLVW 

Plots velocities. I section 


5 

ECHKON 

Calculates equal temperature points 

Called by PLOT 

6 

CONLIN 

Draws the isotherms 

Called by ECHKON 

7 

ENDER 

Writes the values of the tempera- 
ture on the Isotherms 

Called by ECHKON 
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Jocassee 



Figure 3. Lake Keowee 






Fiqure 4. Lake Keowee (region of interest) showing inputs and outputs (for 3-D model) 
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Table s. Meteorological Data for Lake Keowee (February 27, 1979) 


Time 
(hrs from 
midnight) 

Wind Speed 
.Cm /a) 

Air Temp 
<°C) 

Dewpoint 

Temp 

<°C) 

Solar 

Radiation 

tw/m 2 .) 

Wind 

Direction 

(Degree*) 

1 

1.833 

-0.33 

-2. 78 

0.0 

15° 

2 

1.073 

-0.72 

-1.67 

o.o 

75° 

3 

2.325 

-1.61 

-1.61 

0.0 

60° 

4 

1.565 

-2. 22 

-2. 28 

0.0 

15° 

5 

2.056 

-1.83 

-1.89 

0.0 

50° 

6 

1.788 

-2.17 

-2.22 

0.0 

85° 

7 

2.012 

-2.72 

-2.70 

20.94 

85° 

8 

2. 280 

-1.67 

-2. 78 

195. 39 

60° 

9 

0.626 

0.01 

-3. 33 

369. 85 

5° 

10 

1.386 

3.06 

-2. 22 

544. 31 

75° 

11 

1.609 

5. 83 

“2. 22 

655. 31 

15° 

12 

1.788 

8. 83 

-1.39 

725. 75 

40° 

13 

3.129 

11.06 

-2. 78 

746. 68 

80° 

14 

2.593 

12.28 

-5.00 

704. 81 

70° 

15 

1.520 

13. 39 

-5. 56 

579. 20 

80° 

16 

1.207 

13.89 

-5. 56 

383. 81 

75° 

17 

1.565 

13. 83 

-5.61 

146.55 

c 

in 

in 

18 

1.609 

13.72 

-3.33 

20.94 

75° 

19 

2.056 

11.72 

-4. 44 

0.0 

30° 

20 

1.162 

9.72 

-2. 78 

0.0 

o 

in 

PM 

21 

1.772 

8. 33 

5. 28 

0.0 

o 

in 

in 

22 

2. 861 

7. 78 

5. 56 

0.0 

55° 

23 

2. 995 


5.28 

0.0 

50° 

24 

1.386 

5.28 

3.89 

0.0 

60° 
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Table 6. Inflows and Outflows to Lake Keowae 


■MBassssasaassaaas 

Time 

Feb. 27, 1978 

Oconee 

Discharge 

(m*/mln) 

Oconee 
Discharge 
Temp (°C) 

Net i ocas see 
Flow 

(C.F.S.) 

Keowee Hydro 
Flow 

(C.F.S.) 

12.00 a,.m. 

7505.3 

n 

-14395 

48 

1.00 


msm 

-18754 

48 

2.00 


18.4 

-18005 

48 

3.00 

7492. 0 

18.5 

-18713 

48 

0.00 

7491.6 

18.3 

-18698 

48 

5.00 

7494. 3 

16.3 

-18688 

48 

6.00 

7488. 2 

18.3 

-15939 

48 

7.00 

7481.8 

18.2 

3484 

3668 

8.00 

7485. 6 

18.3 

16823 

17540 

9.00 

7488. 2 

18.2 

13503 

8488 

1 8. 80 

/ / * / 

4 A 
1 Os O 

5470 

8096 

11.00 

7504. 1 

18.3 

100 

2680 

12.00 p.m. 

7503. 4 

18.4 

100 

48 

1.00 

7506.0 

18.5 

100 

48 

2.00 

7506. 4 

18.5 

100 

48 

3.00 

7503. 4 

18.4 

100 

48 

4. 00 

7501.9 

18.4 

100 

48 

5.00 

7507.5 

18.4 

100 

48 

6. 00 

7511.0 

18.4 

100 

48 

7.00 

7516.2 

18.4 

100 

48 

8.00 

751 8. 9 

18.3 

100 

48 

9.00 

7520. 4 

18,2 

100 

48 

10.00 

7516.6 

18.2 

100 

48 

11.00 

7509. 4 

18.2 

100 

48 

12.00 a.m. 

7507. 2 

18.2 

-4382 

48 







FLOW CC.F.S. x 1000) 




Figure 8. Jocassee-pumped storage station discharae 
data (February 27, 1979) 

f 





APPENDIX B 

FORTRAN SOURCE PROGRAM LISTING 
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LIST OF MAIN PROGRAM AND SUBROUTINES 


/ljA*sASAll I .THAINN FOR CRC A ICO ON 5 HAY 8(J AT I0;40|24 


1 

i 

3 

4 
* 

b 

7 

a 

9 

to 

n 

t? 
1 ! 
In 
• ( 

it 

II 

* n 

35 

•i 

• » 
• * 

’4 

35 

26 

:7 

Z‘ 

39 

* 

JJ, 

I» 

33 

34 
*S 

I* 

J7 

t; 

|5 

-c 

-t 

-.2 


•.5 


■> > 
53 
•1 
c ; 
s V 


at 
5 7 

t j 

§5 
= o 
£.1 
- ) 

£4 


t 7 

£ 9 
73 
71 
7 ? 
7 ! 
7« 
75 
7» 
77 
73 


c* «*••••«• •.****?«*»»«•#*» J* i>* 2 » •*»***• » 5 I*JS****»* 22 **»****'»**** 

c 

C ’ 

PARAHCTCR IN: | 7 ,ON:20 , kN: S , I VN s U , JMT=I 9 , kNMl : 4 
OlHCNSlOW Hi IN. JN.KN) . Vl IN , JN , kN> , VI IN. JN.KN » ,WHI 1WN« jVN.hNi. 
Cw*IIN, JN.HNI ,WRHl f UN.JkN.XNI ,P IIHN.JWnI ,0 1 IN , JN ,KN » ,C I IN.JN.KNI , 
cChLOtItVV.JWNI .XINTiIn.jAi ,,YlNI IIN.JNI.M I7«, JN.KNI.fil IN.JN.KNI, 
Chi IIN.JNI ,NXI IN, JMI.HYI IN, JN» ,HARI I n.JNI ,#1«HI I WM , JWN I , T H I I „N , JKfJ ) 
COPS* I IN.JNI .OPSY I TN.JNI 
Dl Hf NS I ON A 3 IMMI 

DIHCNSION TUN.JM.kNI ,IPI IN.JN.hM) . 101 IN , J«, KN I . NO I IN , JN , t'Nl . 
CRlNYXI IN.JN.kNI .RlNTFI 1 N , JM, kN I , UD I I N , JN , UN I , T AC I Ut I I N , «'N , NNI 


333 


161 

61 

Ubt, 


DIHCNSION TWI It-N, JWN.KNI .ROW I IVN.JWM.IOO 
DIHCNSION UAIIN, JN,XNI ,VA I IN.JN.kN) 

REAL MISS 


lN*l=Ih-| 

READ 1 , IRON 
READ I , l L N 
READ 1,MST0RC 
F ORHAl I ) 

T TOT 1 :0 * 
nT TOTrO.O 
RE AO 2 , VVIS.A6R 


READ 
READ 
READ 
HEAD 
READ 
RC AO 
READ 
READ 
READ 


i 


AI.Ah.AV.AP 

* I ,0HCGA , ARbP 


EPS » HA K 
QX.0Y.D2 
Ml | T AH t 


V 


W' c 


TAV 


2, EUL.CW.Cb 
2 , A A , CC 


2 , TLL 


A 3 Ell =V V I S 
A 3 1 2 I :v V IS 
A J I ] I :v V IS 
A 3 |4 I :v V IS 
A 3 I S I : VVI S 
B 3:VVIS 
READ 2.TAU 
READ 2, OT 
RE AO 2 ,C T TOT 
READ 2 , 1 S TOP 
FORHAT I 1 
PL2:0X»0X 
TRCF STO 

ftREF=A.6»rO»C*TO*TO 
IFTIRUN.0T.0I GO TO 3 

CALL READ JK II ,J,IN ,JN, IM, JV,Iwrt,JVN,HAR ,HRh* 

CALL INIT1 A | IN.JN.KN, INN, JVN,U,V,U,UH,D,C, 

CP,I,J,fc,IW,jV,ARBPl ’ 

CALL IMlTTr it; J.H, IN, JN.kN, IV, Jk,IWN,JWN,A,B,C,T,HO.HAIt,NNN,mF. 
CRRCF , T4|, NOV, T 0 1 

CALL HejGhTII.J.K . IN.JN.fcN.Hl.HX.MV.CCI 

CALL INLETI II ,J,k,lN,JN,KN,U,V,H,G,T , T 0 , A A , TLL ,DT ,HT TO I I 

GO TO 4 

CONTINUE 

CN AR , NRH| 1 1 )$Hj A P Ioi^6 Y ! OZ J D?I tYuX^ I »U J | W* ■W'fw’RH ' T I AHi T AV ,A, 

F0RHlTl4X.3kl6.61 
h T TOTRCT TOT* T TOT 

CALL INLCTIII, J,K,1H,JN,KN,U,V,H,C,T.T0.AA,7lL ,DI ,H) TO l> 

C Owl I I NUl 
DL 2:Ox*OX 
I STOP: 14 

IF(WTTOT.£O.O.CI GO TO 656 
DO 6 1 LL SE : 1 , 1 S TOP 
READ 2,WS,TSU,T0CW,BA0N,ISGN*,IS6NV 
R E *0 2.AN3LE 

PRINT 161 , VS, TSU.TOCV.ftAON, ANGLE, HlTQT.lSGNlC, I SONY 
FORHAT I IX, 6 1 • 12.61,21 , ,*,I5»» 

CONTINUE 
CONTINUE 
00 6 LL : 1 ,LLN 

READ 2,VS, TSU,TDCV,»AOl,,IS«MX, 1SGNTT 


48 


?9 

1C 

ai 

a 

15 
&*• 
aS 

16 
i 7 
M 
a 
«r 
71 
>2 
73 


76 

77 
?a 

75 

t:o 

131 

1.12 

IL3 

1 

13! 

l!7 
ICS 
135 
113 
111 
1 1 ; 
113 
11 - 
1 1 ! 
lie 
1 17 
11a 

fia 
121 
122 
I. '1 
12* 

125 

126 
127 

i:a 

125 
1 12 
121 
l 32 
I 33 
1 I- 
1 36 
1 36 
137 
I 33 
I 19 
IO 
1-1 
1-2 
1-3 
I-- 
1-5 
1-f 
1 

1-7 
1 - V 
* S3 


m 

153 

15- 

155 

Ira 

.57 


RE AO 2, ANCLE 

PR IN? 161 .WS.TSU.TOEW.RAOM.aNGLE ,H1T0T ,|S6nX, JSGM? 

222 TSU=TIIN*3,JN”5,|) 

TSU=TftCF«< I. 'TSUI 
IloTsTfOT *07 
TT0T1=ITGT1'DT 

TH=ITSU*T0EWl/2. . . . 

CONMENt : THE NEXT 6 LINES ARE USEo TO CALCULATE THE 

C t EQUILIBRIUM TEMPERATURE. 

C 

Fw=9.2'0.KG'VS**2 

6ETA5:0 J 35»Q.0j5*TMtD.p012*>l**2 , . 
a I &S:«. 5*0.05* 1 5U*aETAS*EW*O.OA7«fW 
TAN0 = TOE..'R4ON/KISS 
TE=« TAnB-TREFI/TRCF 

COMMENT : THE SURFACE HEAT EXCHANGE COEFFICIENT IS 

C : MON-OIHENSIONLI2EO. 

ART:KISS»0.001V1 

COnneNT : angle lOEGREES* : ANCLE a Ci.D 17«S£Z9 IftAOI ANS 1 • 

TAUxrlO. IS- 1*5 INI ANGLE t-0 . 0 I 7*. 5 329 1 - 1 SGN* 

TAUY = |Q.15**I*C0SI ANGLE *0.0 171 5 129 1 • I SGNV 
PRINT 161 .TSU.TM.FW, BE TAS.tAMn.il SI, I STOP, LL 
call error < iwn • jwn , iw, ju , ot ,Wm,uhh>i ,kn ,mrh i 

CALL WHTOPI lu.jtr, IUtJ, JWN.X*I,WH,K,MHH» 

CAlL WHAT I JU,J,«,rw.JW.IN,JN.KN, IWN, 

CALL IWTEl I.J.it.IN, JN.KW.U, V ,U.HI.HX ,MT.n«i<. * 

CAM, A V • T AUa , T AuV ,0X ■ 07 , G2 , 0 • C • 0 I , Or SX ,OPs V a AP I 
CALL C CRT N f I I , J , A ■ IN , JN , KN , A#R . U , W , X I NT . 7 1 N T , 02 ,H 1 , M AR 1 


Mn, jWN.U.UH.MARI 
HX,HY .HAR.XINT , 1 INT ,AJ. A] , 


CALL lUNin i i fiiin i irj| jNinn t*iiN.U|V iKini ifini iugi 

call rointxiiIjIk.in JN, KN ox. £v!oi, NO, Ip, eul, Hi . 
Cmar.rinTx.hx.xInti 

CALL ftOINTVI f ,J,K , IN, JN.KN.D* ,Dt .OZ.RO, AP , EUL , MI , M AR , 
CRINTY.HV.YINTI 

CALL DPSXY 1 1 .J.IN. JN, INIs JW.TWN-, JUHiOPiX ,BP S ¥ ,P , OX , 0¥ ,K ARi 
CALL FORCf II , J , I M .JU.XINT ,V INT .UHL 07 , OX « QY , H f ,HX.nV.HRu, 
CDPSX.GPSY.FH.AP.IN, JN, IWN.JVN.AINTX.AIMTV.U, V.EUL.aSr , rtA«.k«( 

CALL PRE1l(EPS,NAxIt,In,JN,P,IIN,0PSX , 0P5y , FH ,0L2 .OMEGA, 
CMRNiI.J.M,IW,JU,OX,OT.CX.!WN. JWN.ARBP l 
CALL Wf I i,J.X,IU,JU, IN. JN,KN. IUN, JWN.U ,V,0,E.H,G .OX, 07 ,D2, 

CRINTX ,RTNf Y .fUL.U.OT, aT.AP. AH.AV. A 3, HI .HX.MY.P ,MAftl 
CALL OYTOP|H,G, TAUX.TAUf.I , J,M ,02 , IN , JM.KW, HI ,MAR» 

CALL OLOuVII.J.K.IN.jK.KR.U.y.O.Cl 
CALL OLOUVI I.U.M.lN.JN.KAi.H.GiU.VI 

CALL RWHII.J.k.IW.JW, IW.JN.KN.IWW, jWW,U,V,H«,HI,0<,Dy,P2,3IK«» 

.. CALL w»ATI Jit , J.K. iw, JN, IN.JW.KN.IWH, JWW.U.UH.MARI 

35 CONTINUE 

00 20 1=1 , IN 
0020 J= 1 , JN 
WO(I.J,l>=W(I,J,ll 
2G CONTINUE 

OC 30 1=1, IN 
CO 30 J=! ,JN 
lA&I.J, 11=0.0 
30 CONTINUE 

IS CONTINUE 

CALL TEMI-II,J,M,IN f jN l KN,U,V,T, TD,0K • 

CC B , 

COY , 02 »W, 0T,TAI,TAH,TAV,Bjf,Hl l HX,HV l HAR.AHT,TRCF,1AHBI 
.CALL TCh6 2lI,J,K,TN,JN,KN,TD.CA f OY,OZ,AAR,C0,HI,AKT ,CW, TAMS. 

CHIT, MY, I, TRFF.TAV.M1.TAH, 03.071 

call olo t i i , j, k , t n , jn , xn , t , t p i 

CALL OLOTII, J,K, IN, JN.KN, 10,71 
CALL TEQKl I, J,K, IN. JN.KN. I, MARI 

CALL DENjjYl!,J.K,Jw,JW,lN,JN,Krt,IW,JWN,A,B,C,MAA.«PMrY,tW, 

CRO,ROu,RREr ,TREF l # 

DO 2000 1=9,11 
K = 1 

C PRINT 9200, ILiltX, IVl 1 , J, K) , J= 1 , JMJ I 

2000 CONTINUE 

C 9.230 FORHaTI/* L=*,X3,3X,*i: > ,13,3x,*M= , ,13/ a W-2EL0C IT Y • / I SX , gf IS.7J) 
00 NO 1 = 1. IN 
DO NO jri.JN 

„ w< I .j,l I=w0< I.J.ll 

-0 CONTINUE 

CALL INLCrilI.JiK, IN, Jw, XN ,U, V,M,G,T, 10, AA.ILE ,nT,WTroTi 
T InErT TCT1.CTT0T 
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t sa 

I 

1 cS 
lot 


US 

»•><• 

\tl 

U 7 
US 
lit 
173 
171 
» It 
1 M 
I 7*. 
I 7i 

!?$ 
I 7* 
l 7 « 
»iC 


NNN 


loco 


ts 

»2 


P«TT=»* ^ 

ir«Tl«^.|^ SO TO N***t 

TTOTUdaO 

CALL OnWHf i .J.K.IW.^V, IN, JN,.<N,UN,JUN,U,V,WM,rtI,MX,HV, 
C CALL y aS»| , l < ?jjS^IIt. JN.MNiW, V b*,N 1,H*,MT,02,*AIU 

irmSTOBf . c t I o » co to ioqo 

CALL STOftC2IU,V.UH.a,l ,3,N. lU.JW.lMiJNtAM, lUMt JWN«O.C|tU|>t« s 
CHI I MAN J MftM , A I ,AH .AV.AP ,DX,D7,67,6f , f MJ X , ! A UT , U .«R , UflH , T A J = 1*11, 

cia{,akI . ca,cu,A,J,c;cuL,TiTM Oo;»ow, re, filter, f«r£F, 10, TArtB.f oil 


* V • * a 

HTlJis^TOTfllOT 


MINT V 
POINT ? 

rORHATl. . 

fORHAIIIl, ‘TOTAL 


TJ/.MT 

>'J,TA 

II*,* 


MTTOI 
TAUX.TAU7 


TAU*=*,m.A,H*, *TAU7 = ‘,f Il.fcl 

fine Thus t ar :*,rs • — 


(M.j , •HOS'./l 
, T»£f ,R»CP,fUL,A,H,C,TE,TO» 


sstt Js?ijisiu:isv* l,,, * CB * ew ’ MT 

CALL pruvi1,j,k,!n,jn,kn,q,v,ua, V a , has I 
CALL TPOINKI I, J,M , I N , JN ,*N , T , NO , TR C f ,MA«,TACTO l» 
UT0PsI5T0P»l 
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CONTINUE 
tNO * 


SO 


ASAANASAUI.CORINT F OR CREATED ON S MAY 60 AT lOtMftsJA 


1 

? 

3 

A 

s 

6 

7 

8 
t 

10 

II 

1 ? 

13 

IN 

II 

17 

18 


THIS SUBROUT|NC ADOS INTEGRA). OF CORIOLIS COMPONENT TO HINT • 

C VINT. • 

SUBROUTINE COR INI 1 1 , J, K , IN? JN ,Hn!ab* , U?V , * I NT t Y IN T . 07 * M I < H *R I 
DIMENSION U| IN,JH,KN| , VI IN, JN.HNI .XXNTi IN,JN| . YlNIlIN, JNI .Hill 
CJNI.MAR1IN, JNI 
00 10 - 


8 

To 


jr»i jp 

i .In 


00 10 JSI'JN 
ir (MARII.JI.LT.llI GO TO V 
DO 8 KS2.HN 


uu a 

XI NT I I .JISXlKTM ,JI-ABH*HIi I, J) • I V( 1 , J, M - 1 
tint(i!ji=yint(1,ji«abr*mi( i,j(*(uii,j,k-1 
CONTINUE ' 


- ! J ♦ V I 1 , J, K I I • pZ/2 
* I »UI I ,J ,K IMM/2 


CONTINUE 

CONTINUE 

RCvySfN 

END 
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NAS A •NASA 
r 

y 

H 

t 

» 

t 

T 


i? 

I? 

13 

I* 

i* 

it 

II 

l» 

’ n 

21 

2.» 

23 

24 

25 

It 

n 

32 

31 

32 


III.OCNSTV FOR CREATED ON IS HA V 1H 43 1 1 1 3b( IM 


iHC.Fgl^Ofclfeg SMORA(J^CALCUL4ICS 3HC OENSITV fltLO FRQH 


U 

12 

10 


i* 

20 


subroutine 

CHAR «HRH, 

C 3 


0£NtTVII,J,K,IW, JW,IN,JM,kN,lUN,JWN,A,S,C, 


T.TW.RO.ROW.RREF ,TREF I 
oiflENSlON R6 (|Ni4N»KNI 
DIhEnSION ‘ROW I I WN • JUN • 


MM1 

HRHI 


3 I |Ni JNiKNl 
‘ | TUI IWN. 
JlifN • JUN I 


DIHENSION N4R < IN • JN I 
DO 10 1 = 1, IN 
00 10 3=},JN 

I r «harii!ji.cq.oi Co 30 12 
00 ll k = I ,KN 
T£H=Tll,J,KI*TRCF«fR£F 
R=A«B*T£H*C*TCh»TEH 

S8Ato2 ,:,R ' RR£r,/ * Rtr 

CONTI NUC 
CONTINUE 
00 20 (U=1,!WN 
00 20 J J = I , JUN 

IF IHRHIIi, JKl.tg.OI CO TO 22 
00 21 K=l,kN 

3£HU=TU(1N*JUiMI«TRCF*IR£F 

RUS*»a*TEHW»C*T£HU*T£H« 

ROUdU, J-,K» = IRU-RRff l/RR£r 

CONTINUE 

CONTINUE 

CONTINUE 

RETURN 

CNO 


Ul IWN r JWH.kNI 
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v. . . . -.■.•s’Z’U.'VSMtVrr* '•o — ■—r 


ASA'NASAI M.airtCSU FOR CREATED ON $ MAY SO AT lllOOtij' 


} 

i 

* 

i 

k 

I? 

II 

n 

!<■ 

I! 

r 

Jb 
13 
• " 
*> •* 

2? 

:s 

.,2* 

jr 

iZ 

:* 

:s 

:■> 

j- 

!1 

I? 

i > 

J 4 

!« 

i*. 

3 7 
» 5 

i i 

4? 

•<1 

4 ? 

4! 

44 

45 
4* 

>7 

**! 
«* V 
* r* 

ii 

* «9 

Si 

f! 

C A 

§7 
:5 
r c 

§5 

e r 

s.» 

tM 

t« 

1 l 

47 

c? 

S'? 

7" 

71 

72 
7! 

74 

7? 

7*> 

77 

7* 


IMIS SUBROUTINE CO**PUtCS UIMuUX .OIMVVV WHICH ARE USED li 
I ME POISSON C3UAT10N 7 OP PRESSURE. 


31 


37 


33 


3« 


15 


3* 


37 


SP 


39 


SUBROUIIMC DINCPUII.J.K. IH, JN ,HN , U , V . H I , OX , OY . OIMUUX .OIMUVH ,MAR I 
0 INC NS I ON 0( IR.JN.HNI.vf IN, JN,HM( .HlilN.JNI .HAAIlN.jMI 

ir ih*nu , ji. io.oi go to so 

IFtHARII.Jl.CO.il 00 10 31 

88 18 tt 

irtHARd .JI.EC. Nl GO 10 34 


GO 10 3S 
GO TO 34 
GO TO 37 

GO TO 3A 

iriHtRii ji.ro. «i go to 39 

^F I MARI I ,JI.(Q.tUIGO TO. 40 


JFlHARtl .Jl.ro. si 
IF tHAAH ,JI .CO. 61 
inMARtl , JI.EC. 71 
lFlNARtl.JI.CO.ei 
iriNAHI* 


OlHUOAStUI 1*4, J,Hl*Ui X • I . JiKMNiml, Jl-Ut 1-1, J.H| 
C«Ull-l,J.K»*HiU-1,JII/t2*0XI . 

01MUVX:iWII*l,J.K)*Vll*I«J.H)*HIIJ*l.JI-UII-I.J.Kl 
C*V|I-t .J,KI*Hl I 1-1 ,JI 1/12*6X1 


GO TO SO 
CONTINUE 


C*U(.I-1 . 


01HUVX3|Jli*l.J.KI*V(I*l iVf K I*HI|I«I,J|-U(1-1,J»HI 
C*YII-1 ,J,KI*M!I J-I.JI 1/11*6x1 
GO TO SO 


DIHUU*=IUt l *1.J,KI*U|I«I,J,K )*Ml( 1*1,01-0(1-1 ,J,K1 
]*Hl |f-l ,JI 1/12*0X1 

Ml 

• K I 

CONTI HUE 

DIMUUKCIUI 1*1 .J,KI*UII*!,J,KI*MII I«l,JI-UII-l,J,KI 
C*U( I -1 ,J,K l*Hl 1 1-1 ,J1 I/I2*fixi 
uiHUW*=tjl I* I. J.KI *VI 1*1. J,K l*H!l I *1, Jl-Ul 1-1, J.KI 
C* W I I-I . J .H I *HI t 1-1 .Jl 1/17*0X1 
GO TO SO 
CONTINUE 


01HUUXS(4*Hl(I*l,JI*UII*l,J,HI*Ull*l,J.KI-3*HHI. Jl *011 .J.K I 
:.UII , J.K | -HI II«7.JI*UI 1*2, J.K »*Ut 1*2, J.KI 1/17*0X1 
01HyvXStM*HIlI«I, JI*UI j*I,J,K 1*W 1 |«1, J,KI-3*MIII, J»*U(I,J,K I 
C*WII f J.KI-HII !*2,JI*U( l»2,J,Kl*V( I *2. J.K 11/ (2*0X1 


.•nil. 

GO TO SO 
CONTINUE 

esavYHi ;s;: y»i i nsvsi tiilUsviu y*si > ii i » rirswi -• • J '“ * 

, 01MUV*=«3*HI II .J|*UII ,J,« 1 *V | | ,J,KI-4*HI II-I,J1*UII-1 ,J,KI 
C*V(|“l.J»i* I *MI (1-2, JI4UII -2, J.K 1 • V 1 1-2, J,K I 1/12*0X1 
GO TO SO 
CONTINUE 


DlHUUXSIN*HIt I*1,JI*UII«1 .J.K l*U(I«l,j,Kl-3*HI<I.Jl*Ut I.J.KI 
C*UII,J,Hl-HII!*2,J|*UII«2 t J,H I*U| I *2,J,K 1 1/12*0X1 

01HUVX;I4*HI i I *1 |JI*UI1*I , J ,« I *V ( I •I.J.KI -3*MI II . “ ' 

C*V(|jJjHl-H 1(1*2, JI*U(I*2.J*H i a VII*/.J.K|l/(2*0XI 


JI*UIZ.J.M| 


CONTINJE 

0 IMUUY! tUII*l.J.Kl*Ut 1*1. J.KI 
C*U(I-1 , J.K l*H| 1 1-1 , J| 1/12*0X1 
01HUVX-IUI 1*1. J.Kl*ViI«l . J.K I 
C*VI|-l.J.KI*HIII-iiJII/l2«0Xl 
S3 TO SO 


K *.! • 1 * HI * 1 * 1 • " u< I" 1 • J*H » 

K l*M! ( X* t,J|-Ull-l, J.KI 


CONTINUE 

01HUUX = |i(*HIir«I, JI*UI 1*1, J.K I *UI I * 1 , J ,H I -3*HI 1 1 . JI*U(I,J,Kl 
C*UII,J,K)-HIII*2,JI*UI 1*2, j.K 1*1111*2, J.K 11/ 12*0X1 
DIHUVXSI4*HHI*1,JI*U(I.1,J,K1*VI J*l t J,KI-3*MICl, JI*U( I.J.KI 
C*V ( I , J.KI -H 1(1*2, J I *U 11*2, J.K I *V( 1 *2. J.KI 1/17*0X1 
GO TO SO 
CCNTINUE 

C 1HUUXS (U( 1*1. J.K I *lll I ♦ 1 , J.K I *HI ( I • I. Jl -U( 1-1 . J.KI 
C*U(I-I ,J.H1*M|(*«1,J| | / < 2*0X I 
OlHUV*=(Ul 1*1. J.KI*VI I *1 , J.K |*HI Jl .1,JI-U(I-1,J,K1 
C*V(I-1 ,J,Kl*Ml(I-l ,J1 >/( >*Oxi 
CO TO SB 
CONTINUE 

01MUUXS(3*HI1I,JI*UII, J.KI*U( I, J.K »-4*HI(I-l, J I *U 1 1 -I.J.KI 
“ ((2*0X1 


C*Ull-l,JcKI*MI(I-2.JI*UI I-2,J.Ai*UtI-i 
0 1 MUW* :« 3*MI 1 1. J I *U( I, U.K l*V(I ,J,K l-4( 
C*V( J-1^,KI*I , I(.-2.J»*UJ I“2.J,K ,*W(I-i 


t KI*Hl (1-2, JI«UI I -2, J.KI *U ( I -2.J.N1 1/12*0X1 

•4*HilI-l.J**UII-l,J,Kl 

• *2 . J.K 1 1 / ( 2*0X I 
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IHJS SUBROUTINE CUCUUTES OPS* »NO OPSV U^EO IN COMP- ♦ 
UMNO THE N.H.S OF POISSON'S £OU*TI0M M mICR (SRIO PO- • 


INIS 


9 

10 


|ya?SSlif 

oo id is; « IN 
Qo lo JzI.JN 
IN: I 
>*• 


fr’/HMII.Jl.LT.m 00 TO t 

OPS* 1 1 1 Jl :|P I IN, JNI-P I IN-1 , JUI «PI Iu> JI4-1 1 -P I Xu - 1 . jw-l I 1/ (2*0* I 


CONTINUE 

CONTINUE 

RCIUPM 

CNO 
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ASA^kasai 1 } .ouv v fou c «mc° # ON s^nay ao *?j uuasoa^ 

1 ^ " IN J$*SUbWOUT INC*COMPUTC5*DIMU¥Y*US£D*BV* SU# • 
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• i 

5 
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<J 

i? 

I? 

S? 

St 

H 

n 

* i 
A M 

21 

i ' 

25 
25 
:t 
2 ’ 
<8 
J 3 

i? 

j: 

\i 

it 

If 

» 

to 

iC 

41 

H 2 

■.! 

4*. 

1.5 

-i 

Hi 

49 

40 
c * 

si 

52 
s ’ 

54 

£* 

St 

57 

53 
jo 

t" 
* * 


INTE 


♦ 


31 


3? 


JI 


34 


35 


36 


37 


39 


3** 


*»0 


50 


SV!S«!IM*«Wi;U:if* 5 ; 47 li!? 3 S?/tSi!rt!!S!iSI?SSlrM 

irtHimiiJiiio.oi ao to so 

IF (NAPll ,J).C0.1 if GO 10 3| 


IF CNARIl ,j|.EO_ 
r (HASH , J I .£0,2} GO 10 32 
FIMAfll!,J.I,EC.jl go to 33 
f (MARI ilwl .fca.41 GO TO 34 

IFtMARIl JI.EO.S* ao 10 35 
i r t n Afl « i ji.eo.^i oo to 36 


PnplMUVJfHARI , 

* " “ i JN I 


IF'MRII'JI.CQ. 


GO TO 37 


!l}R 5 SM:di:lS:;t S 3 !§ 3 ! 

£ lA8v? s!uV|l5|f ! k IS vl[?J»ltK I*HI 1 1 |J4 1 1-UII tU-litfl 

C*V(1 .J“1,kI*mUI.J-IM/« 2*0Y» 

GO TO sg 
COftfY 7 NUt 

01MUVY=13*HI I -4*HI 1 1 . J-l >*U 1 1 , J-l ,h » 
CtVIl ,J»1,K)*HII J,J-2)*U| I,J-2 # «>*VM,J-2,K> I / 1 2*0 V I 

culrUjiHtivirlj » ^miV 1 JjVJ}»u}i*j* 2 V«!^v{ m?*ov» 

go to sn 

r oii 1 1 NUC 

D1HUVYSIUCI,J« | .K I * V 1 I,J.I|K>*KH1iJ« 1J-U«I,J-I.R) 
C«V<I,J-l,K|»HIl I ,J-IM/(2*0YI 
GO TO SO 
CONTINUE 

^DlHUVYS JUU , J« J j r jj« l^H 1 *M 1 1 1 »J« 1 i — 0 ( I |U~ 1 (KI 

GO ib so' * 

CONTINUE 
D1HUVYSJ J*HI II 

GO TO SO 
CONTINUE 

01NUVY=IU<I,J»l.M)*VII,J*l,K>*HIII f U*l»«UII l J-liH) 

C*V I X ,J-1 ,4 l*Ml?f ,J-1)I/I2*0Y» 

GO TO 50 
CONTINUE 

D1HUVYSI4*H1 1I,J«1 l*UII,JM.K»*VII,J«l,KI-3*HI<I.J»* 
CW(I.J,KJ*tf tl ,J,K)-HHI,J*2»*UI1,J«2,KI*V»I,J*2,KI»/I2*DYI 
GO TO 50 
C0NT1NUC 

DlHUVY=IUIIfJ*l|KI*V(I,J.l 1 K)*Hl(IfJ«ll“UlI,J-l,Kl 
C*VII,J-l,K)»HI|{,J-m/«2«0YI 
1 GO TO 50’ 

CONTINUE 

5lMUVY = 14*MlCr,JM I*UU,J*1.H»*VU.J«1.K1-3*H1»](.J)* 

GO (0 SO 
CONTINUE 

DIHUVY=«3*HIII, J)*U«I l J f K»*VII,J,Kl-4*HHI,J-ll*UII r J-l,«» 
C*W(I ,J-1 ,K|.HI (I,J-2l*U<I,J-2,KI»Vn,J-2 t H) »/t2*0VI 
CONTINUE 
RETURN 

CNO 


1 ,KMMlil,J-2|iui I,J-2,NI*V»I, J-2.KI T/I260VI 
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C**FMf5*Su25oUllNE , 'cOH?U1CS*UIUX*02UK t OlUV '• 

C*1I>M«***»*MT1«* •*#*•*•*•*»***•**•******# •***•••••••* 

' SUBGOUIIN* nVISU(I,JiK,INfJM«KNiU»V t HIiDX.OyyDlUX|D?UX f OlUViD?ll 
f HAN I 

W'fe!9* i UllN t> IN,KNI .VI f HJI INfJNl 

ir IN»1I I , JI.EQ.OI GO TO 50 


21 


32 


23 


3« 


35 


3«> 


37 


3P 




OP 


irutim ujj.to.i i 

IMHUI IiJI.C0.2l 


ifiuiSiiMii.t 

iriHiPdji.C. 

i r t n ar 1 1 idi.ro. 


? 


li 


.Jt.C0.3J 

‘ .SI 

. «bl 

iMN*nii,j>.ro.7i go to 37 


GO TO 
GO TO 
GO TO 33 
GO fo 30 
CO TO 35 
GO TO 3b 


r I ,J) .CO.BI GO TO 30 
_r UlPlI J),EQ,9> GO TO 39 
IF (HAT I 1 1 Jl .COclUlGO JO NO 


IUX STUI 


I J I » 

♦1, 


ioT = iuu,j«f>n-uii.j-i,Kn/i 2 *pyi 
2U» = IUI ! M, J,K » -jnul I ,J,HI«UII-l» J,KI l/(0x*0XI 
?UV = lun ,J*l,rt»-2*U( I , J,KI*UI I, J-l, HI l/IDV'OVI 


,j, * 11 / 12 * 0*1 


= I U 1 1 <1 ,J,K I -2*u|}j J,H I *UTI«1 , J,H I l/IOX*OXI 
:l3*yiI,J t K|tU<I,J-2.HI- l «*U|I, J«1.KII/<2*0VI 
Stun , J,K> »u 1 .sWjHl -2*U| l,J-l,H >/<0V*0VI 


HI I/T290X I 


D? 

GO TO 50 
CONTINUE 

010* = IUIT«1,J,HI-UU 

D2UXS-- 1 ‘ 

OlUV 
C2UY 

GO TO 50 
CONTINUE 

0 1 UX = I Lf I I ♦1,J,XI-U<I-J,J,KI)/I2*0«I 
0 2UX = <Ui'l • j,J f M-2*UU.J,KI«UCI-l,J,Kl|/»DXP0XI 
01UV = l*"'Ull!j‘i,H)»3*UiI .J.KI-UM ,J«2,KII/I2>.0YI 
02UV=»UI I , J»2,KI *UT T,J,KJ-3*Ull.J«l,K I |/(DV*OY} 
GO TO 50 
CONTINUE 

01UV = IUI I ,J«l,Kl-UtI,J-l ,HI )/(2*0VI 

02UV = IUI I ,J«1»K l -2*UII, J,KI«UII.J-1,KI l/TOVAOYl 

uiu*!sin*u[i*i:j-k3-3*u! j,H}=u!i»2,JiHi i/iz*ex» 

OZUJSTUI 1«2, J ,HI-2*UII»1 , J,K MUII , J t K i mOX*OX) 

GO TO SO 

CONTINUE 

Diuvrcui 1 ,4« l,H I -U( T.J-1 »K 11/12*0X1 

02UXriUlI,J,KI-2*U(T-i , J, H I «U 1 1 -2 , J,K 1 1 / <OX«DX I 
GO TO 50 
CONT INUC 

DlUysT3*UII,J.Hl+UTI,J-2,KI-N*UCI t J-l, K 1 1/12*0X1 
02UY=|UII ,J,K I «U« I ,J-2,H I -2*U 1 1 , J- I ( K I l/IOY*OVI 
01Ux-lN*UI I ♦ 1 » J ,K I -3VUT I i J|K 1 — U 1 1 < 2 , J ,K I 1/12*0X1 
02UX = IUU*2,J,H 1-2*UI I* l, J,H )«U|I f J,H J l/(OX*OXI 
GO TO 50 
CONTINUE 

DJUX=(UII*l,J,H J-U(I-1,J,K) )/l2*nXI 
0IUV = <UII,J*1,KI-UU,J-1,KI »/< 2*0 VI 
02UX=IUII*l,J,H)-2*UiI,J,K MUII-1, J.l 
D2UVS<y(ItJ*liXI-2*Ull,J,K|«U|I,J-l,l 
GO TO 50 
CONTINUE 

OlUV=IN*y(3,J«l.HI-3*UII, J.HI-Utl, J*2 , H 1 1 / 12*0 VI 
G2UV — ( V ( I . J«2 t K|«U(I,J.KI -2*Utl ,J* 1 ,H I l/(OV*OVI 
01UX = Ul*UT I*I,J,H I-J*UI I , J,H I -UTI»2,J,H 1 1/12*0X1 
C2UX SIUI 1 *2,J,H | -2*UI 1*1 , J,K | *UM , J,H I l/IDX-OXI 
CO >0 aO 
CONTINUE 

B1UXS1UI J«1,J,HI-|JII-1,J,KI 1/12*0X1 

plUY = IUf I ,J*1,HI-UI t . J-l, HI I / f 2*01 5 

DJuxriUI J*l.J.KI-2*UlI,J,MI«Ull-I '.J.Ai /IOX*OXI 

02UV~ TUI l,JM,M-l*UU ,J,H l«UI I , J-« s S I / / I OY*DV I 

GO TO 50 

CONTINUE 

01UVS|f»»yiI,J«l.HI-J*UII, J.H1-UT r,J*2,Hll/l2*0V» 
D2UY = lUI I.J*2,IM«Uf I, J,M| -2*Ull e J« l,HI l/«DV*DV| 
OtUX = <3*yi I,J.KI-«i*U| J-l, J,KI ♦U»l-2,J f K 11/12*0X1 
^2UX='Uf I ,J,Ki-2*UII-l ,J,fcl*Ut’ *2, J,H i >/,OX*OX, 

CONTINUE 


,M I I / I OX*DX I 
,H I I/TDV90VI 
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lie r vrn 
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Malawi n.ovisv ro* created ow s may so at iuiuso 
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“ THIS SWiAflUHKE- COMPUTES UIV7 ,02V7 ,DI V.X C D2VX • 
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7’ 
7- 
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7t 
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T S" 
c 


31 


32 


33 


3« 


35 


3 1 


37 


I? 


SUflNOUl IKE OVISVII ,J,K , IN , JN ,KN ,U , V ,H I • OX ,07 ,0 1 V X ,D2VX , 01 vV ,0?vy 
CHART 

DIMENSION Ul !N,JN,KNI.V|]N, JN.KNl.HlUN.JNl 'MARIIN.JNt 
iriM**<I,Jf ,CQ.OI GO I V 50 

" ■ GO JG 31 

GO 10 12 


If iNARlIiJl.CO. 

, J i • £ 0 i 


IMNARI1 i J > . 1 0 .2 ) 
IT |M*RII,JI.Cg»3l 
triH*flll,JT.r5,4T 
iriMtRt i-,ji.co.5i 

If IMAOI i.Jl.CO.AT 

, - 1 . 7 ) 


GO TO 33 
GO TO 
GO TO 
GO TO 
GO TO 


14 

35 

J 6 

37 


IT IM» R I I tJI.CO. 
ir iN«R|f ,JT.CQ. 8 > GO TO 38 
If IHARII , J I ■ CC . 9 | GO TO 39 
IFIHARII.JT.CC.IOTGO TO 40 
tllVx:n 7 ll«!,j,Kl*V(I-l,J, KIT /12 * 0 x) 
D 1 VYSIVII,J»T,KT-VII»J- 1 ,KII/I 2 *DYT 
D 2 VX=|Vtl*l,J,KT- 2 *Vir,J,K»*VII-l,J,KTT/< 0 X*DXT 
02 V 7 r t V I I , J • 1 , KT - 2 * 9 II, J,K|»V|I,J- 1 ,KT 1 / 107 * 07 T 
GO TO 50 
CONTINUE 

01 VX:|VII» 1 ,J,KI-VII » 1 ,J,KJ 1 / 12 * 0 X 1 

02 VX:|V| I . 1 ,J,KI- 2 *Vt I, J,KI« VII-I , J,K I l/IDX*OXI 

^ 1 VY S I 3 *V ,, ' • ‘ * 

> 2 VY = I VII, 

SO TO SO 
CONTINUE 

01 Vx:iVll«l,J,Kl-yiI-l,J,Kil/l 2 *OXI 


OiVY = l3twlI,J.KI-4*VII,J-l ,K 1 • V f I , J-2,K) T/ I 2*07 I 
D2VY: IVI 1 , J,K I *V| I , J-2,K I -i* V I 1 , J- 1 ,f. I J/|DY*OY| 
GO TO 50 


K I I / I D X *DX I 
K I I / I 2*0 V I 

‘ I D“*OV i 


40 


D2V* = IVIl »1,J,KI -2* VI 1, J,KT«VII-i,J,KT 1/ 
01VY-l4*VII,J.l.K|-S*VlftJ,KI”VII,J«2,KI 
02VV = I V 1 1 , J«2 ,m 1 *V! I , *J,KJ -SeVl ! , J* IsK 1 1/ 

GO TO SO 
CONTINUE 

0lVY:iVII,J*l,K»-VCI,J-l,MTI/l2*D71 

OZVtrivil ,J«1,KT -2* VI I,J,KI«VII,J-1,K I 1/107*071 

OlVX=l4«VlItl,J,KT-3*VII,J,MI-V<I«2,J,KII/t2*0X> 

D27X=IVtI’2,J,Kl-?*VlI«I,J,KI«V|l,JK|i / |0X*0X 1 

GO TO SO 

CONTINUE 

DIV7=IVI I ,J*1,K» -VI I.J-1,K1I/I2*D7I 

02V7SIVII ,J* l,KI-2tVI I, J,KI« VI I,J-1,K I 1/107*071 

01»X = I 3*VI I, J.K 1-4 *V I 1-1 , J.K MVI I -2, J.K) 1/ 12*0X1 

f)2VX = l VII , J,K I -2*VII-I ,J,K» »V|I-2, J,K1 |/I0X*OXI 

GO TO SO 

CONTINUE 

DtV7 = l3*VII,J,Kl-4*VII,J-l,KI*Vl7.,J-2,K»T/l2*071 
DiV7 = IVI 1 . J,K |*VI I , J-2.K* ~2*VII, J"1„K I I/I07*DVI 
D1VX=I4*VII«1,J I KI>3*V|1 I J,KI-V|I«2,J,KII/I2*0XI 
D2vx:rVII«2,J,Kl -2* VI I*1,J,X|«V|3 ( J.KJ|/l0X*0>iS 
GO TO 50 
CONTINUE 

01VX=IVII« 1. /U-VII-1 ,J,KI I/I2*nxi 

OIVy: IVI I, J. l.Kl-VI I, J-1,K 11/12*071 

02VX - I V II*1,J,KI-2*VII,J,KI«VII-1,J,KII/I0X*DXI 

03V7-IVI t»J«l,KI-2-VlI,J,K|«VII,J-l,KTI/l07*07l 

GO TO SO 

CONTINUE 

01Vv:i<««VII,J»l,HI-3*viI, j.ki -VII, j«2,k! 1/12*071 
02VV=IVII,J*2,KI«VII.J,H>-2*VI1,J«1,K)I/I07*07I 
D1V V =I4«VII*1, J.KI -3*V|I,J,KI-V (I* 2.J.KII/I 2*0X1 
02VX=lv< t'2,J,Kl-?*VI 1*1, J.KI •VII.J.KI I/I0X*DX1 
GO TO 50 
CONTINUE 

D1VX = IV I I* l, J.K1 -VI J-l, J.KI 1/12*0X1 
D1VV-IVII,JM,HI-VII,J-1,HII/I2*0V1 
32V X - I V I IM,J,KI-2*VII,J,KI*VII-I,J,K1I/IDX*0X) 

0 2VY : I V I I ,J« 1,K1-2*VI I,J,K |« VI l.J-I.Kl 1/107*071 

GO TO 50 

CCNTIMUC 

01V7SI4*V II,J*1 ,K l-S*VII , J.K 1 -VII , J*2,K1 1/11*071 
DZW:i VI I ,J«2,KI *VI I, J,NI -2*VI 1,J* 1,K I I / 1 07*07 1 
DIVX=I3*VI1,J,K|-4*V|I-1,J,K I«VII-2 t J,K 11/12*0X1 
02VX = < VI I , J,h I - 2* vl 1-1, J.K I « VI I -2, J,K I I/I0X*D X I 
(0 TO SO 
CONTINUE 
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ON 6 MAY 1C A? llSlfctlQ 
COMPUTES OIMW 
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Jt 


35 


26 


27 


26 


39 


HO 


30 


SuflROUMMT DWyi I.J.M , IN . JN ,HN,U. V ,MI ,0T .0 1HV V Y ,M AR» 
DIMENSION y< in.jn.kni ,vi in, jn.mni ,mi iin.jni.nariin, jni 
IFIHAPiI ,J|.EO.O? GO TO SO 
ir IMARII.JI.CO.l I GO TO Jl 
JF INAR I * - ■■ -- • “ 


_ INARl 

lrtMARU.Jt 
iriNAR<i,j),co.si 

IF IMARlI , J I . E 0 • 6 f 
IF I MAR I I , J I , E C . 7 I 
1 F|MAP| 1 ,J|.C 0 .BI 
1 FIMAPII J 1 .E 0 . 9 I 


• J i>»CO >2 I GO TO 22 
.JI.C0.2I GO TO 22 
•JI.CQ.RI GO TO 2M 


GO TO 25 
GO TO 26 
GO TO 27 
GO TO 28 

GO TO 29 

IF I MAR I I . J I . f, 0 • 1 0 I GO TO 10 

}«l,Kl*VII, J«1,K1*HIU.J«1I-VII,J-I,K»* 


OlHWTUVd, J*l,Kt*VII,J«l,KI 

C gS* f U * J-IH/I2*D»» 


CONTINUE 

01MVVir=(I*HI II,J»*VII,J,KI*VI I.J.M l.HUI, J-2I«VTI ,J-2 ,k» 
C*V|I ,J-2,K»-A|*HI 1 1 ,J-1 l*VI 1 ,J-1,K>»VII ,J-1 ,K1I/I2*0VI 
GO TO SO 
CONTINUE 

DIHVVV: IR*NI II, J» J > * V ( I, J*1,K I *V 1 1 . J* 1 ,K I -2«Hl 1 1 . JI*VIX,J,HI 
C* V 1 1 ■ J «A I -M I |] t J*2)*VII.J*?.H I • V 1 1 ,J«2,K) 1 / I 2eQV » 

GO TO SO 
CONTINUE 

DIHVVVSIVII.JM.KI • VII,J«l,K>*HIH,J*l)-VlI,J-I,*l* 
CVII.J-1.«I*MII1,J-1I|/I2*0T1 
GO TO SO 
CONTINUE 

01MVVV5IVII ,J»1 ,K>*VII,J*l,Kl *HIII ,J«l>-VU,J-l,Kl* 
CVII.J~1.MI «MI fi.J-l)l/t2«DYI 
GO TO 5U 
CCNT INbE 

D1HVVV S I 3*111 II,JI«VII. J,K»*VII,J,K)«HIII,J-2I*VII, J-2.M1 
C*W | I , J-2.K l-R*HI I I.J-I l*VI 1 . J-I.K l*VI I . J-l .K I l/|2*0Ti 

'50 TO SO' ' * 

CONTINUE 

OlHVVY:(VII,J« J,KI*VII,JM,MI*HIII ,j4ll-VII e J-l.K)* 

cvi I , j-i ,»u«hi i f , j-n i/i 2 *ori 

CO TO SO 
CCNT I HUE 

Q1HVV7=IR*HIII,JM1»VII I J*1,K»AVII,J«1,MI-3*HIII.JI*VII,J,KI 
C*V 1 1 .J.M >-HI II ,J«2I*VI I ,J«2.M I AVI I . J* 2 . K I I/I2*0V> 

GO TO 50 
CONTINUE 

QlHVVV = |VII,J«l,K»*VlI,J«l,K»*HllI,Jt JI-VIS!,J-J.K»* 

CVII ,J-1.K|*HI|1,J-1II/I2«0TI 
CO TO SO 
CONTINUE 

DlMVVYriV*HIII,J«ll*VII,J»l,KI*VII,J4l,MI-2*HIII,J)*VII,J,KI 
C«V| I . J.KI-HI I I ,J*Z|«V< l.JtJ.K I AVI I ,J«2,K I 1/12*07 I 
GO TO SO 
CONTINUE 

ClHVVv:i2*HIII,JI*VII, J,K)«VI I.J.K |.HIII,J-2I*VII,J~2 ( k> 
C.Vd jJ-ZilCl-RAMI IT .J-llfVU.J-I.M l*VII,J- 1 ,Kl 1CI2ADYI 

RE I URN* 

END 
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ACVNASMtKDIZ? SVM CflCAltO ON IS AUG 7V U IQtCNllO 


14 

s 

b 

7 

S 

9 

10 

n 

i? 

i j 

IS 
1*> 
1 * 
\* 
!-> 
'•i 

ii 
22 
:j 

29 

:« 

i) 

3 ? 

ji 

j.’ 

j? 

J9 

33 

J7 

ii 

10 


C : ; : : ; : r : : : : : i r : : : 1 1 j t : r t : i s s : : : : : : : : : i ; t 1 1 s m : : r : t : : j : J r j : : : : t : :: j i ? ! 
C THIS P90GPAM C >LCUL * T t S THC 2 DERIVATIVES 

C j : : t : : : : : j :: : i : : : r j t ; m i: s s: i : s t r s;» » 1 1 : : ii : : J t j j j i : : : ; : i r : : 1 1 j : ( 

SUBS GUI INI Dill II ,J,K , iN.JN.MN.U.V.U.HlcHX.HV .OX.OV.OZ.DIUWZ,. 
CA3,TAUA,TAUY,niVW2,6l&i } A2ui,olv2,62VZ,filA32»' , ' 

DIMENSION UllN.JN.KN I , V I l N . JN i KN > , U II N , JN , K N I , HI ( IN . JN I , 


b 2 


CMX I I N , JN I . M Yl 1 M. JN I 
DIMENSION AJINNI 
iriK.EO.ltGO 10 bl 
IE (K.EQ.hMIGO to 62 


63 


DlU2-(UIIiJfK«l|>U(I|J|K-l ) 1/(2*021 
D1 V2 : I VI l,J,MII-V(I.J.K-i>l/f2*02) 

D2U2: IUI I , J,k« M-2.*UI J ,J,K I «ul I, J.K-I I )/ID7**2l 
02 V 2 : I V I I .J.K ♦ 1 1 “7 • * V I i.J.K)*vi 1 1 j . K • 1 1 I / 1 0Z**2 1 
DlA32:(AJIK* 1*-*3CK»1 1 1/12. *021 

01UUZ = (UI t ,J,KH T*UII ,J,K* 15-UI I »J,K-1!*W( I , J.K-l > I /( 2 ,*02 I 
OIVUZ=IVU,J,H .11*1/11 , J,H« Il’VII, J,K-ll*«( I.J.K-I I 1/(2. *021 
SO TO 63 
CON I | NUC 

D1U2-HIII . J I * T AUX 
0|V2=H1(I ,JI*TAUV 

0ZUZ=2.*IU|I,J r Mll-U(I,J,K>l/(02**2l-2.*(TAUX*HI(I»J»/0ZI 
C2V7=2.*(V(I l J,K»II-V(I.J.HM/(02**2»-2.*lTAUy*Hl(f .JI/02T 
0IA32:|9.*A3(HMI-3.*A3iK|-A3(K*2ll/(2.*O2l 

0 1UW2:m .»UI I ,J,K • I I * W < I . J.H ♦ I) -3.*UI I ,J,K l*W( I, J,K|-U( I , J,K«2|*t 
Cl , J,H«2) 1/12. *021 

0lVW2:i6.*V(I.J.K«ll*W|I,J,K«ll-3.*V(l.J.KI*U(I.J.KI>V(lci.K*2)*l' 
Cl ,J i('2)I/IIi *02 I 
GO TO 63 
CONTINUE 

OlUZil 3.*UII,J.K)-4.*U(1. J,K-I»*U( I. J.K-21 1/(2. *021 
DlV2s(J.*VI I ,J.M )-9.*V(I. J.K-l I * V « I.J.K-21 1/(2. *02 1 
D2UZ=IUII, J.H-ZMUlI. J, k 1-2.*UII,J,K-I I l/IOZ**?) 

C2VZ=(V< I. J,K-2)*V(I, J, K I >2.* V IX. J.K-l I l/(02**2> 

01A32=(3.*A3(K1-N.*AJ(k-1I»A3(K-2I»/I2*DZI 

OiUW2 3l3.fUII,JK)*WII,J,KI-9.*UI I ,J,K-I»*U(I, J.K-I » 

C *U( I . J . K -2 ) ) / I 2. *02 I 

sfrtir 

END 


62 




asanas a i » 1 it i2 ?••?•***•«****#*** ***«*»*• ****,« 

2 C THIS PR05RAK CALCULATES THE HINT AND HARLOW CORRECTION T£RN AT T»* 

1 C SURF ACC # 

SUOROU*IUc’*C^ROR ( IWNiJWN* IW« JU*OT , mH.WHLOT ( RN»MRHI 
DIMENSION UHLOT I IWN.JUNl ,VHt IUN, JWN.HNI 
9 I ME NS l ON HRHMUN.JWNT 

C mHLOT IS THE Tl« DERIVATIVE OF W AT HALF GRID POINTS A? LID 
DO MOO I W ~ 1 1 1 UN 
.. 00 UOU Jw: I , JWN 

ft ir ihrh(i9,jwi*eq.oi go to sooo 

12 A/HLDt l!|| • JW»=~*Hi IV, JU, 1 1 /OT 

I! 3 ODD CONTINUE 

jiao ccNTi nui 

IS RETURN 

I*. ENO 


63 


MS»‘NA 5 *I It.FoRCC rot CREATES OH S MAy 10 AT 1 1 1 111 29 

i E**iHls"sutRoun 5 c*cOH»*o»rs ihe*r.h.s or*pois$oN$* 

.? * C EQUATION A I HALF OHIO POINTS. • 

SUiBOU t INF F OKCE IIiJiIWiJH»*IHT*TIN? tWUlOT .OX.OV.HI.HX.HV. 

CHRM, 

COPS* .OPS T.FH.AP. IN, JN. 1 UN.JUN, BINT X.BlNTT, U.V.EUL, AIR, MAH. HNI . 

DI HENS ION XINTlfN.jMI ,rIHf CIN.JNI.WHIOTUUN.JMNI .MIIIN.JnI ,HXI I 
CJMI.HYlIM.JNI.OPSXIlN. JN I, OPS VI IN, JNJ ,FH< JUN.JUNT 
01 fit MS I ON HffMlIBN.JWNt.ftlNm IN.JN.KNI ,R IN T V I IN ,UN,KNI 
0 I HE MSI OH UdN.JN.HMI , VI IN, JN.KN* ,HAB I In. JN| 

00 ion 1 : 1 , 1 k 
00 100 j:i.jn 

iriNABlI.Ji.LT.lll 00 TO 90 

OPS* II . j!:DPSXII ,ji-cul*rintxix,j,ki«vii ,j.ki«abr 
DPSVI I .Ji:OPSTI I,JI“EUL*«INTyiI,J,H»-ABB*U|I,J,K» 

9 C CONTINUE 
IOC CONTINUE 

00 10 Urn.luN 

00 10 JW:t,JUN 

1 = IW 
J=JU 

IF IKRHIlU.JUT.CC , 01 GO TO 9 

DPS*M=lOPS*l I ,JMOPi*l 1 * 1 , JHDPSX 1 1 , J «1 MOPSXI 1*1 ,J«S H/ 8.0 
DPSVH = IOPSyiI,JHOPSVII«l,Jl*OP$Vll,J«I)«OPSVtI»l,U»lT >/<(.0 

HVH:IHVI 1 (JI+HVI 1 * 1 iJ!+HV(I,J« 1 I«HV|I*I,J* 1 II/H.O 

8 XlNT=IXINTI|»I,Jl*XINTI 1 * 1 , J« 1 I*XINT II, Jl -XI NT|I,J«I I 1 / 12 * 0*1 
VINT:iVl»TII 1 J«l !*V!NTI I * I , J « 1 > -V IN T 1 1 , J I -VI N T < I « 1 • J I I / 1 2 * 0 V I 
HH:IHII 1 ,JI*HII 1 < 1 ,J)*HIII,J* 1 I*HIII* 1 ,J* 1 II/R *0 
FHIIU.JU J:| 1 . /API * 1 - 11 ./HMI*l 0 *INT*OYINTI-UHL 0 TlIU,JHT-IAP/HHI # 
CIMXH*DPS*M«HVH«.QP 5 YH| I 
9 CONTINUE 

10 CONTINUE 

RETURN 
ENO 


h 

5 

6 
T 
8 
4 

n 

ii 

i? 

i: 

i« 

is 

16 

IT 

1 s 

n 

21 

22 

2 ? 

Z*» 

25 

26 

27 

28 
24 
ir 
it 
22 
» 
u 
15 

1 ? 

tp 
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<m mi ■wmTT fog .i iM MM i 


ASA*NA«JM |l. MC IOHI row CREATED ON S hat >0 AT Ill2|t35 


k 

5 

6 
7 
« 
♦ 

10 

II 

i: 

u 

i*» 

is 

16 

17 
I” 


C IMS PROGRAM PUIS CONSTANT. DErMH FOP CC-I «Q IN IHC DATA 

£*••**•*••••*•••««••«*••••• •••»«••••• «**•*♦*•••***•♦•*••< 

c 

SUBROU I INC MCIGMTU »J.M,IN,JN.KN,h 1»H*»MV.CC1 
3 IHC NS ION HIUI..JNI ,HX<lN»JNI f MV|]N»JNI 
00 100 1=1, IN 
00 100 J=1,JN 

uiii. 1 1 -rr ' 


100 


200 

101 


HIU,JI=CC 
H*U ,J»=0.0 
1=0.0 


PRINI°loi i I I?MI 1 1 , Jl, J=l, JNI 


CONTINUE 
f OHHAT I / ’ 
PC TURN 
CNO 


1=*,I3/, * OCPTM*/5X,9ElN.7» 


.ASA*MAt«f II.1NIMA TOR CRCATCO ON 1* HAY 7* <T UlSllOl 


1 

2 
» 

i 

i 

i 

I 

« 

I*' 

It 

ft 

IN 

I* 

i? 

is 

15 

:i 

22 

23 

2 ** 

25 

>«. 

5 T 


r**.»^'*t»v>CSSS?*««2»!»*««i**«*«*««**««»***«*»*««««»*iSi*i*»**«*«***** 

SUBPOUT IMT INIT1/M IN, JN, KN, INN, JNM ,U , V , M , WH.O .C .P IN tM , I , J.K , IV , 
CARBP I 

0 INC NS I ON UUN.JN.KNT , VtiN.JN.AN) ,MTIN, jN.KNT.WMI IWN, JUN.KNT, 

C INITIAL CONDITIONS ON U aNQ V 

10 00 100 1 = 1, IN 

11 00 [00 JZI.JN 

oo ioa k=i,kn 

UII.J.KISO 
VII,J.KI=0 
W(I,J.k»:o 
oti . j.H i :o.g 
E 1 1 * J.K I » C »0 
100 COnTINUC 

C INITIAL j-OM^ITJONS^ON UH AND PM 

00 loo Jsj 5 1 s JUN 
PINTHT IV. JWI=ARBP 
00 200 K: 1 • KN 
WH(IW.jW,hI:0 
20d continue 

RETURN 

END 
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NA5*»NASMl I.1MUI FOR CREATED ON 15 HAT 1H AT 11U5I25 


I 

* 

i 

« 

5 

5 

1 

• 

9 

I? 
II 
I 2 

\! 

15 

15 

IT 

19 

13 

n 

a 

25 

25 

27 

13 

j? 

ii 

tt 


C THIS PROGRAM INITIALIZES TEHP AHO OINSITV # ( 

C *** , *?unOu*INC*INI»IM?t**KtlNljN*SN*J5tjS.IHN,JMN,A,a,C,T.NO, 

C1W.ROU.TOI 

DIMENSION T I IN .JN.MNI.RO I IN.JN.MNI .VblfWN. JUN.KNI.ROWI IWN. JWN» 
01 HE NS 1 ON MARIlN.JNI.MRMIUN, JUNI 
lODSiTO-TREr l/1R<F 
** r:a.B«IO»C«TO*TO 
R0D=4R-R«EF l/RRtr 
00 10 i:| .IN 
DO 10 J=l,JN 

IF IHARI I , JI.CO.OI GO TO 12 
DO 11 H=|,KN 
TII.J.MlSIOO 
ROlljJ.KCROC 

11 CONTINUE 

12 CONTINUE 

10 CONTINUE 

DO 20 lWsl.IUN 
DO 20 JM= l , JMN 

IF IMRHl3WijWI.EO.OI GO TO 22 
00 21 H - 1 ,KN 
TMIIU.JM.KKTOD 

21 CON{iNul , *' ,X,::W00 

22 CONTINUE 

20 CONTINUE 

RETURN 

ENO 
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IM.CM CL T CREATED ON 5 HAT 10 AT 11)36)12 


I 

S 

« 

5 

) 

9 

« 

1C 

n 

i* 

IS 

\) 

M 

\] 

21 

22 

23 

26 

2 * 

26 

n 

2 ° 

3P 

II 

32 

33 

!? 

36 

37 
36 
3? 

63 

21 

6! 

6 *. 

6 * 

64 
67 

Hi 

69 
c 1 

h 

5! 

5 6 

55 

56 

57 
51 
«7 
fc'’ 

6 1 
62 

63 
66 

65 

64 
67 
61 
69, 

70 

71 

72 

73 
76 
75 

7? 

f* 


C* THIS SUBROUTINE F OR 1MLCT ANO OUTLETS FOR 

C DOHA IN 

SUBROU f INC jNLCTH!tJ(K«INtJN|HN»U«V|HiG*V»’TD»AA*TLLiOt*HTTOM 
DIMENSION HI IN,JN,KN> ,G< IN, UN, INI , Ul I N , JN.KNI 
0 INC NS ION VI In|jN,KNI,TI1N,JN,XN) , TOI IN,JN,KNI 


TSOT ;HT TOT 
I SMI TO! 


t SO T 1 ' 


COMMENT 
C 
C 

c 
c 


PROGRAM HAS 8CCN MRUTEN5PCC1FJCALLVFQR 
KCOUCE. USERS MUST CHANGE TO iutl SPECIFIC 
SUE. UATCH OUT FOR COMMENTS TO START OR STOP 
CHANGC. 


THIS 

LAKE 


*NHl 

c 

comment 
c 




l START CHANGE. 


j »l »3> = < 

G I 1 1 , I , J ) = AA 
HlU l, J|=0.0 
T 1 1 1 • 1 1 3 ) = TLL 
TO 1 1 1 » I • 3 T = ILL 
10 CONTINUE 

5F=0. 00322579. 

PI 

A J 0 £ 

B JOCSE: 

IF I T SOT. GC.O.O.ANO.TSDT.LC.l.O)SV = SF*l -16.395-1 18. 75-16 
C. 3951*1 


lannin 3 

C SC = 1 1 • 6 52S 


lMTSol.GE»l.n«ANO.TSDT.LE«5.0lSV = SF*|-18.756) 


IF I TS0T.GE.5.0.AN0.TSDT.LE.B.0 
CIS. TS6 1 / s. )*ITSDT-5.QI *18.756 I 
IF I TSOT., GE .8.0. A NO. TSOT.LE.il. 
C0.1l/3«|6|fS0r-6.OI-16.B23) 


8.0fSV=Sr*l 1116.823* 
0>SV = -SF«H 1 16.823- 


IFI TSOT. GE .11.0. ANO. TSOT .LI. 23.01 SV=SF*0.1 
IPIlfDT.GE .2 3.©. AN0 s ,TS0T.LE.26.0l SV=-SF*( | 1 6. 5 *0.11 1 


COMMENT 
C 


C*TS0T-0.l ) 

STOP CHANGE 


00 2n K= 1 iKNMl 
00 20 J=8. JNM1 
Ull . J.K iriv 
Hi 1 1 J.H l=SV 

Vll,J,K)=g.g 




TI2.J.K) 


20 


T0I|jJ|JK) = T o I2,J,(|) 


C 

COMMENT 
C 


CON 
00 30 K 


l.KNMl 
START CHANGE. 


GIT, 

C 

COMMENT : 

C 


UI7,5,K) = 2.*Ulf ,S,1>- 
H|7,5,1 )=2.*UI8,5,K). 
V I 7 1 5 .H 1=0. Q 
' -,5,11=0.0 


'U I 9 | S . 1 ) 
'019,5,1 ) 


STOP CHANGE. 


30 


65 


«0 


CONTINUE 

51=0.05789526 

IF I TSOT. GE .0.0.ANC.TS0T.LE,6.Q|SVI=0.068+SX 
IF I T5DT.GE. 6. O.AND. TSOT. LE.B.OISV 1=51*11 117. 56 -0.068 1/2.1* 
CITSDT-6. 01*0. 0681 

IF I TS0T.6E.8.0. AND. TSOT. LC. 12.01 SVl=Sl*< 1 1.068-17.56 1/6. 
CI*ITS0I-8. OMIT. S6) 

IF I TSOT. Gl. 12.0. ANO. TSOT .L£ .26 .01 S VI =51*0.068 

00 60 1=1,3 

Ult3,7,K)=SVl 

MI13.7.M t=SVl 

TI13.7, 11 = 7112,7,1) 

T 113,7, 6) = Tl 12,7,6) 
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iS 

u 

n 


i»# ,m u " ,w 


tm 


..V- 


69 


UA‘NAS/»tl 

? 

<4 

* 

A 

T 

A 

* 

1C 

ii 

u 

IS 

ir 

IS 

1*5 

Vi 

22 

21 

24 

2*! 

26 

IT 

zr 

in 

ir 

31 

32 

33 
• 

is 

i*- 

17 

3b 

34 

nr 

4 1 
42 
4 I 
44 
4t 
4* 

U 7 

41 

44 

53 

51 
• 52 

53 

54 

55 
55 
57 
I! 

5" 

or 

cl 

52 

53 

64 

65 

66 
c»7 

6* 


i.imc r on created on s hay ao *r iiiUtiT 

c »*••***•*•••.•••♦••**•****•♦*♦* n* /»«*♦**♦♦•♦*•••*•*« 


1U.U 


1 6 J 

:cj 


OIMENS 
OIHCNS 
OIHENS 
OIHCNS 
OIHCNS1 
00 200 . 

DO 200 JSl.JN 

1FINA4J ?,j|.CO.OI GO TO 20U 


,HIINiJNtKNI,HARtIN,JN),HIUN,J 


SUBROUTINE INTEII,«M,IN,JN,XN,U,V,U f HI.MX,HY,HAR,XINMINTiAJ 
C.AI.AH.AV. TAUX.TAUY 

c i?s ssiSii'SiJs *j4j«!:Rrc:K!M.M 

ON HXIlNjJN| f HTIIN'J*) 

ON 0^5 1 JK’uNlInPSVIlNi JN* 

ONjOIIN»iJN»HN| ,t| |N t JN,MN| 


A r I r* » | a iw p * b4 

tinMI.jMo.q 
xiujll.jiro.O 

ALL OWWYII l J»K»lN»jN»KNtU»V#Hi*I)T tDlHVVV |HAR 1 _ 


UVVTM(J l n»ji4 l jrfinr<fU > Vini«uTtuinv«r»n4Ni 

0V1S0II,J,H,IN,JN,KN#U, V ,111 , OX , D V , 0 IUX , 02UX «Q l UV ♦ D2UV , HAgl 
OW I SY 1 1 ! J»K } IN# JN.KN*U, V .HI (OX,OY i^l VX.D2WX |0 IVY ,02VY »HAR I 
D12Z(l,J«,H,lN.JN,HN.U.V.U|Hl.HX,HY,nx l 0V«0Z ( 0lUU2.A3, 


ioao 


101C 


call 

CALL 
CALL 

CALL 

cHtt.?i8?!{lttS:4T«K«lilWIH:S}K:HIfi5 

ir IK.CO.il GO TO 1000 

IF IH.CQ.KNI GO TO 1010 

XSUN:i*l4|0]HUUX«01HUYV4HII I . JI*0IUH2I 
C-AH*l02UX*HlU t J>*02UV*HIIl,Jli 
C-AH4I01UX4HKII,J|«01UY*HYI I.JI I 
C-AW4I l.O/HIli.j) )*IA3IM*D2UZ«01A3Z*0IU2I 1*02 
YSUM5I AX4|01MUWX«01HVVY*HI I I . JI40IWU2I 
C-AH4i02VM*HIIIiJM0?VY*HII I.Jll 
C-AH»(Q1VX*HXI I ,J1 ♦01YY6HYI 1 ,JI I 
c >AV*n.0^HM!|Jn4««3IKl*02V2«DlA3Z«DlVZI l*OZ 

CONTINUE ° 

X SUNS | A I * 10 JHUUX*D 1HUV V*HI II. J )40|UWZ I 
C-AH*IQ2UX*HI 1 1 , J H 02UV *H 1 1 1 ,JI I 
C-AH*l01yx*MXU,J>i01UY*HYH,JI I 

C-AV4ei.fi/HMI.JI |4|A3IM*02U2*01A3260lUZ||40Z/2.0 

crr*!iS};i8ivvi!5;?^ tf Y:«!{!^r o,vwn 

C-AH4I01VX4HXII . J I *0 1 V Y6HV 1 1 « J I I 

C-AV6I1.0/MIII .JII4|A3lKI*0ZVZ*DlA3Z6DIVZfl40Z/2.0 
01uT=lUII,J.KI-0t !,J,K I I/O I 
01iiT = (ve! 1 j,KI-CII,J,K(l/0T 
Q £2.0/132 

&PSXIX,JI=I1./AP)*II./H1II,JI I6I-XSUH4Q-HM I, JlcOlUTl 
tlP5YU f J) = ll./AP)4ll./Hl 1 1, Jl I*e-YSUH*0“HIII,J)401VTI 
GO TO ltoo 
CONTINUE 

XSUH=IA14I01HUUX«01HUVV*HIII,JI401UW2I 
C-AH6 1 02UX*HI I I , J I *0?U V *HI 1 1 , J I I 
C“AH4 (DlUX*HXlI» JI«DlUY*HYjI , Jl I 

C-AV4I l.fi/HI II, j| } 61 A3IK 1*0202 *0IA32*01UZ| I40Z/2.0 
YSUHS I AI* IDIMUVX»01HVVVHI C I , JI*D1VUZ I 
C-AH4I02VX4HII1 ,J I 402VV4HI II , Jl I 
C-AH*I01VX»H* I I, JI«OIVV*HV< 1. Jl I 

C-Av*l 1.0/HIIl , Jl 1* I A3 I K 1*02 VZ«DIA3Z*01VZ 11*02/2.0 
CONTINUE 

XIMTII,J|:XSUH«X1NTII ( JI 

VINT I I ,JISV5UH*VINTII , Jl 

CONTINUE 

CONY INUE 

RETURN 

C NO 


70 


If 


Jl"f SVM CRCMCO ON 19 NOV 7? AT OytS7tOO 


4 

1 

4 

5 

6 
T 

• 

* 

10 

11 

1 ? 

i 

16 

IT 

18 

I " 

2 C 

21 

27 

23 

24 
2 ! 
26 

H 

sc 

11 

12 
3? 

34 

14 

it 

37 

38 
34 
42 

41 

42 

43 

44 

45 

46 

47 

48 
44 

if 

52 

< i 

54 

55 
5b 
57 
59 
59 
6 " 
fci 
62 
67 

64 

65 
65 

67 

68 

69 

70 

73 

74 

75 
75 
77 
79 


V . : 


.UO1QI.O.70 


l:W;? 83 fi? 8 i!!S; , -° 

0.0526,0.0526,0.25 


_ ^0.0526,0.2! 
1.0*0.0092,301.56 
1. 150^18,-0. 000019,-0, 0000044 


. 742.22,0.0,0.0 
0.2? 717,1.0' 


S:8:u;8;u:3:I: i t^i? i ! :f ;!i ,11*11 .ii , 


Q , a ,’ 2 ; n ; ij , 11 , 11 , 11 , 


liillil iliil iliitiWift 


u,u , u , u «<, jj , 1 1 , 1 1 , j * , 1 1 , 11 , 1 1 , 1 1 » J I , M » * i , i , o » u , 
u.o,o,o,9,i,n,!i,ij,J,h, ji,ll,i,,ii,irq,g,o 
O.CiO.g.g,?,!!,! ,1{, , 1, *11, • ,1 ,6, 1.3,5 

o : o ;?; } , i : , , , i !, j }. H .} i .! f . M . M , 6 , i , t !, 6 ,iD 


I'M 

2 : 8 : 8 : 8 ;S:S: 4 U j ii i;S* 1 ; Ulh h:i:S: 8 :S t0 

o!u!o!oJo!ol4ji6!l!o?io!i6?lol?o?lo!?g?ia,iQ,io,3 

0,0, 0,0, 0,0, 2,9, 9*9. 9*9, 9. 9,9*9, 8,8*7 

8 :S: 8 :S:i:i^ s’sw&MfsfjiiWi 0 
8 : 8 : 8 ;g:S: : : : 

0,0, 4, IQ, 10,9, 9, 9, 9, 9, 9, 9, 9,9, 9, 8,9, 9, 7 

n.D.^,9 *9«9»9»? »9 • 9 • 

4;io;9i9,9,9;9,9,9,9.9 

fe ' S ' n ' fi ’ n ' f , s * f *2 

3, 0,6, 8, 9, 9 


0 , 0 , 0 , 0 , 0 , 0 , 0 , 2, 9 

0 , 0 , 0 , 0 , 0 , 0 , 0 , 2,1 
a, 0,3, 0,0, 0,0, 6,7 
1.633.-0.33.-2.78 
JS,0 

1.073,-0.72,-1.67 


75.0 ’ 
7.325,- 


1.61,-1.61 


2 

60.0 

1.565,-2.22,-2.28 


15.0 

1.056,-1. 


83,-1.89 
2.17,-2.22 


2 

50.0 
1.788, 

85.0 

2.012,-2.72,-2.78 

85.0 

2.260,-1.67,-2.78 

□? 4 § 4 , 0 . 0 t , 

5 #U 


9,9 .9.9.1 ■ U . a . b . 7 .0 
9,9, 1,0, 0,0, 0,0 


9,8,8,9,9,10, 3, 0*0,0 

1.0. 0.2. 9.9.1 .0.0.0 

7. 0. 0.6. 9. 9. 1.0. 0.0 
0 , 0 , 0 , 0 , 2 , 9 , 9 , 5 , 0,0 
0 , 0 , 0 , 0 , 6 , 8 , 8 , 7 , 0,0 
0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 . 0,0 
0 . 0 , - 1 ,, 1 


0 . 0 , - 1,1 
0 . 0 , - 1,1 
0 . 0 , - 1,1 
0 , 0 , - 1,1 
0 . 0 , - 1,1 
20.94,1,1 
195.39,-* ,-l 
3.33,369.85,1,-1 

1.386 .3. 06, -2. 22, 544. 31, -1,-1 

1 .609.5.83, — 2.22 ,655.96, 1,-1 

15.0 

1.788.8.83, -1*39,725.75,1,1 

50.0 

3.129.11.06, -2.78,746.68,1,1 

80.0 

2.593,12.28,-5.0,704.81,1,1 
1.520,13.39,-5.56,579.20.-1,1 
1^207, 13. 89, -5. 56, 38 3. 81, 1,1 

1.565.13.83, -5.61 ,146.55, 1,-1 

1.609,13.72,-3.33,20.94,-1,1 
15. U 


71 


I 

Si I?U2,«. 72,-2.711, 0,0, *1,-1 ' S 

SI j?m,i.j],s.2«'0.0ii,-i 

*9 2?i£l,7.7«,5.56,0.0,l,-l 

St |:«°5,7.00.9.2S,0.0 f l,<l 

*Q *0*0 I 




72 



73 



N'ASA.NAIAI ll .OLBMV roo CRE HTfTft ON 1* **»V 7M AT 151*9:33 


I 

j 

j 

.. 5 

5 
9 
T 

6 
9 

13 

11 

12 
11 


C IN ORDER TO RETAIN VALUES OE U AND V AT ONE TINE STEP LAC 

C*««««****«**«#*«<*TM«'*9**« + «**t**««*e* + *«9»9«#*9*+*«**#«t*+« + *4« + + *M 

SUBROUTINE OLOUV H , J.M , IN , JN, KN.U , V.n ,[ 1 

01 HENS ION U UN, UN, AN J , VI IN, JN.HNI ,DIIN, JN , KN 1 ,£ UN, JN ,HN » 


9 31 


DO B31 «:l,«N 
00 A3) 1 5 1 J i'N 

00 eSi usiJjN 

01 i ,j,a i run ,j,ki 

END 
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l I.PRCIL FOR CREATED ON 5 HAY ID AT > 11 * 312 * 


I 

; 

■ i 

5 

6 
T 
8 
9 

1C 

11 

I ] 

It 

15 

16 
17 
IS 

19 
2 C 
21 
22 
2* 
29 

25 

26 
27 
2 E 
29 

20 
11 
12 
11 
It 

15 

16 
17 
13 
29 

to 

tl 
M? 
t ; 
it 
t5 
1 6 
tT 
tl) 
t< 

SC 

51 

52 

53 
St 

55 

56 

57 
59 
59 
6 P 
61 
62 
62 
69 
65 
06 

67 

68 

69 

70 

71 

72 

73 
7t 
7E 

76 

77 


C THTS SUBROUTINE COMPUTES pACSSUaE FOR FAR FIELD 

c 

SUBROUTINE PPE IL ICPS.HAXl T , IN» JN,P, I TN.OPSX .DPSY .FM.DL2 .OHEGA, 
CMRH. I. J.K, lW.JW.DX.DV.CX. INN, JUN.ARflpI 
DIMENSION P ( 1 UN , JWN I ,F M 1 1 UN , JUN I , DPSX I IN , JN > ,0PSV I IN, JNI 
DIMENSION HRH ) 1 UN , JWN I 
1 7 N = 0 

1 r xsO.O 
DOSARBP 
I TN- ITN« 1 
00 1U IUO- 1 , IWN 

00 In JW- I , JWN 
I W : I IHN • I l-IWO 

1 = IW 

'iMPHIIM. JUI. £0.01 00 TO 57 
<HRHI!MvJWI.EC.l> CO TO II 
IHRH<IU,JWI.EQ.2> 

IHRHilb.JWI.EC.3l 
IHRHIIM.JWI .EO.tl 
IMRIMIM, JWI .EQ.5I 


51 

57 

10 


20 


l'r 

IF 

IF 

IF 

IF 

IF 

IF 

IF 


IMRHI IM.JWI .CO. 61 
IMRHIIW.JUI.EQ.Tl 
IHRHIIM.JWI. CO. 8) 


GO 

GO 

GO 

CO 

GO 

GO 

GO 

GO 

GO 


TO 

TO 

TO 

TO 

TO 

TO 

TO 


12 

13 

1* 

1* 

16 

17 

IB 


if imrhI!u;jui.eo.ioi GO TO 10 

PN=.25*IPTIW-1 , JM) *PI IW* 1 i JU 1 


. . _ . «PIIW,JU-1 l«P|IU t JU«ll-DL2*FM<IU,JI 

GO TO 50 
CONTINUE 

PN=. 25* IP I IU-1, JUI *P<IU«1, JUMP 1 1U,JU-II«P I IU,JU)«IOPSV< I, JMI 
C* DPS VI 1*1 ,J«ll l«0V/2.>DL2«FHI IM.JWlI 
GO TO SO 
CONTINUE 

PN=.2S*!S| !U»!jJuj .oi JH*1- jin 4p? IU, JU*iMP(IU, JUI 
C-lDPSVM, Jl*Dp£v<IM.J>l«6v/2-DLZ«FHl IW.JUII 
60 TO 50 
CONTINUE 

PN=.25*|PIIU«l,JUI»P|IU,JU-ll«P|IU,JU)-inPSXir,J)*DP5Xll,J»lll 
CtOX/2*PITM,JU|«IOPSYII,J«ll*OPSVII«l,J*lT)*OV/2-OL29FH< IM.JWI I 
GO TO 50 
CONTINUE 

PN = .259IPIIU«1,JUI«PI IW.JW«1)*PIIU,JU|-|0PSXII, JMDP5XII,J«1| |t 
CDX/2«PIIW, JW l-IDP5VII,JI’OPSYM*l»JII*DY/2“OL2*FHI IW.JUII 
GO TO 50 
CONTINUE 
PNiARBP 
GO TO SO 
CONTINUE 

PN=.25«IPIIU,JU«1I»PIIW-1,JUI*PI1U,JMI«<0PSXII*1,J«1I*0PSXII»1,J 
C*OX/2*PIIW,JUI-lOPSt(I t JMOPSVIlM.J» I *0V/2-DL2*FHI IW.JUII 
GO TO 50 
CONTINUE 

PN=.25«|P<IU>1,JU)*P(IU,JU-1I«PIIU,JWI*<0PSX(T*1, J|«0PSXII«I, J*l 
C*OX/2*P|IW,JU|«|OPSVII I J»1I ♦DPSV 1 1« 1 . J* 1 1 1 *DY/2"DL2*FH I IW.JUII 
GO TO 50 
CONTINUE 

PN=.25«IPIlU.JU«ll*PllH-l,JWI»PIIW,JU-l)* 

CPIIU, JUI* I0PSXII«1,JI*DPSXI1« 1, J*lll*0X/2- 
C0L2*FHI IW.JUII 
GO TO 50 
CONTINUE 

PNS.25*lPIIU«l , JMI*PIIW, JU-1 l«pl IU,JU«1I«00-0L29FHIIU. JUM 
CONTINUE 

PNEU = OHEGA*PN*M-OHEGA|«P|IW,JUI 
IFIAB5IPNCUI.LT. I 1 0. *•1-16111 GO TO 51 
OIFF = »BS I IPNCU-Pl IU.JU I I /PNCUI 
IF IOIFr.LT.Exl GO TO 51 
E X SDIFF 


PIIU.JUI 

CONTINUE 

CONTINUE 


= PNEU 


IFIEX.LT.EPS! GO 

IFIITN.LT. MAKIII 

CONTINUE 

RCTURN 

END 


TO 

GO 


?g 


7S 


AS*»NASAM I.PjnEj^FOB^CRCMfO^ON^tJ Hit 7H A| 15 j«V»N2 ^ 

? C ■ 'IMIS PROGRAM PRINTS UUI f HE VALUES OF NUMBER OF. ITERATIONS ANOF 

l 4 ^RESIOUAL^ERROR IN SOLVING POISSON 

5 SUBROUTINE PRI 1 T K < MN,EX»****** m *********** ** M ••••••** * 

4 _ PRINT SSOO.ITN.IX 

T 5500 FORMTIF* 'lTNr»,l«,5* t ' EK5* f E15.|» 

” ENO 


76 


•NASA • 1 1 

? » « 

3 

« 

5 

6 
T 
9 

• 


PfiUV STN CREATED ON S HIT 10 IT 11:9I:U 


this subroutine minis thc value 

OHIO POINTS. 


or U INO V IT ill HUN 


n 


15 


13 


1* 

C 

19 

C 

It 

9100 

IT 

19 


1« 


2C 


21 

I 9 C 

22 

153 

23 


29 


25 


21 


2T 

191 

29 

IS1 

29 

IDS 

32 

107 

31 

10b 

37 
• » 



SUBROUTINE PRUVI I , J.K, IN, JN,KN,U,V,U* , VI, MIR I 
01 HE NS ION UIIN.JN.NNI , V( IN, JN.KNI , MARI IN, JN», 
CUMIN, JN,KN| ,yM IN.JN.KNI 
DO 9100 K=1,KN 
00 9100 J=I,JN 
00 9100 1 = 1, IN 
UMI f J f MI=uit f J.N|*SU. 

VAIl,J,Kt=VtI.J.K|»30. 

IF I MAR I I, J I. EO. 0 I U* = 1000000. 00 
IMHlR,I d J ,.E0. 0 , 91 = , 000000.00 

00 150 K=t.NN 
NRITEI6. lOSlH 
DO 190 1 = 1. IN 

WRITE ( b . 1 P* l (UM I f J.MI , J= 1 , JN l 
CONTINUE 
CONTINUE 
00 1S1 N=1.KN 
WR IICI6.107IK 
00 I 9 1 1 = 1, IN 

WRITE lb l IU6liVIII,J,HI,J=I|JNi 
CONTINUE 
CONTINUE 

FORMAT, 'I'.'U-WELOCITT FOR 9 = *15» 

FORMAT I *1 • • • V-VELOC 1 1 Y FOR HS'IS, 

FORMAT I// ,}2F6.2I 
RETURN 
END 


AS«*Msti i i.rcaoi ran created on i mo rt at izuatj* 


i 

? 

i 

« 

'5 

6 

T 

8 

is 

IJ 

i l 

u 

II 

•I* 

2 ? 

21 

22 

23 
2 * 
;s 

24 

ii 
2 « 
IS 
II 

i'i 


c this proorah rcaos jape for oata i tor. the variable ocnsitv cas 

SUBROUTINE RE A02 I U , V , UH ,P IN TM, I , J , N , IN, JU , I N ,JN,KN , I UN . JuN ,0,C , 
CHXfMYfHI 1 MARfHRH 1 Al 1 AM,AV,AP 5 OX # Ot t 6i {OT|IAUI 1 TAU3 1 U|UR 1 HRHi 

‘ iT0»TARB 


jj.JdU.U.J.', 

0 I HENS I ON HVI1N.JNI ,HVI IN,jNI r HZIIN,JNI .NARI IN,JNl tMRHIIUN, JUNI, 
Cltf IN. JN.KNI .NRI In, JN.NNI ,URH| 'UN, JUN.KNI 
oIhcnsiSn Tl IN,JN,R n! .ROIIn,, »,KNI , TNI I UN, JUN, AN I ,ROU<iUN,jUN,AP 
RE AO (7,CN0:II I HUI I , J ,H I ,‘i = I ,KN I , j: 1 , JN > ,1 sj , INI , 

C M l V I I i J * I ; , A - 1 , A N I , JS I , JN I • I : 1 , IN i , 
Cin0U,J f Kl,x:i,»(N»,j:i,JN»,i:l,IN», 

C( IIC( I | J«HI,K=l.MNI,j:i f JN|,E=1.INI, 

Cl I INHlIu,JU,N I ,K = | ,KNI ,JU=1,JWNI ,Itf=l,IUNI, • 

C<| IUtI,J,KI.K:I } «N>«J=i, JNI ,lr mInm 
C lfl UR tf.J,AI,M=l,NN|,j:|,JN 1,1:1, INI, 
CIIIURHI!U,JU,Nl,ft=l,AN> ( JU=l, JWNt. !U:|, IWNI, 

Cl (PINTHIIU.JUI , ju : 1 , jun i .lull ,]UNl 

C,|IHIII,JI. jr l.JNI.isl.INI, II HA I If JI.J=t.JNi 1 i:!,INI.<IHrlI,JI..' 
Ci, JN I|I-I ( INI,IIHARII,JI|JA|, JN I — laZNI ,1 1 HRH I I U , JU I , JUS I , JUN I • 
C I u si , I uN i , 1 1 if II ,J,hi,hsi , hN i ,j:i,jnI , 1 = 1 , INI, 
ci iiroi i,j,ki,ksi,nni, j:i, jni , 1 = 1 , Ini. 

Cl IITUIIU«JU,AI,RSl,HN>,JN:i,JUH),IUsI,lUNl ( 

Cl I IROU I I W, JU.,K t .HSl.MNI , JUS I , JUN I , INS I, I UN I , 

CTAl.TAH.TAV.AKI.CfMCU.A.B.C.CUL.T, I u , r6 ,ROU , TC .RACF , TREE , TO , T AH* 

cai,ah,av,ap, ox, ov,07,ot, Taut, taut, ttot 

CONTINUE 

RETURN 

ENO 


I 


ASfNASAI I I.REA03K CL T CREATED CN i HU BO AT Hi 50, 20 

c this subroutine reads amo prints the mar l mar matrices 
c 

SUBROUTINE RCA03H I £.J, IN. JN, I U , JU , TUN , jMN, «AR , mRHI 
01 MENS ION MARI IN,JNi,MRHtlHN;JWNB 
DO 300 1=1, IN 


?-r 


to 

11 

300 

i: 

13 

!*• 

«co 

15 

2 

ifi 

IT 

1 ■ 

3 


READ 2.IMARH.JI ,J = 1,JN) 


00 <400 t US I • I HA, 

RE AO 2 j (MRH I IU , Jw I ,JU = 1,JUNI 
PRINT^J^IWtlMRHI IW,JUI ,JU=1,JUNT 

rOHMATO 

FORMATl/15*, SaiH 

RETURN 

END 


79 








< 




j 

£ 

< 

i 

7 

a 

a 

1C 

M 

1? 

I.' 

1 N 
IS 
IN 

n 

i* 

:i 

2 2 
23 

£ 

27 

a 
2<s 
10 
Jl 
12 
12 
lo 
I! 
35 
37 
Si 
3* 
0 0 
01 
02 
<*3 
*o 

<4« 

07 

09 

<4 5 
5C 

51 

52 

K V 

SO 

55 

56 
« 7 
St 
5® 
6? 
fcl 
62 
6 3 
60 
65 
6‘. 
67 
t* 

69 

70 

71 


70 

75 

It 

77 

7* 


C 1HIS SUBSOU TIME COHPUICS AP IN INC POLSSONS COUATlON 
C «•••*••«• •••*••••**** *•**•*•• •»«••*••••••••••••**•*•••*•*• 

SUOROUTINl R0INrKII,J,K,lN,JN,*N,0X,0T,02,R0,AP,CUL,Ml f 
CHAft RIN7X MX XiNTt 

OlN^SlON'niNJkUN.JN^Nl .(JO I IN,jN,KNI , X I N T i IN , JN I >HI ( IN. JNl, 
CHARIIN,JNl,HX|Ili,J*l 
no tog isi, in- 
oo ICO j= 1 . jn 


|F^IM*RII, JI.C0.01 f.o TO 


00 

if 

ir 

IF 

I< T 


*11. J. I 1=0.0 
1 In k=2,KN 

(INARI I , Jl .CO. 11 00 10 
t MARI I ,J) .CQ.2I 00 TC 
IMRII.JI.rq, 31 00 TO 
ikiri I ,ji .rc.ot oo ig 
INIKI I ,JI .EC. 51 GO 10 
IhlRI ! , Jl .EQ.bl GO 10 

iHARiiiJi.ro. 7i no io 
fHIRI i ,ji .rv.n oo 10 
IMRII JI.C0.9I GO 10 


101 

n 

13 

IN 

15 

lb 

17 

ia 

19 


11 


12 


13 


IN 


15 


lb 


17 


1* 


IN 


20 


102 


UC 

101 

IDO 


210 

201 

200 


IF 1MRit.JI.CO.10l 00 10 20 

Ax:02>IR0ll*l,J.HMR0l2«t,J.N-ll-R0ll-l,J,HI>R0li-l,J,N-lMF<N*D 
00 TO 102 

CONTINUE , „ 

RX=02*IROII.l,J,K»*ROII*l,J,K-ll-ROII-l,J,KI-POII-l.J,K-lM/IN*» 
00 TO 102 

RXSIDZ* IRO |I*l,J,XI»ROIl»l,J»H-ll-ROII-l,J,K I -ROI I-l.J.M-ll 1/ IN •• 
GO TO 102 
C ON T IlluC 

AXS02*IN*IROII«l.J,KI«ROII*1<J,K-lM~l*IROIl,J,K)*ROt1,J,K-lll 
C- (ROII *2, J,K 1 ♦Roll*?, J.K-ll I 1/16*0X1 

RXSoi*V^*lRO|l.J.K|.ROII t J,M-lll*IROII-2,J l K|«ROIl-2«J,M-1l I 

C-N* IRfl 1 1 - | . J,K l.fiOl !-! , J s K-« | ! ! i * S=>DX 1 


00 TO 102 
CONTINUC 
RX=D7*I«*|R0II* 


J.KCROI 1*1, J.H-1II-3*IR0II,J,KI«R0II,J,K>1M 


RXS02*|N*|R0II*f ,J,K|.ROI 1*1, J.H-1M 
C-IRO 1 1*2. J,K | * ROI t*?,J,K-l) I I / I N*OX I 
OO TO ID* 


CONTINUE 

RXs02*lR£ 


__ - 1Ol1«l l J t K>* , ROII«l.J.II-ll-ROII-l,J > KI-ROII-lfJtH-lll/!N*0 
00 TO 102 
CONTINUC 

RX:2)2*IN*|ROII*l,J,Kt*ROII*I ( J(K-|il-3*iROII,J,Kl*ROII,J,N-III 
C- 1 ROI I *2, J ,N)«R0(1 *2, J.K-ll 1 1 /I 6*0X1 
GO TO 102 
CONTINUC 

KX=D2«|R0lt*l,J,H|.R0l I » 1 , J ,M -1 1 -RO 1 1 -1 , J.H I -KOI £-1 » J.K-1 1 I / 1 6*0 
00 TO 102 ’ 

CONI INUC 

RX=02*l3*tfl0!l ( J,KI*R0II,J,K-lVl«IR0II-2,J,K1*R0II-2.J,M-l 11 
C— N* < RO 1 1 -l , J ,K 7 *N0 1 1 — I , J , K — | 1 M /I 6*0X1 
00 TO 102 
CONTINUC 

9X=0Z*I3*IR0H , J.K l*RO(l,J ( K-lll«IROII-2,J,KI*ROII-2.J,K-I)l 
C-N*IRO(I-l ,J,KI*R0II-1 ,J,H-1M I/I 6*0X1 
00 TO 102 
CONTINUE 

RINTXII,J,KI=RINTXII,J,K-1I«RX*HI I I , Jl « IRO 1 1 , J ,N I *ROI I, J.H-ll l*P 
CM X I I • J 1 / 2 *0 
CONTINUE 
CONTINUE 
CONTINUC 
00 200 1=1, IN 
00 200 J=1 ,JN 

IF I HA RlI.JI.CO. 01 GO 10 701 
DO 210 H =2 , KN 

RlN?Xll,J,Kl=RINTXII,J,fcl-IK-l I*02*I1XII.JI*IR0I1, J,KI«ROI I, J,K-| 
C / 2.0 
CONTINUE 
CONTINUE 
CONTINUC 
DO 300 ISI, IN 
30 300 Js 1 , JN 

IF (HARM jt.ER.OI GO TO 301 


80 


?• 

in 

11 

a? 

4* 

it 

it 

ii 


310 

301 


ftSOHXStff infill ,J,K »«MIMTjMI f JtK-1) 1*1 M/2) • AR*CUL*Ml ( I i J) 
K iNT JI ! J ) SKINT (l Jjmrsuhx 
continue 
CONTINUE 
CONTINUE 
RETURN 
.END 


AI‘*NAS*|I I. POINTY FM CREATE© ON 5 HA7 10 AT HtSAiO) 


« 

4 

7 

4 

9 

j: 

ii 

i? 

ii 

Is 

It 

it 

i» 

2" 

21 

22 

2! 

?? 

* 6 

27 

29 

24 

»•» 

«* 

It 

32 

I| 

3* 

J* 

34 

37 

ie 

34 

4P 

41 

42 

43 
*« 

45 

46 

47 
41 
«4 
r«* 

51 

52 

53 

54 

55 
54 
57 
u 
«3 
6 ? 
61 
62 

63 

64 

65 

ti6 

67 

69 

69 

70 

73 
72 
7 ! 

74 

75 

76 

77 
75 


C THIS SUAROII T INC CONFUTES 7P IN THE POISSONS EQUATION 

SUfcPOUf 1NE R0INT7II,J,K,IN, JN , KN , DX ,07 , 0? ,R0 , AP.EUL ,HI ,MA* , 

"IsrNJlON^HINm IN.JN.KNI ,ROI IN , JN ,MN I , 7 1 N T f IN , JN I ,HI t IN, JN I , 
CM7IJN.JNI ,H«A(1N,JN1 
DO | 00 I = j , IN 
J; 1 , » 


DO 

IP 


DO 
IF 
I r 
IF 
IF 
IF 

If 

if 

IF 


100 J=1,JN 
IN»R| I, JI.EC.OI 
-•"* ‘ 11 = 0.0 


RINT7 < ! , J, 

i in K=i,4N 


IN4RI ( ,JI .1 
IN4R I I ,J).I 

IH»RI I ,JI.L 

I H *R I 1 , J I .C Q « 4 1 
INARI I , Jl .CC.SI 
INARI I ,J» .£0.41 
IN *R I I , Jl .£0.7/ 
IH4RI I .JI.CQ.RI 
NAHI I I JI.E0.9I 


GO TO IUI 


0 12 
GO TO 13 
GO TO 14 
GO TO IS 
GO TO 16 
GO TO 17 
00 TO (• 
GO TO 19 
GO TO 20 


II 


12 


IS 


14 


IS 


16 


17 


16 


19 


271 


102 


1 1 U 

mi 

100 


210 

201 

200 


IF INARI I. Jl .E0.10I __ 

ft7=02*(ROi 1 ,J*I,K I *ROI J ,J*1 ,K-1I -ROI I ,J-1 ,KI-R01I »J-1 ,K-1 i ;.'(4*0 
GO TO 1Q2 
CONTINUE 

N7=02*«3*IR0II. J,KI*ROII f J t K-lll*(RO(I,J-2,K»«AOI I,J-2,X-1 II 
C"4*IR0lI,J“l,Hl*R0II,J-l,H-llM/l 4*07 I 
GO TO lOi 
CONTINUE 

R7=02*IN*lftOII, J«I,KI*R0I1,J«1,K-11 1-3*1001 I , J.K I «ROII , J, K -1 1 1 
C-IR0II.J«2,KI*R0II ,J«2,H-1I I 1/14*07 I 
cn in Ins 

CONTINUE 

R 7:02* IRO II,J*1.RI , ROI1,J*I,H-1)»ROII.J“I,K1“ROI1,J"I,K”IM/I4*0 

GO TO 102 

CONTINUE 

RTs02*IR0l!,J*l,M*P0U,J«l,K>ll-R0tI,J~l,KI-R0<l,J-!,K-ll 1/14*0 
60 TO 102 

f?SJJSVf • IROII.J,KMROII p J,»l-IM*lflOII,J-2 l KJ*ROII,J-2,H-lM 
C“4 * iro 1 1 , j-i ,Ki«AoiT,j-i,K>im/m*o7i 
GO TO 102 
CONTINUE 

R7:D2*IR0ll,J«l,Hl«R0(I,J*l,K-l)~R0ll,J-l,KI-R0(I,J-l,X-in/I4*P 

GO TO 102 

CONTINUE 

R 7=02* 14* IRO 1 1, J* I ,K I *R0 1 1 . J« 1 ,K - 1 1 1 -3*IR0 1 1 , J ,K I *R01 1 , J,K -1 1 » 
C" > IR0II,J*?,K)*R0(I,J«2,K-1)I)/I4*DTI 
GO 10 102 
CONTINUE 

R7=O2*IROII,J*l.K!*»0II,J*l,H-l»-ROII ,J-I .KI-ROII. j- 1 .X-ll 1/14*0 

GO TO 102 

CONTINUE 

R7=02*|4*lfl0II, J*l ,H»*R0l I, J* l.K-l I l-3*IR0(I, J,X I *ROI I , J.K-l 1 1 
C-IROI 1 , J.2,K»*R0I ! , J«2,K-1 I 1 1/(4*071 
GO TO 102 
CONTINUE 

RY:D 2*13*1 R0IIiJ,KI*R0II,J t K-|)l*(R0ll,J~2,KI*R0( I, J~2,K“1 1 1 
r-4*iR0II.J-l,HI*R0< J,J-1,«-1M1/I4*07» 

GO 70 102 
CONTINUE 

R INI 7 I 1 . J,K ISPINTV I1,J,N-1I*RV*H1 I I ,J 1 « IRO 1 1 , J,K I *R0l I , J, K -I I 1*0 
ChYIIjJl/2.0 
CONTINUE 
CONTINUE 
CONTINUE 
00 200 I:l,IN 
00 203 J : 1 , JN 

IF INXRCI. JI.EC.OI GO TO 201 
DO 210 K-Z , KN 

ffINT7fI,J.K) = RINT7IIyJ,l<l a> IH~tl*D7*H7|I,t)i*IR0iI,J,KI*R0II,J,K-l 
C/2.0 
continue 

CONTINUE 
CONTINUE 
CO 300 1=1, IN 
00 300 J=1.JN 


82 


If 


so 


•i 


•2 


■ I 

310 

a*. 

301 

IS 

300 

St 


SI 



ir iH«nn,ji.C9.oi go io iui 

SsUHV = |SlfWvn ,J|H H R 1NT V ( 1 * J » K — l * » * I 02/2 i *AP*CUt *MI< 1 » M 
VINT(IiJ|:v!NT(ItJ>*RSUNV 

coHrl»Juc 


CONTINUE 

RETURN 

CNO 
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ASA'NA$Am.R»H r OR CRCATEO ON 5 HAY 10 AT lil56>lt 


*1 

J 

H 

s 

6 

T 

I 

* 

I? 

12 

U 3 

IS 

I* 

‘T 

2 1 

22 
23 

it 

26 

u 

!3 


2 IMIS SUBROUT INC COMPUTES VCRTICAl, VCLOCIUCS AT HALT 

c«« 5 » ,D * p 2 i! a , 5 .*** aa *•••*•**•*«•*»**«•»* *•*»»*•••••*«**« 

C SUBROUTINE RUMII,J f K,IU,JH,IN,JN,KN,IWN f JWN.U,V,l*H,HI,OX,OYiD7, 

gi^ENsioN HRHllMW;tjMN| V,,Nf ‘ ,N,MM * * MM< * W N t J“N • HN t * HI 1 1 N , JN » 
KNMlRNN-l. 

00 IQ IW- I | IWN 

00 ID JWEl JUN 

$5 raM*' 0 * 0 ’ 60 f0 • 

k:hn-«D«I 

1 = IW 
JSJU 

C-m II ,J|i|UtI , J,M«UI I, 4 ,K-t M 1 / 11 * 0 X 1 
01 HVY=|Ml iIO.j'tl'ml’UJH |HI*VU« If JM,K-m«HtlI,J«t I* , 

Cl VII .J*I»K)*yt 1 , JM.K-n »-MI( 1 * 1 , J I * I V 1 1 « 1 
c-rfii i.ji*m i,j,k) *vi ijj.K-n n/Ii*oyi 

UHIIU,JU,H-I)=WMllw, JWfKMli.O/HH^IOlMUX + OlHim^OZ 


9 

a 

to 


CONTINUE 

CONTINUE 

CONTINUE 

RETURN 

END 


84 






JL! J5CAIC0 ON 5 MY JO^AT J 1 « 5§ 1 57 

- J 
•» 

5 “ 'SUBROUTINE BURI I 

OINCMSION U 1 1 N » 

CHUIN.JNI.HjnN, 

00 to 1 : 1)111 

12 .1- I ! 1M 


* 

B 

n 

to 

it 

i! 

i! 

!• 

I* 

■P 


£ U *^M^^SUBB«^J^^*^CoSplJTC^RU^*VEB^^^5ctOC^^t5 M * 

C INI (.SB At GRID POINTS* 

•dS t *b! J5i ! 8*?5fiSg t jK i J§ i 3i ;SS 1 w* < > • 

WB I ? »5«H I 5*IM“ 1 1 *OZ*l Ul I » J»K t *HX 1 1 » Jl * V 1 1 » J»K t 1 1 » Jt I *H1 1 1 # Jt 


B 

I 

10 


C*U I T 1 J • H I 
CON 1 |NUE 
CONTINUE 
CONTINUE 
RETURN 
END 




I !•!! l8fSf¥ ,IMC c ° M<ims * l,t ’■'"<“ ukci'io •> «*ir 

_tOBf IN, JN.HK.lMNtJHNtUiWtWMiMI, ««♦«»» 

1 1 WM| iMNtJUWyNNI »NI HMf JNI • 

Vtvjwni 

iJMAfKNl 


j* 


k&jH.s} towwH&iSrty ’ 1 

DIMENSION URhI IWNi jun.kni 
KNHISKN-I * 

“ lo ■ 


t 

CH ni 

• MNM 

oo 


1 MN 

00 10 JUSt.JNN 

ir ihohiim. Jut.ca.oiio to • 


f 

10 


HIUtfSIHttl 
MflVSIHT t 

hi*v:inI I 
CO 9 MS| « 
i=rv 

JSJH 
UAV 
¥*¥ 

UOH 


JW I •LO • . _ 


St.MNM 


I I V, JU,K IS|K-n>DI»(uU*HX«V*¥AV*NV 4 VI«H.t< V*WHI IU»JW,KI 


i! 


jm* 

INU 


CON 
CON 
CONTINU 
RETURN 
END 


ASA-MSM II .ST0RC2 FOR CREATED ON $ NAT 80 AT >2102130 


v 

« 

' A. 

7 

8 
9 

t? 

{] 

13 

I? 

16 

17 

IF 

19 

s? 

22 

23 

29 

25 

IS 

28 

29 

JO 

u 


t” THI$*PR0GfiAN*ST0RE5*TMC*R£lEVAN T *DATA INTO*f ILE*I ************* 

*** SUBROUTINE S T ORE 2 1 U , W , WmJp 1 N 1 H , I, J , K , 1 U , JU , IN , JN • AN , I UN , JUN , 0 » E , 

|Tj ilTA«Ii3«Ai;??i4 UWiJg :2:<E£ i8ie8'#j£“?4 t*«* 

cS37?BS;38«V3iT:s!?sT£!!i«A:56riT**' ,, * 0,,H - J 

DIMENSION Ht I IN, JNI ,HT I IN,JNI .HI IIN.JNI .MARI IN, JNI ,HRHI TUN, JMNI , 
CU I IN, JN.KNI ,WR I IN, JN.KNI . URH II WN, JUN, KN I 
OIMENSION rliN, jn.kni ,ro|in,jn,kni , TUI IUN, JWN,KNI ,R0UI INN, JUN,KM 
WRITE 187 1 1 1 U I I «JiK),k:1,knI«J£|,JN) ,1:1, INI, 

Cl I I V 1 1 , J , K I , K 5 I ,KNI , J" 1 , JN 1 ,1-1 , IN I , 

CH iOII,J.KI,K=l,HNI,J21.JN»,I=l,IN>» 
CUIt!l,J,«),K:i,(IN),j:iJNI,!:l ( iNI, . 

Cl I IWHlIU, JW,H | ,H = 1 ,HNI , JW S 1 , JWN f . 1 W= 1 , IV 
CIIIWU,U,M.K?[ f KN|.J:I,JNIlrl.]NI. 

Cl I IWfll I ,J,Kl ,K 2 I ,KNl , jsl ,JNl , I = f * INI, 

CIIIWRM||w,JW,h 7,K=I (KNlfJWTti JUN I>1W-1,IUN), 

Cl IPINTMII U, JVI , JWC1 ,JWN) ,IWi| »IWNI 

C, I IMUI,JI,J=I,JNl,I2t, INI.I I MX 4 I . Jl,J=l, JNI,I = 1,INI.I (HYII.JI.J 
Cl ,JNI,I=l,lNI .IIMARII,JI,j:i,JNI.f :i,INI,(IHRHIlW,JWI,JWr|,jUNI, 
C I W- 1 , 1 WN I ; ( I I flI,J,K»,K:I ,HN | . j=j , JN J ,1 = 1, IN I, 
CIIIROI|,J,KI,M=l,KNI,J=l,JN)|f-l,INI( 

CIIITWI1W,JU,KI,K-1,KNI . JW = 1 .JWN I ,IU = I .IUNl . 

CU I ROW I IU. JU.H 1 ,Krl ,KNI . JUU . JWNI , IW=I , IUNl , 

CT»I.TAH,TAV,AKT.CB,CW.A,B,C,EUL,7, TW.RO, ROU.TC ,RREF,TREF,TO,TAHO 
CAIfAH]Av,AP,OK,OY,0?,DT,TAUX l TAUT,TTOT 
END FILE 8 
RETURN 
ENO 


EWNI , 
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7.UMB2 FOR CRCA TCO ON 5 MY 60 AT 12j04|2Q 


.1 

i 

3 

4 
c 

6 

7 

I 

9 

IT 

11 

I? 

14 

is 

is 

17 

IP 

]9 

i'1 
21 
22 
2 3 

24 

25 
2b 
27 

3? 

31 

32 

33 

34 

15 

it 

37 

J5 

39 

«n 

4 1 
•12 

43 

44 

45 

45 
47 

46 
49 
tn 

li 

5! 

54 

55 
5b 
£7 
59 
59 

2 ? 

62 

63 

6* 

65 

6b 

67 

b« 

?? 

71 

72 

73 

74 

75 
74 

77 

78 


?.** THIS* SU3ROUT1NC*COmPU7CS*b6uNUARY* TEMPERATURES ** ** 


n 


12 


13 


14 


15 


17 


SUBROUTINE TCM|2I I . J,K , IN, JN.KN.TO. 
CTAMB ,HX,IIY,T,TRCF IV, I Al.TAM,B3,Dl 


UX,DY,D2,HAR,CU I HI,AR1,CW. 


. i •nD,n«,n, , i , mir • i * » , mi, i • n , u j , u ■ < 

DIMENSION T(IN,JN,KN) , TOI IN, JN.HNI ,MARI1N,JN> ,MXI IN.JNI ,MY( IN.J' 
CH I 1 I N , JN I 
KNNl = KN-I • 

do lun K = 1 , NN 
00 100 l = 1 , IN 
00 100 J= 1 , JN 


DIHTUXSQ.O 
D1MT VY=0.0 


01TWZS0.0 
I f IMRI I , JI.CQ.O) GO TO 300 
If IMRlI JI.EU.lil GO TO 300 


if 

ir 

if 

IF 

IF 

if 


IMM|,J».EQ.1I GO 

" I , J1 .£0.2! p 

GO 
GO 
GO 


(MAPI 1 , J! .£0.21 
I MR I 1 , J I •£ 0. 3 1 
IMRI I , Jl .C0.4I 
IMRI I , Jl .FQ.5I 
IMRII.JI.CQ.^I 


TO 11 
TO 12 
TO 13 
TO 14 
TO 15 
TO 16 
GO TO 17 


(MARI I ,J) ,EQ. 

IF lMRl!,JI.EQ.8| GO TO it 

If imriiJji.cq.ri co to 19 

jg | r .i;- J l.e,.lDI GO TO 20 

UlTX IWm, j.k, -Tll-1, J, Ml 1/12*0X1 
02TX = ITmi t J f K)«TtI-l,J,K>-2*T(!, J,KII/IDX*OXI 
02T2=(T( I ,J,K*11 i'TII, J,K-1 1-2*111, J,M 11/(02*021 
01 T Y =0 .0 

02TYs2*(t(I,J-1.KI-T(I,J,K||/|D2*027 

IF IK.CQ.lt GO TO 1 1 Q 

IF IK.EO.KNI GO TO 120 

GO TO 200 

CCNT1MU? 

01TX = lllI«l,J,KI-T(I-|. ,J,H I 1/12*0X1 
02TXr|ToNl t J,KMT(I-l ,J,K |-2*T| I, J,G»I/10X*0XI 
D? l2 = I T I I , J» N >* i J « T 1 1 , «J , K - l I "2 *7 1 1 , J .K I l/IOZ*DZI 
D 1 7 Y=0 .0 

P2TY = 2*1TU,J*1,MI-T1I,J,K1 »/IOV*OY> 

If in.eo.ii co fo no 

IF 1N.E0.KN I GO TO' 120 

GO TO 200 

CONTINUE 

OlTXsOvO 

0ZTX=2*|TlI«l.J,K»-7|I,J,Kn/10X*0XI 

0 2 72= 1 71 1 , J,«* ll *T II, J,N -1> -2*1 II, J,K 1 1/(02*021 
. 01TY:ITI!,J«1,KI-TII,JM,K> l/!2*0Yt 

D2!Y = I11I.J«1,NI«T1 I, J-l.KI -2*711, J,KI I / 1 0 Y*OV » 

if iK.ca.ii go to no 

IF IK.EO.KNI GG 70 120 
GO TG 200 
CONTINUE 
0 1 TX=0»0 

D2TX:2*lTlI-l l J f KI-7(I l J,KII/(DX*OXI 

D2TZ = IT(I,J,HMl«T|I,J t K-ll-2*TH,J,Kll/IOZ*OZI 

OlTVrlll J,J«I,K»-TII,J-1,KI l/(2*DYl 

C2TV=ITII,J*I,HI»TI J, J-1,NI-2*T«I,J,KII/I0V*DVI 

IF IK. CO. II GO TO ISO 

IF IK.CQ.HNI GO TO 120 

GO TO ?00 

CONTINUE 

0 1 YX =0.0 

C2TX:2*| Tl I«1 ,.!,K l-TII , J,K I I/I0X*0XI 

0272 = 1711 , J.K* U <7 |I, J.K-1 1-2*7 I I , J.K I l/IDZ*02 I 

0 1 7 Y =C .0 

IF IK.EO.KNI GO 70 120 

GO TO 200 

CONTINUE 

0178=0.0 

8 2TX;2*|T|I*1,J,KI-TII,J,K | I / I0Y*DXI 
212=17(1, J,N •11*711, J,K-1)-2*I(I,J,NI)/(DZ*D2 I 
Ql TT*0 •() 

1, J,KI )/(DV*OV> 
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IS 

u 

ii 

44 

15 

16 
<7 
41 
If 
SO 
91 
97 
9 J 

94 

95 

96 

97 
91 
9f 

i cn 
101 
107 

184 

105 

106 
1 07 
101 

109 

110 
111 
112 
111 
l 14 
I IS 
I It 
I 17 


If 


20 


110 


120 


200 


16 

It 

3QU 

100 


fr !8»Io.*i<iii G 8o T to l }?o 

81 T? s I T I 1 i I ♦ { ll| jIk-i! ^ItII^JiO l/ID2fD2l 

QITYSO.O 

D2TY=2*»T II»J«1 1 KI»TII, J.KI l/IDV*DYl 
I r u. co . ii oo fo no 
IF IK.CQ.KNI CO TO 120 

00 TO 200 
CONTINUC 
01 T* : 0.0 

02T*:2»(T 1 1 - 1 , J ,K I -T 1 1 , J,Kl| / IOX*OM 
O2T2=(TII,J,K<ll<T|I ( J,K-ll-2*TtI,J,K>»/l02*02l 

01 I V SO • 0 

D2IY = 2«ITII».)-1.K|-TII,J,KII/I02*0*.I 

if ik. co , ii oo to no 

IF IK.CO.KNI GO TO 120 

30 TO 200 

CONTINUC 

CT = *KT*I I T SI • J , 1 I’TRCr «TJtCn-TANa /TftCF 

C T "C " *H III J I 

gjT|s2^7({, J ,2l- CT *l)2-TII,J t lll/l02* D 2l 

CONTINUC 

02 TZ:2*IT|I,J,K-1I-TII.J|HI I/I 02*021 
GO TO 200 
CONTINUC 


TOII .J|K 1 = 11, O/IM II.JI I * I n T 6 I * 1 0 III TUX * 01HT V V *11 1 1 I I Jl f 0 1 TW2 I * T 6M 
•<QlTX»HXlI.Jl<D2TX»HI II |JI <01 TY*HYI l , Jl <02TY*HI I i , J I T » I T AV/MI I l , 
ll*B3*02T2l«0T«TII|JiKI 


C • 

Cl l*BJ* 
CONTINUC 
CONTINUE 
CONTINUC 
CONTINUC 
8ETURN 
CNO 
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roP created on s may #d at i ? t 07 1 1«> 


c» 

c 

c 

c* 

c 


THIS 

OOIO 


SUBROUTINE 

POINTS. 

«••••***** 


COMPUTES TEMPCOATURCS at interior 


1 « 


II 


12 


11 


IN 

IS 


15 

IT 

100 

IB 


19 


2 ? 


21 


22 


21 


I" 


25 


25 


27 


/X 


29 


SO 


SI 


S 2 


S 3 

20 U 

J 4 


3 S 


Jh 


S 7 

20 

3 * 


39 


90 


Ml 

“1 

41 

24 

4 N 


MS 


46 

21 

47 


48 


49 


SC 

a 

51 

9 

52 

n 

51 





SUBROUTINE TCHmi'J'H'PtiJN.lVN'U.V.TiTOoDX' 

CCB « 

CUT.O?.W.OT.TAI,TAH,T/lViBJ.HI.HX.HY.HAR.AKr t TRCr t TAHei 
DIMENSION U(IN,JN,KNI r V< IN.JN.HNI.HIIIN.JHI ,HXIIN*JNI »HY(lN*,tNl. 
CNXfll IN.JNl ,T| IN,JN,KN1 ,TOI In, JN,KN),W( IN,JN,KN» 

*NMl=XN«l • 

00 10 1=1, IN 
00 10 J : 1 1 JN 

IT lHXRII 1 JI.CO. 6 l GO TO 100 

ir'?A«!lTjl!?l!ii?°G0°T0 0 T 

CONTINUE 

8lHTUXr(t5l j«J jJjK(*TII*|,J t Kl*MI(IfI,JI-U(I-l l J,KI*TIT-lf J «K I 

C 01HTVr:!viI.J»I.Kl«TII»J*l,Kl»Hl(I.J*l)-VII t J>l,HI«TII l U-l ( KI* 
CHI(I,U-1»|/U*0TI 
OlTx:<TlI«l,JtMl-TII-l,J«Kn/l2’'QXI 
OITy3ITII,J«I.Kl-T(|,J-l,Hn/I2*llVI 
D1TU2:( T l*utf . J t N* ll -Tl t,J'H-l)*WI I.JiK- 

D2TX5(T(I«I»J,K)»T(I"I|J(KI-2*T(I f J|KIl/(0X»0XI 
02TY=(TI I,J«I*KJ «T(I, J-l ,« 1-2*I( I , J,K ) 1 / ( OY*OY 1 
Q2l2=(f (J, J,H*ll«TU,J.K-l 1-2*1 f I J.K ) 1/102*02) 

IT IHAR 1 1 .JI.EQ • 1 1 1 GO TO 200 
01HTUX=0.0 
DIHTVVSO.O 
0 1 TW2;0 .0 
CONTINUE 

IT (K.EQ.ll GO TO 24 
IT IK.EQ.HNI GO TO 20 
GO TO 21 
CONTINUE 
D1?W2=0.0 

D2T2=S!'*I Tll,J.K-l)-TII,J,K)*C8*H!l I , J » *02 1 / ( 02 *02 I 
G 0 TO 2 1 
CONTINUE 

CTiAKT*fIT(Ijj*|i*TRET»TRErl'”TA M BlT7REr 
C T r C T »H I I I , J [ 

riiTu»-nn’ • 


,J,K-I1 1/(2*021 


01TU2=0.0 

02T2r2*ell I, J, 21 -CT*D2 -T ( I , J, 11 1/(02*021 
CONTINUC 


TO 1 1 , J ,K I : • 1 .0/H I ( 1 , J I 1 * I •! A I * ( D1HTUX *0 1HT VY«HI 1 1 » Jl *0 1 TW2 1 
C«TAM*(0ITX*HX(I.Jl*02TX*HIII, J I * D 1 TV »HV 1 1 • J I *D2T Y*HI < X t J) I 
C* ITAV/HII I ,JI J*fl3*D2T2 » *0 T ♦ T i I,J,*I 


CONT INUE 

CONTINUE 

.CONTINUE 

RETURN 

END 


90 


5**ps*m 


i? 

i! 

i! 

it 

19 

20 
21 
22 
23 

n 

26 

n 

2B 

29 

30 

n 

a 

35 
3b 
3 T 

jS 

90 

u I 
12 
1 1 
19 

• 5 


t nOB r OR CREATEO ON 5 N»V 10 AT 12l09tl* 

C THIS SUBROUTINE ALLOWS FOR VERTICAL MIXING AT A 
C PARIlCULAR GRID POINT. 

C ■•••*» ***•••••••*•♦ ••***•♦•♦♦ 

" SUBROUTINE TCOBIltJ.Kf IN.JN.KN'TtHARI 

0 ! H|N5^0N T | IN, JN.KN) ,MRR( IN, JN1 

00 10 I: 1 , IN 
io j=ii 


no 


00 to JM.JN 
I r ‘IHARII.JI.CO.OI 00 TO 9 

continue 

M ARK SO 

00 I KSi.KNNl 

TXSO. 01573 

TTsTl I.J.KI'TX 

IF lM0.lt GO TO 7 

IF IK.CO.KNHIi 00 TO A 

IF nf ,GE. Til, J, Mill GO TO 111 

HARK "1 

HW"”" 

S? ^ ( T T^GE . T 1 1 i J.R ♦ 1 1 1 GO TO 111 
M ARK S 1 

avT:|TT»2*T(I,J,K*1M/3.0 

TIIiJ.KIsAVT-TX 

TlI.J.KillsAVf 

GO TO 5 

CONTINUE 

IF |TT,GE.TII,J,K«1I> GO TO 111 
MARKS! 

AvT=i2*rr*Tii|j,K*i)i/3.o 
HI,J.K> = AVT-fx 
nif^KHisAvi 


5 

111 

CONTINUE 

CONTINUE 

■ 

CONTINUE 

3 

CONTINUE 

10 

CONTINUE 


RETURN 

END 
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AS* ‘NASA! n.T»» INK row CREATCO ON t« NOV >f AT 


H 

. !> 
A 
7 
i 
a 

1C 
It 
I? 
I J 
in 
15 
1 A 
17 
IK 
:« 
2V 
21 
27 

23 
2« 
25 

24 

27 

28 
29 

ii 


C THIS I'WOORAH PRIHTi IHC VAIUC5 or r,TV,RO,ROW 

c for uahc HCkocc 

SUIROUTINC tpRlHKII.J|K/IN,JN.KN»I,RQ. iref »mar,tactuu 
D 1 H 1 NS 10 N 1 ( IN I JN| NN 1 ! no ( IN,i/N,HNI,MA/ 2 |lN,JNl f f AC T UL 1 1 N , JN , KN 1 
If IkNiLCiti 00 TO lOt 
00 10 k ; I ,MN 
DO 10 I : I • IN 

, PRINT U*K.-1.ITI!.J.KI .JSIfJNI 

ro PRINT 12,(R0I I,J»KI,J=I,JN1 

11 T OHHAT I / • M S • , 1 3 1 jxj'ls' .13/* TCNPCRAniRC , /l5)l.«C15.7ll 

12 rORHATI • DENSITY V|5X,BE15.7I I 

101 continue . 

go loo h:),mn 


inn 


mo 

153 
1 05 
10b 


0 JQO JS1|JN 

K I : T 1 1 1 JjK > 


00 lOO J=ltl« 
ITUL ( ! , j J 
LI * ‘ 

continue 


tactuli 

TACTUL I 
IT INARt 


i i jii\ li i i | iJ<n i 

IJMSli .♦TACTUll I,J,K 11* TREE 

r.jf.EO.OI TACTULI i!j,k mioooo 


000.00 


00 ISO k: 1 .kn 

WRITClb.lOil K 
00 mo IriyIN 

CONTlAui 10 ^' ,TACTUt,1I * J,,U ' J^1,JN, 

CONTINUE 

FORM ATI* l 1 .* TEMPERATURES FOR KS»,I 51 
F ORTJAT ( // , 22 rb.it 
RETURN 
C NO 
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IJlAS,U 1 1 


7 

n 

1 

I! 

u 

i* 

it 

17 

1! 

1«5 

20 

21 

5 ? 

:» 

2S 

2* 

27 

2 ! 

29 

1? 

11 

n 

a 

15 
!*> 
!7 
»e 
29 
<*? 
m 
• 2 
•»! 
44 
4! 

«5 

47 

49 

*9 

C 4 

si 

C 4 

53 

£4 

54 

it 

c • 

si 


.UV f r OB CRCAUO ON 5 

c* f rHirsuSjioSMNrcgpuiEry**^* 


MAY HO AT 1 2 I I 3 1 03 

*»***«******4»'>4**44**«44*4S>4AlVAe44 

AND V /OR VARIABLE DENSITY 

***»**o*****«*»*4**t**»r«*** v 


8 

9 

IC 


Sl'OR.OUTINt 


II.JfKi IW, JW, tw, JM,KN, iWN,JWN* 0, V,0|C .M.G.DX .Dy.OZ. 


CRiMfX, MINTY 1 , 1 Ini] 

col, AI 1 AP,AH.AV,A3,HI.HX.HV,7',MAR» 

DIME Nil ON ul IN » JN , KN| , V 1 1 N , JN ,KN ) , 0 II N * JN , KNI , E J IN , JN , IjN J» 

CM I IN, JN.HNl ,01 IN, JN,HN»,HII IN , JN ) , HX ( I N , JN » , M Y T I N , JN» t R 1 1 UN . .JWN > » 
CHARI IV, JNI 

dimension wiin.jn.kni 

DIMENSION A 3 t m M i 

DIMENSION RINTxUN,JN*HN> iRTNTYI IN.JN.KNI 


ION 

NMM I *NN“ I . 
A SOT «AM* I 1 


5=1 


/IOX*OXMl/IOV40Y7l 
• IN 
,JN 


00 10 
00 10 
iw;i 

jr'fMARII.JI.tT.m 00 TO 9 
CO e H= 2 ,KNM 1 

DIPJCIPI IN, JMI-PIXW-I .JUMP I IW, JW -1 »"PlIV-i,JU‘l 1 1 / 1240 X 1 
C«CUL* r )INTxll,J,Kt 

CIPYSiPI IN,JNI-PIIM»JU- 1 I*P« JVM, Jtt)-P(iw -1 ,Jtf-l»l/i 2 *DY> 
C*£UI*R 1 NIY 1 1 7 j ,K I 
|S 0 T«AV 9 A}(K j/l 02 *D 2 » 

C:<HllI f JMA*HI(J r J)«a/HI(l.JI» 

0 I MUUXS I U 1 1 ♦ I * J • K 1 9 U 1 1 ♦ I » J ».H I' *MI I X ♦ 1 » J I *Wl I "I » J»K I 
C*UCI -1 ,J,K | • Ml t f- 1 |JII / li*Ox I 
DIHUVYSCUI I * J* I ,H livi I, Jfl.K > AMI 1 1 ,JUI-UU, J- 1 ,K I 
C»V 1 1 »J-t *K I *HI 1 1 , J -1 1 1 / I 2 * 0 Y I 

C*V ll-t ,J,K I 4 HIU -1 ,JI 1 /I 2 * 0 X 1 
0 IHVVY 5 IVII,J«i f KI*VI I ,J*i,K l*HTI I I ,J-l,K »4 
CVII,J-l,Kl*HIlij-ll|/ 29 DYI 


U'UM 

DlUX : IUlI*ItJ,KI"-UII-I,J,KII/j2*0XI 
01UY = IUI I, J«l, Kl -UU,J- !,«(»/( 2*0 Y> 

OlVXSIVI -VI I-l.J.KlI/l 240X1 
0iVY=iViI.j«l.HI-VII,J-liHI>/l2«0YI 
01UN2=IUI I , J,K «1 > *WI I , J,K • i I -U( I , J,7 

OIVU2=IVII.J,X«lJ*UII.J,H«I>“VII,J,N-ll*HII,J ? K-in/l2*027 
OOUXSIUI l ♦ I , J , X I 4UI I - I , J,K>-OfI,J,K»W 10X90X1 
DdU V = I U ( I » J * 1 » K > * U I I , J-l , K I —0 I I . J , K I I / I DY *0 V I 
C0U2SIUI I ,J,KM! «UII,J,M-II>OI I ,J,MM/I02*D2I 


■Utl,J,K-l)4WU,J,K-l 1 1/12402) 


CDV 7 :|V(I«l t J,KMV(I-I,J f K)-E(!,J,K)t/(QX* 0 XJ 
OOVYStV(I,J'l,KI«VII|J-l,Mi-C(l,J f H)|/tOY*OYI 
00 V 2 =I VII. J,H* 1 ) ♦ V| I.J.X-l I -EM ,J,M )/ I 0 Z* 0 Z» 

H 1 1 , J , k I S 1 0 T/C 1 • I -A I • I 0 1HUUX « 01HU V Y«Hl 1 1 , J 1 401 UUZ 1 -HI 1 1 , J I 4 AP 
C«0iPX«AH*MlUtJ)9l00UX4D0UY>tAH9HX 1 1 ,U J 40 1 UX ♦ AM4H V | Z , J > 40 1 UY 
t‘AV4AJIK I A0DU2/HI I 1 , J I > • HI I 1 1 J I *U I I , J • K I / C 
0 I I ,J,rt » = 1 0 7 /C » 4 ( -Al4|01HUVX*0lHVVY*HII I, J I 4DI VUZ I -HI I I, JIMP 
C401PV*AH4HI I I,JM>IOOVX«DDVY IMH4HX II»J»40IVX* AH*HYI I* JI401VY 
C*AV*AJIKI •D0V2/HI (I,J))»HII1,J)*V(I,J,K)/C 
CONTINUE 
CONTINUE 
CONTINUE 
RETURN 
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NAs*>*MSMH.uvtop elt cmeateo on 19 nov 79 at »ho 9 j 19 

l C THIS PBOORAH CALCULATES U ANO V VELOCITIES AT THE SURrACE US 

i ? t t *.o*.***»22222*5j*i°#*i**2*2****t«i ♦♦*♦*♦****«*•*****«♦•♦•♦♦ 

*'****SUBROunNC*UVTOPIHtriTAu5*TAUT,I?*tHtDE,lN,JN,>(N*HI.MA«I 
. DIMENSION MniN^JNI.MARfJMiJHl.HIlN, JN,KN| ,G»1N, JN.HNI 
do aao I - 1 (if. 
do eno j:i.jn 

If IflAlfdiUliLT.tll CO TO TOO 

M=f »UX*HIII,JI 
■ TVS TAUT *HII I * Jl 

CO 1 f Vf K C = 1 I J, k <,K«2)-2*PZ*TYJ/J. 

CONTINUE 
TETUftM 
ENO 


S 

« 

5 

6 
T 

a 

9 

1 ? 

I! 


II 

19 

IS 

It 


too 

aou 
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I 


■N»*Am.WHAMJ FOR CREATED ON \H HAT 7M At IStSOllO 


jj *"jHIS PROGfi AM*CAL CUt A?ES*THE*V A*UC*Of l »J*FRON*VAt!uCS*OF*WH*AT 

oTO5i$M T 

IF mi.yirfr.ll. GO TO 35AO 

nn Kin «■ - i • u u 


00 36 10 K=1VKN 
lusl 

C/T • 

CONTINUE 

CONTINUE 


IS10 CONTINUE 
IS40 CONTINUE 

ish $mr 





* 


M* S* 

. ? 


it 

5 

6 
7 
a 
<> 

IP 

n 

15 

t" 


* 5 A ( 1 I * LHTOP FOP CPC AT CO ON IP HAT 7 R AT ISlSOiOO _ , 

C H'IS PB 00 P/IT 1 SC T S IHC VALUE OF WH EQUAL TO 7 ERO AT THE SURFACE 

o jSt S " " !«» ! IS^ltS l“ N ’ J “" 1 “* 1 

DIMENSION MPHUWN.JPNI 

" H88 iHSljMB 

OO TO JOfJO 


3000 
3 300 


DO 

00 ... 

ir (M*K/i;w,UW>.EQ.DI 
uh(!v..JuVi 1=0 
CONTINUE 

continue 

RETURN 

END 


LISTINGS OF PLOT PROGRAMS 
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9 C«H*t 65 &VI 


;4 

SB 

SI 

i t 

S 3 

So 

45 

56 

i? 

Sfl 

S? 

6 B 

61 

t 

o3 

64 

o5 

6 b 

6 ’ 

6i 

09 

73 

71 

73 

73 

7 «i 

75 

76 

77 
73 


it.PLotrckfr enema jh as auc «i *r izts&in* 

v ...«.«••.* >••■••<)* • 75 * •* 

C 


»»*« 


ThU faoGAAn PLOTS me 5UAF&CC ISOTHERMS FOR IMC REGION 


M C 


e v v 2 tol^» * RC Rt *° wl ™ ** 0|,CN ropMr *” 


PB 2 DISCHARGE JCnPCRATU»C. 

Ip, 11"' 


ro t*i. s 

,t ThE Ot HER 
*- USfRS MANUAL. 


"TO 0!HrNSI0N«LUe HHC 10 HOURS. 

IflULA TED I IHC I HOURS I «^»f _ HITS IS , NOT_RE A U ** 


VARIABLES HAVE UCLW EWSCRIft-EB IN THE 


P 

i 

/> ARAMETCR IN:i7,JH;Z0,IWtMU,JUM:i9 f kN=5,MH12^ 


5 1 



Wl 8 *« 1 QD 01 

acAD^sLpi^^pto.pRiPs.cTTorfNnHE 

CA?L*Jt 6 TSUBUF, 10 « 0 .a» 

NPLT=0 

MTIHLS* 

88ci°fit2V^3*AjVv |wH,PINTH.I , J.K.IH. JWt IN. JN.HN, 
CtVN> JUN.O.E, fix I h Y Jh 1, H AR , mAh. Ai ,*h7*«, a p. ox, 
~CO*.,OZ,Df,TAUX, fAUr»W»Wfi,WRM,TAI , T Ah , T A V , ART , 
C(JB.gwjAjB,C.EUL, T »rw,RO,R°u, T ErPREF,TREf ,70, 

cowtInue 7 

BALL. PLOTS II Bl/F, 1000, 1 1) 

CALL PLOTin c n s t!^ = 3 M* 

F CRN A I t I 

k=i 

P»zn 

P6=crior.noT 

Nl:O.D 
W 2 :o.O 
STARTSlO.S 


JUN, 


l ? 3 


03 335 J J: 1 , JH 

Hx'ait jjlzbcF'ii.tTen.jj.in 

ifihabcii.jji.eq.ogo ro 333 

IF)HX?13,15).GT. 10» SI START = 11.0 


CONTINUE 


CALL PUpn7.fl,OtO,-3l 

-CALL ECMKOy<Rx tAN.^M, t , IN , 1 , JN, « . 75, 4 .0 , 0 . 04 , S TART.q . 5 , 
C IH, JR , START, 19. 0,(5, 0,0.0, a. OfO.OiNl ,«2,0.Q7,I.0,NPLT) 


RCA 10 , 


C ThE NEXT 66 LINES ARE TOR DRA.ING THE BOUNDARIES Or ThE 
C C3HAIN. THESE LINES MUST BE CHANGED F OR ANY OTHER DOMAIN. 


CALL 

TALL 

CALL 

CALL 

Call 

call 

CALL 

c ctii 

m 

CALL 

CALL 

CALL 

CALL 

CALL 


PLOT ( »D .21^, “0 , 2$,\ "3 ) 


FAC TORI o 

PLorr 

plot i 

PLOT <13.0,4. 
PLD T I 1 3 ,8 . 
i'L 01117.0,8. 
PLOT (17.0,10 
PLOT ( 14,0,10 

e. L 2I ! i**:o 1 

PLOT 1 12.0 ,n 
PL0T(12.0,13 
PLOT * i%.0 ,13 
PLOT Ml. 0,14 
PLO T l It .0 , 1 4 


0> 3 I 
» 

I *) 
0,21 
0 , 2 ) 
. 0 , 2 ) 

: 8 ; I | 

. 0 , 2 ) 

. 0,21 

. 0 , 2 ) 

. 0 , 2 ) 

. 0 , 2 ) 


98 




/V 
4 *5 

kl 

*5 

hi 

6t 

«T 

<§ 

si 

s: 

hi 

94 

96 

hi 

sa 

99 

u: 

i ci 
1C.' 

111 

ICS 

ice 

1C7 

ica 

1 C 9 

u? 

112 
113 
1 in 

ill 

HI 

1 19 
122 
121 
1 Z2 
12 3 
124 
W5 

i 2’ 

124 

129 

132 
131 
13? 

133 
1 3<* 
1 35 

136 

137 
1 33 
1 39 
1h3 

1 4 1 
142 

142 

>45 
1 4C 
It? 
14« 
1 **9 

isr 

151 

111 

II? 

Ui 

lt» 


CALL rL0T(l(.Q,l«,t),2) 
cm monin.0ttti.oi2> 
C4LL PLOTIllt.Otl7.Ot2> 

C ILL PLOTlHtO.ti7.ui2> 
C1LL PLOntl!o,U,Ot2l 
CALL PL0Tl9.0.jS.0t2l 
CALL PLOn*.0,i6.Ot?> 
CALL PLOTtS-D, 16. 0,2) 

call plot la. o; it. o!5i 

EItt ct81i!B:8:H:S:s( 
citt StSTIf:3:Ig:B:i! 

CALL PLOT H, 0,20. 0,21 
call PLOT)*. C, 17.0,2) 
CALL PLOT 3.0.17.5,21 

SStfc PMl! :8ji8:S:!l 

CALL PL0TI1. 0,17. 0.2) 
CALL PL0T11. 0, 17.0,3) 
CALL PLOTU. 0,15.0,31 
CALL PLOtli. 0,15.0,2) 
CALL Pf.OTIl, 0,13. 0,2) 
CALL PI OT (1.0, 13.0,3) 
CALL PL3TI 1.0,11.0,31 
CALL PLOni.O, lt.0,2) 
CALL PLOTI 1.0, 9, 0,2) 
CALL PLOT I t .0, 9.0 ,3) 
CALL P.L0TI1. 0,7.0,31 
CALL PLOT! 1 . 0 , 7. 0, 2 ) 
CALL PLOTI3.0,7.Q,2> 
CALL PLOTt!. 0,5. 0,21 
CALI PLO Tl £ .0 , 5 . 0 , 2 ) 
CALL PLOT 15 ,0 ; 6sO § 2 1 
CALL PL0n7.0,b.O,2) 
CALL PLOT)]. 0,3.0. 21 

|itt itsfllBiB: :S:1I 

CALL PL0TI10.C,1.C,3) 
CALL F AC TOM 0 • 25 ) 

CALL PLOTlC.O,0.O,3> 
CALL FACTOPll.O) 

CALL PL0n0.0.“0.5,-J> 
CALL PLOTl-i.O.-l. 0,-31 
CALL PLOTlCi.fi, I. I, 51 
CALL PLOT tl.O, 1.1 ,21 
CALL PL0T1S. 0,9. 0,21 


call PLono.a.9.a,2i 

Eibt BtB?iB:8!i:i;ii 

THE NEXT*25 LINES ARE*Fo5*¥R I TING THe’caPTIONS^Of'the' 
PLOTS. THE SPECIFIC USER MUST SCRU TINIZE THESE LINES 


AND HAKE NECESSARY CHANGES. 


CALL Syn|OLI0.0,a.f.0.1K,23HFIG ISOTHERMS AT A- .0.0.25) 

call H«M6i:R(999..m.,fl.iK f pa,a.0,0i 

CALL STHBOLiO.O,a.3,u.lK,3oH LAKE KEOUEE- ( R IGIO-L 10 H0DELI,0. 


zz 


125 


CO, 36 ) 

I?(M.ir .2S.O)GO TO 22 

fALL.5»ffeOLU.0, 0.0,0. 14, 2ollilMULArlONS F OR F Ed, 26 1979, 

CQ.Cj.’c) 

GO 16 123 

CALL,, SYHBOL. (1.0,0. 0,0. In, 26HSI Mill A T IONS FOR F EB n 27 1979, 

C ?o8»?We 

call StNAOL ( 1.5,8. 7,0.1,! 2HR0N I/O; LCO ,0.0,12) 

Call m«h|c«i999. ,999.,o.i ,Pio,o.o,oi 

CALL STWff0Lll.S,a.5,Q.l,33«OiSCHAR6E VELOClU : 

C3| 

CALL NUMBER/ 3tB,S.&,o.l,Pl,O.Q, *2) 

CALL 5 v KI2Mi*5t2 ,3 »ft*i* 2 ? HDI SCHARGE TEMPERATURE: 

CALL SYHB0LiK. 2,R.H, 0,07,1 HO, 0.0, 1 1 
CALL NUMBER! 3.«,8. 3, 0.1 ,*2,0. Q,*ii 
CALL SYHR9Lll.5 v a.l,3.1,32H6lND SHE EG (MAX) s H/SEC, 0.0, 32 


CM/ SEC .0.0, 3 
C, 0.0,29) 
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IS 

no 

U\ 
I il 

us 

m 

us 

u* 

I 70 

1 J 1 

.psf.i 


no 


*'H. V 


t , 

‘iftt rawft:!:WftlJWibs# lOucAiet fi-oKn 
|{fct »« ■ 

cScl PL2T««!c;2ii:o;-i? l,iH * ,s *' 0 ’ 2 ' 

CONI I Hue 
CUD 


cn/ see, o.o,s 

HRS, 0.0,311 


100 


*CCW»SbW" 1 1 

'2 

3 

4 


i.fluv eir on rw *y to *r iiijoux 


ft 

7 

• 

9 

13 

u 

IS 

1? 

u 

17 

1 ft 

19 

ro 
2 1 
22 

Si 

25 

2ft 

27 

28 

29 

S3 

51 

5 2 

5 3 
34 
if 

JO 

37 

38 

39 
H 2 
Hi 
«2 
43 

4 H 

HS 

Hft 

H’ 

48 

49 

50 

51 

15 

SH 
< c 

si 

ST 

S3 

59 

63 

tl 

c.2 

c 3 


ftft 

t7 

6* 

fc9 

73 

73 
7 ? 
7 ! 

74 

75 

76 

77 
7 > 


THIS Pao&RAil PLOTS TriC U - V VELOCITIES FOR THE REGION 

n ft* J ^ y |* # *> 

THC F OLLouING VARIABLES APE REAO WITH AN OPEN FOPHAI 1 * 
S DISCHARGE velocity. 

= gischarge temperaturi 

* = S/nO^SPeEE * I MAXIMUM 1 1 
= CUPRENT. 


PI : DISCHARGE VELOCITY. 

P 2 : GISCHARGE TEMPERATURC. 

P 10 ; 

PH ' 

PS 

CTI3I = USED TO UlMENSIONALWC TIME TO HOURS. 

NTIMC = THE NUMBER OF HOURS TO BE PLOTTED. „ 

P6 S TOTAL SIHULATCO TIME (HOURS!. ** THIS IS NOT RE AG ** 
ALL THE OTHER VAR IAQLtS HAVE BCEN OLSCRIBEO IN THE 
USER' MANUAL • 


51 


99 


-S 


PLOTS U ANO V ON CONSTANT DEPTH SECTIONS I FEBRUARY 1979 MISSION! 
PARAMETER IN r I 7 , JN = 20 , IWN = 16 . JWNS 19 . KN = 5 . ft NH I =4 
DIMENSION UIIN.iN.KNI . VI IN. JN.KNI .OllN, JN.KNI . C I IN f JN , KN I , 
CWHUWN.JMW.ltW) .WtlN.JN.KNI.MRtlN, JN.KN 1 .HRMIIWN, JMH.KN!, 

CKNI.Tl IN. JN.KNI 

DIMENSION 1 &UFI 1000 I .. „ 

READ 51 , PI t P 2 ,Pl 0 .P 4 ,P 5 ,CTT 0 T,NTIME 
FORMAT! I 

uscale sio.o 

VSCALE = 10.0 
NT INC = 3 
ARM IN - 0 . OH 
ARn a x = 0.15 
DO 11 1 NEWS 1 . NT 1 ME 

COY ,D? .ot.VAIJX.fAUY t W.WR , «£5H .TAiItAH.TAv.AhT^ 

CCA.C..A.S . C.CuL.T.TWiRO.ROW.TE.RRCF.TREF , TO, 

'jswisir 1 

DO 99 1 = 1, IN 
DO 99 J=1 , UN 
HJU.Jin.O 
CONTINUE _ „ 

CALL PL 0 TStIBUr. 1000 .il> 

CALL PLOT ( 0 . 0 , 2 . D.- 3 I 

eorhat, 1 1 „ 

00 10 K= 1 ,KN 

P§=k 

Pft=CT T OT* T TOT 
call F ACTOR I □ . 25 I 
IFIK.GT.ll GO TO 20 
30 30 1 = 1 • IN 
00 30 4 = 1 » JN 

If IMARII.JI.EO.O! GO TO 35 

A 1 = 11 - 1 ) * I .0 

A j:| j-x 1 * 1.0 

A A I = AJ tU 1 1 , J ,K UUSCAlE 

AAj:AJ*V<I.J,K(*VSCALE 

YW= 0 . 29 S 0 RT ( IAAI-AI )*« 24 ! AAJ-AJI** 2 ,I 

TW=A| 8 A X 1 ( ARN 3 N / 0 . 2 S.AMIN 1 ( TW» ARlftftX/ 0 . 25 1 I 

CALL JWtOHDtAI, AJ, A A I , A A J , Y W« 0 . 0 , 121 

wm 

GO ro 100 
C CltTIMUE 

DEPTH? 1 1 . 3 /KliHl I * (K-I 1 

00 no in , in 

30 48 jrlljN 

IFCHIU jJ.GT.OCPrH, GO TO 45 
CONTINUE 

LD^=( 0 EP f H/HU*?| j)) «KMN 1 
IF! LD 1 .EC . 0 ! GO TO 55 
L 02 =LD 1«1 
LD3=l5i«2 

01 FF=! 0 EPTH-LD 1 * 002 » 
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oifSfE53,?Eif t " NP ° iCKl 

U2:uil JjJlO*' 

U 3 — U 1 1 , J , L 0 3 i 
V1=V<I,J,L011 

!iiv! ( Ijlf-fei! „„ 

iSllilS-uir^SJiiSSSESSW 

CU=U1 

AV='VJ-2*V2*V1 1/(2*002*6021 
BV=(A«W2-3*V1-V3l/(2*DD2I 

cy=«i 

a2=do2«oiff, _ , „ 

U0iPTH=AU*A2*A2*gU*AZ*CU 

Y0EPTH=AV*A2*AZ*BV*AZ*CV 

SO TO *0 

CONTINUE 

AJrDEPTM 


U^£p Th - * iz»A2 t BU*A2*CU 


C v rv 1 1 ( J, 1 > 

V0LPTH=A¥*A2*A2*BV*A2*CV 


CONTINUE 

l*i. 


A J = l J-i 1*1. Q 
A A I z A I •U0EPTH*USCALE 

aaJsaj*vdepth*vscale 

TH,ARH**.'C,25i5 

CALL .A30M0 < A I ,AJ*AA1,AAJ,Y*,0>-Q,12I 

CONTINUE 

CONTINUE 

CONTINUE 

CALL PLOri-l.Qi- 1 . 0,-31 


C THE HEX T H9 LINES ARE FQ* DRAWING THE BOUNDARIES OF THE 
C DOMAIN. THESE LINES MUST BE CHANCED FOR ANY OTHER DOMAIN. 


CALL PLOT'12.0, 1.0.3) 
CALL PLOT'12.0,6.0,21 
CALL PL0TU3. 0,6. 0,21 
CALL PL0T(13»0,8.C,21 
CALL PLOT 117.3.8.0.2) 
CALL PLOT'17. 0,10. 0,21 
CAlL PLOT t I* .3 ;io.c,2i 
CALL PLOT'l". 0,11. 0,21 
CALL PL0TI12. 0,11.0,21 
CALL PL0TI12. 0,13. 0,21 
CALL PLOT ( IH.O ,13.0,21 
CALL PLOTlm. 0,1*. 0.21 
CALL PLOT (16. Q, IN. 0,21 
CALL PLOT'lb. 0,18. 0,21 
CALL PLOT(1<<.0,1S.Q,2I 
CALL PLQT(l>t.g,17.ai21 
CALL PLOT'll. 0,17.0, 21 
CALL PLOTIH. 0,15. 0.21 


CALL PLOT <9.0, lfi.0 ( 2 I 
CALL PLOT'S. 0, 16.0,21 
CALL PLOT If. 0, lb. 0,21 
CALL PLOT'S. 0,17.0,21 
CALL PL0rilO.3,lT.£,2t 
CALL PLOT'IC. 0,1?. 0,21 


l>Lk rLUl IIU.U, 1 V.U,< 

CALL PLOT'S*"* 19.C.2I 
CALL PLOT '1.0,20.0,21 

2tt PL0Tl?:?*,i9:li:l{ 


Call plot tb.c; i7.n;2i 
CALL PLOT'3. 0,17.0, 21 
CALL PLOT (3.0,20.0,21 
CALL PLOT'l. 0,20. 0,21 
CALL PLOT'l. 0,7. 0.25 
CALL PLOT 13.0,7.0,21 
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1U 

Itl 

lb? 

HJ 

lbn 

ID 

1 oS 
lb* 

!?i 

it: 

itj 

IN 

ITS 

1Tb 

i;; 

IT* 

113 

111 

M2 

113 

Mb 

111 

11b 

it. 7 

II* 

l*" 

1*1 

143 

1*3 

m 

Mb 

1*T 

14* 

14* 

2C3 

2C1 

•* P» *4 

5:i 

::b 

2C7 

?ub 

:c* 

:i: 

’ll 

212 

213 

214 

215 
21b 
21T 
21* 
:i* 
:;r, 

•PIT, s 


c«u 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

Ettk 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 

CALL 


PLOT < 1C.Q , 3 ■ 0, 2 I 
PLOTUO. 0,1. 0*21 

FigTlHieal* 0 '** 

PLCTI-MC,-MO t -M 
Plot io. c. Ml ,31 
PLOT lb. 0,1. 1, 21 
PLOT (6. 0,4. 0,2) 
PLOT to. 0,4. 0,2 I 
PlOTIO. 0,1. 1,21 
Plot to. q,o.o,3I 


THt licit T 25 LINES ARC FCR WRITING ThC CAPTIONS OF THE 
PLOTS. ThC SPECiriC USTR HUST SCRUTINIZE THESE LINES 
AND MANE WTCESSAPY CHANGES. 


VELOCITIES AT K = ,0.0,241 

LAKE KEOWEE- IRIGIO-LID MOOfLl.O. 


23 

123 


CALL STHB3LIQiQf0f4i0.14.211HFXG V 

CALL NUMBER 19 94 . ,999 , , Q. 14 .PB , 0 . U, 0 > 

CALL SVHBOLtO.O, 0.3, 0.14, 3bH L 

C0.3CI 
IF IP 6. GC. 25. 01 GO TO 22 
GO TO 23 

CALL STHaOLll.O, 0.0,0. 14, 28HS1HULATI0NS FOR FEB. 2G 1979, 

C ?6 a lS‘I 2 3 

CALL SYMBOLI1.0,0.0,0.M,2BHSIHULATIONS FOR TEft. 27 1 479, 
C3.C,281_ 

CiL’l ‘ SY' ROil 1.5,1. 7,0. I,12hRUNN0: LOO ,0.0,121 

Call nuhberi99 9. ,99 9.,o.i ,p 10 , 0 . 0,01 

CALL SVHB0Ltl.S,1.5,0.1,33HDISCHARGE VELOCITY 
C3 I 
CALL 


call ^S!g?!?:5 a :S:S:8:!:2 p 4ftB!?e«^E temperature: 


10 

u 

S.P LOW 


i* j ( oi jfui iib Tnu 2 awnonui 

CALL SYMBOL 1 4. 2,8. 4,0. 07. 1 HO, 0.0, 1) 

CALL NUMBER! 3.1,8.3,0.1,P2,0.0,*Il 
CALL SYHg0L<l.S,8.1,0,l t 32HlllN0 SPEED IHAX) : 

CALL NUMBER ( 3.8, A. 1,0. MP4, 0.0. *21 
CALL SYMB0L(M5f7.9,3.1,33HCURRCNTl JOCA SSE FLOWI : 

‘c’JlL NUMBER 1 3.1 , 7. 9, 0.1, PS, 0.0. «1) 

CALL SVHB0L(M5 t 7.7,u.l,31HT0TlL SIMULATED TIME t 
CALL NUMBER I 3*8,7.7,Q«l,Pb,0.0,.2l 

flit 5S«*b : l: » : l:iA“: lifsaailt-iaK • 

CALL SYMBOLI 1.S, &. 7,0.1, 22HVCL0C1TY SCALE I CH/SLCI, C.O ,221 
CALL AXlSI4.1,b.7,lH , ♦ 0, 1 .0, 0. . 0 . , 12. » 

CALL SYMB 0L(0*2,.j. 7,0.21, 2H N,4$.Q,2I 
CALL PLOTIia.O,1.2S;-3| 

CONTINUE 

CALL PLOT 10.0, -3.0, -31 
CONTINUE 

CALL PLOT* iO.0,-2.0,-3! 

END 


CM/SCC,0.0,3 
C, 0.0, 291 
M/SEC t 0.0,32 
CM/SEC, 0.0, 3 
HRS, 0,0, 3! 1 
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> c =.* iUHC « U ). ni!W eit * cseArco on ? maY eo *r moot a $ 

t £**.«*,•••.«*»•*»•»•■. a..**.**! **•***••♦♦•..»»*. ».***«»♦*♦••♦ 

i C This PflOSRAH PLOTS T«£ 0 - . VClOCITUS TOR INC REGION 

t THC WIKI* VAAIAOLCS ARC RKaO WITH AN OPEN fOR^T J« 

C PI 


3 

« 

« 

t 

y 

t 

i 

10 

ii 

u 

n 

15 

U 

17 

II 

!! 

21 

£2 

*1 
2 * 
2 b 
2 ’ 
25 
21 
» m 

11 

12 
33 

1 « 

IS 

i* 

IS 

15 

:? 

52 

••1 

55 

i5 

<tl 

57 

5* 

51 

50 

H 

« 5 
55 

< r 

*L 

47 

51 

S» 
t: 
bl 
• 2 
t> I 
a* 
65 
«ii 
NT 

*1 • 
71 

7? 

7a 

75 

7b 

77 

71 


pI i SI iiSiSH u«°t 5 l!s»t. 

Sl c ! JMHtS 1 ,...!*,... ' 

PS s CUNRlNT • 

Ct TOT = usro TO 0 I HE MSIONAL I 2 C • INC TO HOURS * 

Si’TSi.IE’i.WWKS K.HSJuHrs!'T».» ,, »r .»* 

a LL TME 0 THC R VARIABLES HAVE bCCN DC 5 CRIM 0 IN THC 
SERA HAVIML* 


51 


H 

20 

C 


ss 


31 


8 wms > wiftrAjSM ? twias 5 K ; 5 «:^ 5 i 3 * s « inR * N. J# -iw 

C mh7i w« . jwn. n M l : u i I n , j& . nmi . UR i i n . jw, hn . u A h < i * g. J wg. ** » . 


MW I , 


LWNIl *#N > JHIfltll 

CM I I In* JN 1 * HX I in* JN 1.MVI IN, JR) ,HAA|1N, JN I ,NRHI IMN* JMNI 
( .DTHfYj|ON j TM j lJMNj t j^|K(jl,ROlIM, JN.KMJ iP I ft T H 1 IVN * JMN 1 ,ROU 1 IVN , JWM* 

DIMENSION liliF 110001 

USCALC -10*0 

VSCALC = 10*0 

USCALES2Q.0 

HBVCsO.SJ 

NTlHEsJ 

DO II INEM:1 .NUNC 

CALL REAP? |U Jv ,MM ,P1NTH . I « J ,K .1 M . JV. I N, JN.KN , 

cimn.jvn.d.c.hx .hv.hi.naAiNrh.ai .ah.av .Jp.ox, 

COV.Oi.Of ,lAUX,!AUV.U*WR.-RH.f AI , T A h. T A V . AK T , CD , 
gAj^cjJujjT.Ta.RO.ReVl.TE.RREF.TREr.t 0 , 

CALL PL0TSIIBUf.10a0.lll 
ACAB 51.Pl,P2,pio,P5.P5,CTTOT,NTlhC 

ARRXM: ~ — 

FOftHAT 


CM. 


ARRIM=0.D5 
o* IS 
11 

UP 10 1-7.13 
CALL F ACTOR 10. 25 I 
CALL PL0ri0.0»lb.0,«3> 

00 20 J=1«JN 
P|=I 

~ CTTOTPTTOT 

lHARII.JI.LT.lll CO TO 


20 


uu < 

;|ij 

Aar-IK-ll iHIlf.JI 
AAJ=AJ*V II .J.KlaV SCALE 

fi>i: 5 :?UiNii:si 5 t l . c .;;?ii.. „ 


CONTINi/C 
CONTINUE 

draws bottom SURFACE 


nw =2 
go io 


ill 


. JN 

.Jl.Cfl.Ol 


ir imaI. 

HH-ktfl 
I r INN.GT.ll CO TO 
AAJ=I J-l lal.o 

HK - HHiUMKl , 


GO TO J2 
33 


CALL 
SO Ti 
CONT 
AAJ: 


_ > 0.11 

FLCTIAAJ.AAK.21 


rHifeiU 


CALL 

JD=J 
AJD- tJO- 1 1 * 1.0 
CONTINUE 
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io 

•i 

l\ 

i« 

as 

at 

17 

;i 

H 

t? 

4! 

99 

9« 

9* 

91 

9* 

12? 

122 

182 

10S 

IS! 

1C9 

1C9 

US 

\\h 
11! 
1 in 
l is 

1 it 

\\i 
1 1* 
Iau 

i : i 
12? 
123 
1 29 
125 
i;t 

III 

129 

130 

ill 

lit 

13! 

136 

Ul 

139 

192 

1-1 

Ul 

Us 

196 
• PRT,S 


30 MVMi MO0.0.Q.M1 

AAKs-Hlll . J01«K**H1 . 
CILL PLOT UJD,A*K,a) 


tmc meat q7 limes anr row muwino tmc ^IOunoaries or the 

DOMAIN. THESE LINES MUST bE CHANCED TON ANY OTHER 222fl!JZ 


CILL FACTORll.Ol 
CILL PLOT <-0.5,"2.Q,"3) 
CILL PL0Tlb*L.0.a,2) 
CILL PLOT 16.0,1.0,21 
CILL PLOT <0.0.8. J, 2) 

Silt Tt8M:8:"-3?j;!.. 


C THE NEXT 25 LINES. 

C PLOTS. INC SPECIFIC, 

l AND HIKE NECESSARY CHANGES. 

c * ■ 

i 


IRC FOR WRITING ThE CAPTIONS OF THE 
C USER MUST SCRUTINIZE THESE LINES 


CILL SYHB0LI0.0,1.2,0.19, 36H 


VELOCITIES IT 1= .0.0,291 

LAKE KEOUEC-IRIGIO-LIO MODEL 1,0. 


22 

41 

12J 


CO • 36 1 
IF IPt.GE. 25.0100 TO 22 

CALL°S?MB0Ltl.0,0.9,0.19,2«HSIMULlTX0NS FOR FEB. 26 1979, 

c 8*°»SM 23 

CALL SYMBOL 11.0, 0.9, 0.19 ,26H SIMULATIONS FOR FEB. 27 1V7V, 
CO .0.28 1 
continue 

- LOO ,0.0,121 


CALL SYMBOL* 1. S, 9.5,0. 1. 12MRUN NO: 

CALL NUMBE&12. 6,9. 5,0. 1.P10.0. 0,01 
CILL SYMB0LI1.S,9.3,0.1,33H&1SCHARGC VELOCITY : 

C 3 » 

CALL SYMBOL ll.sJvIljSllj^iMDisjHA^GE TEMPERATURE: 

call STMBOLII. S, 8. 9, 0.1, 32 hMI NO SPEED (MAXI t 
CALL NUMBER 1 3. ■ ,8. 9,0. 1,P 9,0. 0.-2 1 

CALL SYMBOL(I.S,i.7,Q.l,3lNCURACNTlJOCASSC FLOW): CM/SEC, 0.0,3 

CALL NUMBER(3.S,8.7,0.1,PS,0.0,»1> 

CALL SYhBolI1.S,8.S,q.1,31HT0TAL SIMULATED TIME : 

CALL NVISEh I3.8,8.S,0.1,P6,0.0,*2> 

CALL SYMBOL (1«S,|.Q,0.1,19HSCALES 1 HORIZONTAL I , 0 . 0. 19 8 
CALL SYMROLIl. £,7.5,0. | ,2 OHLC NG Th SC ALE I HE TEAS) ,0.0,201 
CALL I VI SIN. 1,7. 5, IN ,«U. 1.0,0. ,0, ,610. I 

CALL SYNB0Lll.S,7.a,C.1.22HVCL0tlfv SCALE 1 CM/SCC I , 0.0 , 221 

r 


CM/SEC, 0.0, 3 
C, 0.0, 291 
M/src,0.a,32 


HRS, 0.0 ( ID 


ALL AXn<9.1,7.0,lH ,«0, 1 *0,0. fO.,12.1 

eltt T!S8ti|:!:l:8:8:i:isim?& lK!HS{W«8i!s!i 

CILL lxr,S | 9.I,6.0,lH t »0.1.0.n..0,,2b.l 

.5,S.5,0.1,22HVClOCITV SCALE I CM/SEC I ,0.0,221 
*0, 1 .0,3. ,0.,6. 1 


L 1 71 
. 0 , 20 ) 


10 

II 


S.PtV. 


CILL SYMftOUlc . 

CILL AVI., 'i, 1 . S.S< 1H . 

cill plot Io. 3, o.5, -3) 
csNriA »r 
CCNl lh\J * 

END 
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•tCcV*i.UMOMUi .nyW CLT UMiS on 2 AUG, 1 0 AT" 


this, program plots rHt v - w velocities for thc region 
or INICRCST. . . . 

rwr following variables arc rcao with an open format i- 
pi : El5C«A»tl ^MPtRAtiRC. 

PlCr RUM NUMBER. 

PH = WIND SPEED (MAXIMUM), 
pi : CURRENT. 

CIIOT S USED TO 0IMFNSI0NALI2C TIME 10 HOURS. 

UTIME = TOTAL NUMBER OF HOURS SIMULATED. 

Pfc : TOTAL SIMULATED TIME (HOURS). ** THIS IS NOT RCAO *• 
ALL THC OTHER VARIABLES HAVE BEEN DESCR1BC0 IN THC 
USERS MANUAL. 


PARAMETER IN=17 fJN=20.IWNS)6 , JWN = 19 . KN = S , KNMl =H . 

DIMENSION Ul IN , JN.KNl.il I N,JN,KN» ,01 1 N, JN , KN I , C I IN > JN r RN 1 , 
CwWUWNfJUNtKNlJPINTHliWN.JttNI , 

DIMENSION MX IIM, JNI.HrllM, JNI .HTllN.JNlfHARIINfJNI ,HRH«IUN,JWNI, 
CWIINtJNiRNI .W*l I Ni jn .KM 1 »¥(?H < IN • JN.mN > 

DIMENSION T l IN • JH *pNi ,R0 5 IN. JN , KN) , TNI JWN . JNM, 
CKN1.R0WIIUH.JMN.KNI 

dimension laufiioooi 


USCALCslQ, Q 
VSCALEslfl.O 
W1CALE:I0.0 
N0TL =C .33 
NTIMCrl 


DO II INEWsl.NlIHE 

CALL AEAD?iU;v .UH.PINTH.I , J, K . IH , JW. 1 N . JN.KH. 

C I Ml. , JWf), 0 .t .NX (HV.HI.HAh.mRH.AI . AH. A V . AP .ox. 

COT lOl.OT. i AUK. iAUTfW.MII ,VRM,T A I , TAH.TAV ,AK T,CB ,CM. 
CT i^l^jT.TWfRO.ROW.TE.RRCF.TRSr.fo. 


CALL PLOTSlIBUF. 1000, III 

READ 51 .Pt ,P2, PIO.PH.PS ,CTTOT,NTiHE 

FORMAT! I 

ASH1N=0.0§ 

ARMAXsO.lS 
FCRHAT 11 _ 

DO 10 J=1 
CALL FACTOR IQ.2S I 
CALL PLOTIC. 0,16. 0,-3) 
gO. 20 1=1, IN 

P*=CT10T*TT0T 

IF (MARIIfJI.LT.il) 00 TO 20 

Airjl-i 1*1.0 

DO 30 n =1 ,kNH1 

AMs- I m - 1 1 •Hill ,J A 

A A I :AI •0(1 , J ,M t •USCALf 

wu,i,kjio!a’ 

AAK = A«C-Wl F> J,K 1*MSCALE»HSYL 
TW=0.2*59RT( IAAI-AI1**2*(A Ad-All 1 »«2) 
l! SI * ARNIN/Q. 25 .AMINl (TM.ARMAX/0.25M 

CONTINUE ^ I , A!l»iAltAAK,Ty,0.0,121 
continue 

DRAWS BOTTOM SURFACE 
NN=0 

QQ IS1 ,IN 

IF IMARlI.JI.ro. 01 GO TO i. a 
NM=NM*1 

IF (fM.GT.D GO TO 33 
A A I s 1 1 -1 1 • 1 .0 


CONTI * m 
*»I=(I-l 1*1.0 

vstitai itwm.*, 

10=1 
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2? 

tt 

II 

IS 

i\ 

ci 

tt 

S' 

tt 

■iJ 

St 

V5 

St 

97 

M 

133 

1 Ll 

IB 

ICS 

let 

1C7 

ICS 

IC« 

113 

III 
113 
III 

11 ? 

lit 

117 

111 

l i! 

131 

13? 

133 

13*. 

12! 

lit 

131 
1 29 
133 
131 
133 
133 
] 3* 

1 35 
lit 
137 
131 

US 

1*1 

It? 

Its 

1st 

1 S5 
lt& 
1ST 

Its 

.P«T,S 5 


. AiOi^iOrOti.O 
1? 


a A* 

CALL 


0 ,J< 


: -h I < I .jlvKini 
L PLOT IAID, AAH, 3 > 


C THE NEXT 07 LINES ARE POP DRAWING THI g C ON n A « X p 5 „0 f* ^ ThC . 

C COPAIN. ThCSC LINES MUST 6 C CHANGED FJH ANY OTHER DOMAIN. 


CALL rACTORll.OI 
CALL PLOTj-O.Si- 2 . 
CALL PLCTIfc,- 1 ~ 


•5*“2.0.»1« 

0,6.6,51 

0,1.0,21 


CALL PLOT ( 6 . 0 , 

CALL PLOTtO. 0 , 8 . 0,21 
CALL PLOTIC. 0 , 0 . 0,21 
CALL PL 0 T(Q. 0 ,- 2 . 0 ,-J» 


C TiC VCXT 25 LINES ARC FOR WRITING THE CAPTIONS OF THE 
C PLOTS. THE SPECIFIC USCR MUST SCRUTINIZE THESE LINES 
C AND HARE NECESSARY CHANGES. 

C , 


VELOCITIES AT J= , 0 .tl, 3 tl 
LAKE KEOUCE-IRIGIO-LIO MOOCLI.O. 


2 T 
125 


CALL 5YMB0LIQ.0, 1.5,0.1s, 2SHfIG 
CALL NUHBER(5.S,1.5,0.1S,P8,0.0,OI 
CALL SVHBOL(0.O,1.2,O.IS,36H 
CO, 36 | 

IF IPi.CE.2S.QIG0 TO 2? 

GO TO 23 

CALL SYHBOLI1.0,O.V,C.lt,2tHSlHULAUONS FOR FEB. 28 1979, 

{23 

CALL SVHBOLll .0 , 0 .9,0 . 1 t , 28H5I HULA T I ONS FOR FEB. 27 1979, 

C Co2f INUE 

CALL SYHP0Lil.S,9.5,0.1,12MRUN NO: LOO ,0.0,121 
CALL NUMBER 12.6,9.5,0.1, PI 0.0.0 ,01 
^C ALL SYH6OL<1.S,9.3,0.1,33 HuI SCHaRGE VELOCITY : CM/SEC, 0.0, 3 

LL NUHBFRI X. «. 9 . 1 . 0 . I .PI . n.fl. .91 

C, 0.0, 291 
H/SEC,(J.D ,32 
cn/scc ,0.0,3 

HRS, 0.0, 311 


CALL NUHBER« 3 . 8 , 9 . 3 ,Q. 1 ,P 1 , 0 . 0. *21 

CALL SYHB 0 Lil.S, 9 . 1 , 0 . 1 . 29 H 01 SCHARGE TEMPERATURE: 

CALL STHBOLIt. 2 , 9 . 2 , 0 . 07 , 1 H 0 , 0 . 0 . 1 ) 


C) 


SaLL SY 5 tS?Shl!i*l 5 Blll 5 lH 5 i «5 + I^EEO 


(MAXI 


JOCASSCC FLOW: 


CALL NUMBER I 3 . 8 , 8 . 9 , 0 . 1 , PS, 0 . 0. *21 
f CALL SYMB 0 Lll. 5 , 8 . 7 , 0 . 1 , 33 hCURRENT 

CALL NuMBERI 3 . 8 ,B. 7 , 0.1 ,(* 5 , 0 . 0 , * 1 1 
CALL SYMB 0 L{l.S,S.S, 0 . 1 , 3 lMT 0 TAL SIMULATED TIME 
CALL NUMBER ( 3 .B, 8 . 5 , 0 . 1 ,P 6 . 0 .O, * 2 ) 

CALL SYMS0L(1.5,8.Q,0.1,19HSCALrS (HGRJ ZON T AL 1 , 0 . 0 , 1 9 1 
CALL SYMBOL 11.6,7.5,0. 1 ,20HLE NG TH SCALE (ME TERS I ,0 . C, 20 ) 


CALL AXXSIS.1,7.5,1H , *6, 1 .C . 0 . . L . , 61 0. 1 

CALL SY MB OLll.S, 7. 0,0. 1.22H VELOCITY SCALE I CM/SEC 1 , 0. 0 ,221 
CALL AXISIS.1,7.0,1H .*0. 1.0,0. ,0.. 12. I 
CALL SYMB0L(1.S,6.S,0.1,17 HSCALeS l VE RT ICA L I ,0 .0,1 71 
CALL SYMBOL!}. S, 6. 0,0.1j20HLtNGfh SCALE f ■ * - 
CALL AxIS(S.I,6.0,Ih , *0. 1 .0.0. .0. ,20. > 

CALL SYMBOL! 1.5, £.5.0.1 «22HvE L0C1T V SCAt 
CALL AXTSIS.1 .S.S.1H . *0 , 1 .0, 0. ,0. ,6. I 


SCALCfCH/ SEC), 0 . 0 , 22 ) 


I? 

fcCHKorl 


CALL PLOTI 10 . 0 , 0 . 0 , - 3 ) 
T Nu 
I NU 


continue 

CONTINUE 


END 
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„Cc.*|UMWO) 


j 

* 

3 

* 

4 

« 

13 

w 

r 

I* 


la 


ii 


is 

24 

27 




3'. 

U 


!! 

R 

;a 

39 

:i 

•2 

93 

44 

43 

«4 

47 

•*# 

S3 

31 


II 


34 

33 

34 
37 
33 
54 
*•' 
ol 
• 2 
41 
4 ■* 
45 
44 
4? 
44 
44 


73 

74 

73 

?* 

7a 


• E<W0K a, c.e.nt. « e,V*N|o 

,B,! ,s t «SIcSS3*SB'A?tcS 0 T,|!S e j.f8¥. r ?S5** <0 P,,K '"‘" 


T NC COMPLETE P ACKAGE COMSltTS 
ALL 3 ARE CAM LOGUCO TOGEIhCR 
AND DECKS ARC NOT NEEDED. 


?C f«l u fiS 0 ^Zi7Lt’..£gl’8°KA 0 Sfg L iSAc 


ANO END 
C CtIKOA 


ARY RECTANGULAR GRIMED SCALAAFlCLDCAN fC^CCMTOU^CD QNHIIGO 

or calcohp type pi,pmR_Bv 55 .ttin6 qp_pkoper .calling ^Rcuhents and 


PROCEDURES AS INDICATED BELOW ANO THEN CALLING tCHKON. 


• CALLING S T ATEMCNT IS AS fOLLOWS- 


CALL ECMKONCHH,lN!,IN2,N£Xl,Nrx2,MfYl ,NrV2,HI,WIO,ptTlWC,SAHCOH, 
aCONlNT.RGRXD.lN 3 , lN4.ZLIT.aBIG, ANOA TH , A30U * At AST , AWE ST i NOA 5H0 . 
3 NO A SHU iXCABCLiSHOOTH.IRCCcYJ 


—DESCRIPTION OF CALLING ARGUMENTS— 


HH 


GAINING GRID DATA TO BE CONTOURED. ITS 01 flFNS IONS 
]N2. DIMENSION HMLIN1.IN2D. POINT 1>| JS .LOWER LC 
HD. INI IS DIMENSION IN X DIRECTION *AND IN2 IS 


left 


IS ARRAY CONTAINING GRID DATA TO . 

ARC INI ANO ]N2. DIMENSION HHL IN 1 . 

CORNER OF GRID. INI IS DIMENSION 
0 1 MENS ICN IN Y DIRECTION. 

CX INCREASES FROM WEST TO EAST AND T INCREASES FROM SOUTH TO NOR TH3 


NCXl. NCX2 1 NEY1, ANO NET 2 DETERMINE THE PORTION OF HH GRID TO 

BE USEO. NCXl AND NEX2 ARC THE FIRSTCLEFTHOSTl AND LASTCR IGhTHOSTJ 


COLUMNS TO BE USEO. NET 1 *«D NEY2. ARE-ThE F IRSTl BOT TOM 3 AND LAST£TOI‘ 
nQWS TO b£ USED. ETHUS any SECTION OF HH CAN BE USED! 

FOR FULL GRID— - 
. NEX1 > 1 

I { Nl 

NEY2 > IN2 


Nl IS HEIGHT IN INCHES OF CONTOUH MAP BETWEEN LIMITS NIY1 ANO NCVZ 
u 10 1:1 WIDTH IN INCHES OF CONTOUR HAP BETWEEN LJHITS NCXl AND NEX2 


PL T INC IS STRAIGHT LINE PLOT INCRCHCNT IN INCHES TO BE USEO 

ALONG CONTOUR. GOOD VALUE IS .04. Buf CAN BE VARIEO UP OR OOWN. 
SINCE LARGER VALUES CAUSE PROGRAM to RUN A LITTLE FASTER. IDEAL VALUE 
. IS LARGEST THAT HILL STILL GIVE SMOOTH LOOKING CURVES. 

. — DO SONC EXPERIMENTING WITH IT. START WITH .03 OR .04 ANO INCREASE. 


SAHCON IS ANY SAMPLE CONTOUR VaLut. rT TS USED AS A STARTING POINT 
FOR COUNTING UP ANO OOWN TO GET OtHEfi CONTOUR VALUES. 


CONINT IS CONTOUR INTERVAL TO BE USEO. 


RGRIO IS AN INTEGER*? STORAGE *RR*Y USED JNTERnALLV IN PROGRAM 

ANO NEED NOT BE INITIALLED. IT IS INCLUDED AS AAGUAENt TN ORDER 


TO TAKE ADVANTAGE OF VARIABLC DIMENSIONS. DECLARE AS INTCGCH*2 
BEFORE CALLING. 


INJ ANO INN ARC X ANO Y DIMENSIONS OF RGRIO. DIMENSION RGfllOC IN3 , IN4 3 
IN3 MUST BE AT LEAST AS LARGE AS NCX2-NEXK I 
I AH MUST BE AT LEAST AS LA»GF AS Nf Y?-NEVl< 1 

ETHUS RGRIO MUST BE AS LARGE AS PORTION Or DATA ARRAY HH BEINB USED 


CUT ANDJMIG ARE LOWER ANO UPPER CONTOUR CHICK LIMITS. NO CONTOUR 

Be drawn below value or zlit or above value or zoic. 

UVENT DRAWING FOR ANY COMPLETELY WILD OATAJ 


I « HW 4 


UL TO PRI 


ANCRTN, ASOUTH. AE AST . AND IWFST CAN BE USED TO ELIMINATE ANY 

nundER or inches from any side or final orawing. 


FOR 


WnkiBMKTa 1 JSSSYMtiKtW? ¥i D Ig s I. w, °- 


FCR EACH OF THE ABOVE *#Th POSITIVE VALUE. THIS HANV INCHES 
WILL *6C ELIMINATED ON SIDE TO WHICH IT APPLIES. 


THIS allows US TO riT. ANY RECTANGULAR GRID TO ANY MERCATOR 
0* OTHER MAP LIMITS WITHOUT ACTUALLY ADJUSTING THE GRIO. 


NOASHO HID NOaSHU CONTROL TYPE OF CONTOURS (SOLID OR DASHED L1NCS3 
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7- 

AS 

11 

12 

13 

M 

IS 

IS 

13 

If 

93 

w 

u 

95 

9k 

97 

**2 

99 

1 u3 
1C1 

I S 3 

ics 

l =«. 
137 
133 
139 

li? 
112 
I 13 

1 19 

IK 

117 

U« 

1 19 

123 

121 

122 

12? 

129 

125 

1 2’ 
129 
129 
133 

tt\ 

1 53 
1 39 
13; 

13b 
1 J7 

III 

1 93 
1-1 

i:i 

1*9 

1*5 

1st 

19’ 
1SJ 
199 
I 5” 
:ti 

152 

153 
159 

*I 5 

15b 

147 


C 

t 

C 

C 

c 

c 

III 

r, 

L 

C 


IP £ 1THEE Of. 0OTM Ale Zt F0 O-l L£SS CONroUlCS A^E-SOUO LlNT£> 

IF BOTH ARE POSITIVE, CONTOURS WILL BE OASHEO AS FOLLOWS—' — 

PEW DOvM SF.CTIOM LENGTH > N0AS«O*PL1 INC CPLUnC IS INCREMENT LEMGI 
PEN UP SECTION LENGTH > N0ASHU9PITINC ...... 

C THUS LENGTH 07 DASHES AND SKIPS IS FUILV VARIABLE] 


XLABCL CONTROLS LABELING Of CONTOURS . L lUES ... 

ONLY IP XLABCL GREATER THAN ZERO. VALUE Or XLABCL. 

IS HEIGHT IN INCHES OF LABEL NUMBERS. LINES ARC LABELED 


ARE LABELED 
Of.XLABCL 


’with nearest whole number value or contour, if special 

hfESklCS iS 5SSSU26.°Shl»S5?i.rci:SSS c 8uS"r 1? IBSiioto. 

shoufw is a control roR varyinc contour smoothing. 

*“ VALUE BT TWLEN 0.25 AND T.5 


1NZTIALI2E SMOOTH TO SOME 

CANY VALUE OUTSIDE THIS RANGE IS SET I NTCRNALL Y TO 1 .QJ.. . „ 
LARGER VALUES GIVE SMOOTHER CHART WITH LESS DETAIL. WHILE 
. _ smaller values give LESS, smoothing and more detail. 

C NORMAL VALUE r*R MOST SIMMS SWUD BE ABOUT 1.S3 , .. 

ANYTHING LESS THAN AOOUf 0.90 OR LARGER THAN ABOUT 3* IS . 

PROBABLY NO GOOD. BEGIN WITH 1.5 AND EXPERIMENT UP OR DOWN 
TO DETERMINE MOST DESIRABLE VALUE FOP TOUR NEEDS. „„ 

CINPUT GRID OATA VALUES ARE NOT ALTERED IN THIS SMOOTHING! 

. IRECCY IS :*LOT TAPC RECORD COUNTER. 1NITXACI2C TO NUMbER 

OF PLOT RECORDS WRITTEN BEFORE FIRST CALL TO CONTOUR SUBROUTINE. 

all OF THE ABOVE ARGUMENTS EXCEPT ARRAY RGRIO MUSI BE OEFINEO. 
ARGUMENTS ARE NOT ALTERED WITHIN PROGRAM, ANO RETURN INTACT. 

PLOTTER BUFFER SPACE MUST BE SET UP AND CALL TO PLOTS 
MADE BEFORE FIRST CALL TO THIS SUBROUTINE. 

PLOT TAPC MUST BE CLOSED OUT AFTCR FINAL CALL. 

ANY HUMBER OF SUCCESSIVE CALLS CAN BE MADE TQ CONTqUB ' 

SUBROUTINE ECHKON. EACH MAP BECOMES A SEPARATE PLOT RrCORO. 

NO INTERNAL MAP SPACING IS PROVIDED. WITH PEN RETURNING TO 
,, ORIGINAL OR IGINCLOuCR LEFT CORNER] AT COMPLETION OF MAP. 

CTWUS IT IS SIMPLF TO PUT MORE THAN ONE SET OF CONTOURS OM SAME MAp] 
aWYSPTCIAL MARKINGS OR LABELS THAT ARE OCSIREO MUST BE DONE 

this subroutine oraws contours only 

With INCOMING ORIGIN BEING LOWER LEFT CORflER OF CONTOUR CHART. 


SUBROUTINE EChKCN I HH . IN 1 , 1 N2, NC X 1 .NEX2.NEY 1 ,NC Y2.H1 ,W I D.PL T IWC. 

3 S i? C 2S , f 1 2? i ^* B ^{9»fy 3 r w ^»* LIT l* BI6 »* N0RTH * A!i0uTH » AE A£ T » 4tfES f * 

JNOAsMD.WDASHU.XLACEL, SJlioTM.IRECCYl 

SEE ABOVE COMMENTS TOR DESCRIPTION ANO USE OF ABOVE ARGUMENTS 
^®J9N0ly y S TX7C ON ✓ S MM Ij.SMWI.X.Y.XGR10,YGRID.CUTOF.SOH1.SoWI.TMAX.xPP, 

2.Y pp.cgig.u.v.nxux.jooo.nuvx.nuvy.yortk. SOUTH .CAST. West. CL IT. CBIG. 

rLCLX.LCLY.iNdROS,0lNC,dLOSIt,PVAL.PVOL;NENTE4.HlNUM,NMXl,NMYl, ’ 

*d5la 5?*<Ju I * VALLIN.HlNC. HAXCRD.WHAT .L0ASH1 > LOA SH2 . D AS HER , 

"COMMON '/QENDEQ/HI XCM i MODE jHOCh . Xitl, AST , Y YL AST 
LOGICAL DASHER ,U0LAB5, OUTS 

OIMCNSION HHIIN1.IN2I , RCR1DI I NJ, 1H9 1 » »£N < 9 1 
INTEGER RGRIO 
DATA IMAP/O/ 

MMC SUBROUTINE FOR CONTOURING .SCALAR FIELD ON CALCO^P OR, MILGO TYPE PL0‘ 
I:1M4P;1 


VRITCft.lO)! 

10 .2X.11NCONIOUR MAP, IU, iX,2iHIMC0MlNG XRfiuMeMTS rOLLONl. 

WHifClJ.ianMI ,lH2 ,MCX 1 . HE JT2.NE Y 1 , HEY 2 .HI , W10 *PC TTMC , SAM CON, Cum 



1 

HT 


WT*2 Ht Y1 NEY2 HI WIO PLTIffC SAMCON CO 
9.3F 12.?., 21101 

ASOUTH, AEAST, AWES! , NO A7M0 , NO A SHU , XL ABE 

- AHORTM ASOUTH AEAST AWES’/ NDASHD NOASMU X 

2CABCL SMOOTH IRECCY ,/ 1 ,2X ,Zf I 2. 3, 9F 10.9 ,2U, 2f 10. J, lb, / 1 
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in 

up 

lil 

147 

til 

til 

UY 

tit 

it* 

173 

ill 

17* 

175 

174 

1 77 
17t 
17* 

Iff 

lil 

113 

16* 

liS 

m 

14S 

14* 

1*3 

1*1 

i*r 

HI 

1*6 

197 

l*a 

t** 

233 

731 

$s; 

2 l* 

2 3 r 
2i.4 
237 
231 
2 3* 
21 ” 
211 
212 
?13 

111 

214 
217 

2 1 5 
21 * 
*. iU 
221 
22 ? 
221 
22 * 

225 

226 
22 ’ 
223 

! !* 
2l! 
2 32 
223 
2 2* 
2 3*. 
2 24 


BttfKB. 

coninciconimi* ^ 
HjCONlMC,,NC,Q.)<iO TO 3 
white < 4(211 ; . 

F0flMAT»/.2X 
60 TO 1X9 
ItCl! - - 


,3H/1AP.13|tfH UKO IMrCHtfAU 
7*1 


lBaS^ 2CN^A5HU 
?nLD5i»Ja.OT < O.AM).IOASHl .«r.0> DASHERS. TRUE, 

WUffiflftt:,... 

IF I CONI tic . L I . 0 . > C ON I NC - -C 0 H 1 NC 

ftb^?8 5 * C0NINC 


30 


32 

« 


VALUN = -*t*t9t,*I9 

WUVX:NEX2-h|xj*l 

!l^]fMu9x!GT'.3!ANO.NUVX*LE,IN3.4ND.NUVY.Ct.J.ANO(.NUVV.LC.lN4)r.O 
WHnrU,7lNEXl,NEX2,NUVX,WCYl»NEir2.NUVV . . , 

FORMAT t/. 10X, 2 3H&A0 ARRAY LIMITS. S« IP . / 1 QX , 3 110/ 10X i 3 1 10 1 

co to tin 

SKIP IF mm OR NUVY CESS THAN 4 
YORTHIHl "A NORTH 

Yr (wcst“t.o.>we.«tso. 

ITlEAST.OT.MIOlEASTrWlD 

IfiSOuTHsLT.Oi !S Out MIC. 

ir I vor th.ot .hi iycrthihi 
hixershinuh 

muMW-iaiH 

XXLASTS99. 

YVLAST=99. 


10 ft 


Q INC iPLTXNC 
CLlTSOINC/I , 
CBICSQlNC/2. 


9* 


TMAXI4.CU1 YORTH-SOUTM» EAST -WEST I 
A6RTO=WIO/FLOATINUVX-1) 

VC«IO=MI/FL0AT l|fUVY-ll 

hircsxcAio , 

IF|YCRlo.LT.XCRJO)HXNC:yC«lD 
X: SMOOTH 

IF1X.LT..Z6.0R.X.GT.7.5IXS1.Q 
h INC IX »HIfiC 

CUTOFISOR T (XGRI06XGRI0 *Y6RX0*YGRlpi * .01 
CLOSITl.04 

CL0SI7 IS VALUE FOR CLOSED CONTOUR CHECK 
nmxi =»<*:*, i 
NNY i liliE Yi 
NHXlUAIftXl-1 

MHIrMUI* 

reySincyI-I 


NEXT DETERMINE MAX 
ZMAXiMMlNEXl.NCYll 
2HIM-2HAX 
DO 30 1=A1CX1 i NCx 2 
^ 30 J=NCYl,HtY2 


AND MIN V* LUC S IN SCALAR FIELD 


lHHII.JI.6f.8MAX IZMXK IHH 1 1 * J1 
LI*U».I.T.*HXN» XHrKSHt* ( I » J1 
coat i Hut 


NEXT Of TERM INC SOTTOH SfARUWC VALUE FOR CONTOUR LOOP 

WftSftm.J.IGO .0 „ 

P.it IPV*L*C0NINC 
GO TO 32 

IFlPVAL-CONINC.LT.ZMlNKiO TO 35 
PV4L IPVAL-CONXNC 
CD TO 34 


110 


?:? 

242 

2h3 

?•*** 

ft! 

24’ 

lit 

249 

- I * 
«■ 0 

Hi 

2 54 
251 

IS! 

css 

$tl 

2e2 

263 
249 
2 dS 
244 
247 

St« 

264 
27H 

271 
1 11 

2 7? 
2 74 
2 75 
276 
2 77 
278 

2 74 
?«3 
281 

2a j 
2H 

Si 5 
246 
2*7 

ill 

29" 

2«i 

292 

553 

295 

246 

297 

249 

249 

*33 

ICt 

IC2 

JGJ 

3C4 

3GT 

*;< 

Sc* 

313 

111 

31 ’ 

313 

314 

3 IS 


M2sN0Vf-l 

CONTOUR LOOP*SIABtS BcLOV hT JCTAriMtA/f 34 _ „ 

This LOOP DETERMINES WH£Rt TO 5TARf * NCfc> CONTOUR, TMfN CALLS 
SUBROUTINE CONCIM TO ORA* EACH CONTOUR , EXIT IS MADC WHEN 
AU CONTOURS CCNpCCTEO. 

WERE ARC 2 SCANS FOR EACH CONTOUR VALUE. flRM WITH VARIABLE OUTS AS 


FALSE SCUCTS ONLY CONTOURS ENTERING CRIO rROH OUTSIDE (ROES. 

SECOND SCAN WI Th OUTS TRUE SELECTS REMAINING 1 NNCR CONTOURS, 
STARTING POINT CLOSEST TO PLOT PEN POSITION IS SELECTED IN EA 


EACH CASE. 


36 


600 


IF IPVAL .G£ . 2 HAXIG 0 TO 110 
OUTSs.FALSC. 

00 37 |=I, N 1 
00 37 J=l,N 2 

JSSMW t0 

DO 103 1 = 1 , N 1 

?? < i8?S Voi *¥S 600 

IFlT.£0.1,OR,J.£Q.1.0R.I,tQ.Nl,Ort.J,EO.N2tGO TO 600 
GO. TO 100 

IF (RGRIOIX ,dt.EQ.l»GO TO 100 

IF(RGR1D(I,JI.GT.1.AND.0UTS|G0 TO 100 

II=NEX4*I 

JJSRCY4* J 

MEW 1 > =HM( II ,JJ) 

HCN(2»=HHUI,JJ«1» 

HEN< 3 »=HH( 1 I 41 ,JJ* 1 > 

NCNj 4 t=HHIlI*|,JU) * 

Do *, /i C K = I ,4 

- ------ “- 400 


TflSfiiK 5 ?l 8 !rtIfMeS 5 s!?fttt 5 f{ 8 i' 

If (HE NIKI .U.PVAUHEN(K)=PVAL-PVOL 


4 DO CONTINUE 

IF (OUT S ICO TO 2 S 0 
NENN =1 

IFtI»C<|,l.AN 0 ,HENtl).GT.PVAL.AND,tiENI 21 .I.T.PVAUG 0 TO 601 
NCNN-3 

irU.£C«NX.ANO.HENI 3 I.GT,PVAc.AND.MCNI 41 .ET.PVAL)GO TO 601 

if IJ^C.l»AND.hEN(4).GT.PVAL.AN0.HENtll o Lr,PvAUG0 TO 601 
NENN =2 

i? l ^ t ?s!! 2 * A,, 0 * HEN, 2 »* CT * , ’ w *E.ANO.rtEN( 3 ).ET,PVAL>CO TO 601 
GO TO 602 
DO 410 K=l ,4 
TlsK 
l2=K*l 

IF<K.EQ .4 112 = 1 

IF(HEM(in.GT.PVAL.AN 0 .HCN(I 2 ).LT.PVAUG 0 TO 400 
CONTINUE 

Mr 602 

ZF|fi«RIO|I,JI.£O.O.Oft,OUTS>GO TO 640 
I1=RGR10|! ,JI/10 
I2=PCRI0(I i Jl-10*Il 

IF IT 1 •£ 0 «AICNN*QR ,12 •£ C ,NE NNtCO TO 100 

r = Y«I^«lFLOATU-l»HPVAL-MCNll H / IHCNI 21 -MEN 11 1 1 > 
Xr«GRIO*rLOAT( I-ll 
60 TO 45 

Xg^glg^rLOAr^X-O^PW/IU-HEM.ail/.HEN.SA-nCNiai.l 
CO TO 15 

V=YGRID*FEOA T f J-l I 

d=ix-xpp»»ix-xpp >♦ t r-rPP)«iY-rPPi 
iriO.GE.02»G0 TO 100 
D*=D 


250 


4140 

405 

601 


6«0 

602 

340 


440 

3¥l 

SUi 

15 


51 ? 

ill 

hi 

m 

m 

3 21 
32 * 

m 

) *j« 


131 

1 31 
33* 

1:1 

St* 

!<*5 

344 

$5$ 

349 

Say 

3S1 

J?I 


raiif 8 " 1 "" 

fill ?' 

!ff^eRioacL^iiuiJq.o?^otiro(LCL^i.cLV)»i 

urn 


WHITE lL, 13 nPVAL,pat.A 6 S, 0 UTS 

roBrt NCxT X int*li5aRiuliNE conlin to actually draw contour with value pv 

CALL CONL!NIHH,lNl,IM 2 ,nGRXD,IN 3 tlM) 

NOV GO BACK TO INNER LOOP TO SEC If THERE ARE OTHER PVAL CONTOURS 
TO BE DRAWN. 


IDS ^?<c|t$JOO TO 612 


OUTS ... 

QO To.38 


4 TRUE. 


«12 PVALSPVAt ♦CONI NC 
INCREMENT C“" 
Go TO 36 


ontour and co to top or loop r or next contour 


110 CALL PLOT (0, ,0 * ,~3l 
I M AP S I HAP4 1 
7RECCY = TRE CCY ♦ 1 

Write ig»* im imap , ircci .ijjeccy 

“^SnS^eId^’Wi} 1 | H f OMTOU ” HAP,X3,2KH BEGINS WITH PLOT RECORD, I 3, 1«|W 
WRITE! &il 16 >m6|?NC,VALLIN,HAXCR0.WHAT 

i i a 0*^4. a • < a kiuua** « fur f kif nc up it f f i r t mu aa. aaaJ* Aim 

t WRnAi i 1 1 a (6 Annuo i uin niwnLiiuiN tutun wn wVfiiWwrr 

12Q 2 a1?HH9 st SCUARCS , I* , 12H ON CONTOUR ,F 10.2,/i 


»P ft! S.CONLXV 


' iH **" 


112 . 


-•SuMJt.tl » 
\ 


i 

9 

15 

i! 

H 

l* 

16 
17 
la 


?* 
• r 

zi 

27 

:s 

33 

31 

3? 

33 

l« 

13 

it 

17 

33 

39 

Hi 

42 

m 

•**» 


*b 

47 

•*9 

t s 

f T 

C 1 

n 

:4 

55 

56 

57 
51 
£9 
t" 
bl 
ts 

63 

tl 

64 

b 7 
ta 

69 
^T. 
7 S 
77 
77 

74 

75 
» ( 

77 

7a 


oash£r* 


,<Oi*UMHT C|?6Artb OK 2 ivu 74 A.r ir.li.2R 

:M1A WhRMr V £ ««««^ • 

ZYPP.CGlC. J .V'MKJA'JDOOttfWVXiNUVV .fORTM, SOUTHEAST .VEST .CL IT. CD 10 , 
JtCL X .LCL Y . i 
4 WHK 1 liMliyil . 

*p?»ic8£loN r ^MlIll! .1M21 ,ClOCt^b2> ,xxPi:0T(275,ZI,HAXI*<i,LC)<Em . 
2(03014 00.2 I, Hi P0 SI 400 7 
INTEGER RSfilbf 1W3 .IH4 I 

LOGICAL INCS .DOLANS .OA5MCR , CLOS ( OUTS ( DASH 1 X 

THIS SUBROUTINE 15 CALLED TO ORA* EACH INDIVIDUAL COM! OUR 

17 DOLAftS ENTERS AS TRUE. LA0CL CONTOURS fcJlTU HEIGHT HI HUM 
DAJMlV=DAJWr« 

LASLIT:? 

IF (DOLAAS1LA0L it=o 
Jncs:.ML 


use. 


1DPLOT si 

«HABD:LDASM1 

NSOFT = LOA Sh2 

NUOCrO 

xxsx 

rr=r 

xeiqrxx 

TBlGiVY 

LUSICLX 

L2Y=LCLY 

1PER=2 

* 8 : 9 * 

I |/~U • 

T 0 T :Q • 

MVPTOT =a. 

Mconoso 

CLOS=. FALSE. 

CO TO ^00 

Clio SETUP. REPIN Loop THAT PICKS tHACl straight line sechenteo jravlrsi 

2547 I F (NCORO. L T.NOO I GO TO 252 
WRITE <C,2Sl)NC0ftD,pvAL 

251 roBKAH;,2x,X5,11H SQUARES M«E ,PlO.S,2X,7MSHUTOm 

“TM5 r JS s i8« S58S.5??'S!3f s T !!82<,’Ste STMi-MWI-KC! «„.<» 

“^.T.McSM'clWPclfe. SVfl.W.r.Kiitt" *»«”* 

ao TO 730 

252 NCORO-NCORD* 1 

HVPTOlsHYPTOTbHVPC * 

HIPPS I NCORO) =HYPE 

CORO (NCORO.l t XXX SO 
CORDINC3RO,27=YYSO 
XO=XG*XXSC 
Y 0- YO* Y YSQ 

IT INEXtT.CQ. 1)L2X=L24-1 

IF INEX£T.£Q.4IL2Y=LZY-1 

IF(L2X.LT. JU0P.L7X.Gt.MAVX)CD TO 730 

XF ILZY »LT * 1 .OR .L2Y .G£ . WUVf l CO TO 730 

GO TO 700 
701 CLOSr. TRUE • 

-° ,L 1 Tr “ 1 

^C^NexCKGTl^HCNTCRzNCXSr-Z 


7 oo 


F {OtKl T . _ 

4(30 NtW©ur=3 

&XJJX8I6-X&»ID.FL0AT(L2X-1 ) 
IF (0*1 .LT. 0.70X1=3. 
irioxi mCT • XG RID) 0X1 =R GRID 
OVI=r«MG-VCRID«FL0ATIL;Y-l» 
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45 
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r 9 
it 
49 
»T 
51 

>r 

v*t 

*•5 

va 

59 

k: 
i ci 

in 

1=5 

act. 

u» 

i=? 

1=9 

ur 

!!' 

?, • o 

i a* 

life 
i : r 
t is 
1 i« • 

a:* 

t:i 

a;: 

121 

III 

i:t 

i;’ 
i:a 
1 = 9 
14 = 
l ji 

i n 

!*j 

a j- 
i j« 
145 
147 
1 38 

12 : 

14 : 

1*2 

\ll 

145 
lib 
1 «*7 
lie 

u ; 

113 
1 SI 
1 * : 
in 
it* 


tiu n 


j:L?V*tfNVU 

SUST EXIT POINT LOOP 

MAXI l liHMl IrfJl 
NAXIf , 

MAX 1 31 =HH IlM.JMl 
MAX ( 9 1 5 HH 1 1 * 1 f JI 
DO 901 H = 1 


If 111 ! IWS 11 !;«£«»« rfT?fev«»f?of 0 ‘ 

inMA*mi.lT.PWAUMAX»Ill:pYAL-PYOL 


410 

422 

4s* u 


424 

426 

43C 

435 


432 

*41 

H36 


CONTINUE 

mlifiHiMi. 

pil!! a, “ 4 * 

I?<III.CQ.4>I2=1 

statement below selects sides that have exit Points 

ir tMAXUl 1.LT.PVAL.AUD.HAXI121.0T.PVAUCO TO 920 

inBuMOUT.CO>2 1NM2SIII 

CO T 01 AJ 22 , 929 . 92 b, 92 »l, III 

OVE= I ipval-hax <1 11 / (MAX C 2 1 -MAXI 1 1 1 l*TCfiID 

SSff 8 -«» 

0X2; j ,-PV'At.=*f4ASi2S i/SHAXS 3 i -HAX i ? i i » *XGRIO 
OYifVSRIO 

$Y2=? IPVAL-HAX (41 >/<HAX(3)-HAXI9i 11*YGRID 

8 o 2 !o G S 5 o 

0X2;|»PVAC-HAX <1 I I / (ffAX(9 I -ItAX I 1)1 1* ACRID 

0X2=0. 

cior iiii.i 1 = 0 x 2 

c iBflWi ,a,=0V2 

Snless we h) wc null point square, humout shod u i with our at 0 x 2 , 
ir iNUHOwr.Aic.i ico to 932 
NCXFTsmN] 

rro 441 

<MUHour.eq. 2 »co to 93 s 


GO TO 444 
ir<AJUngur.„ 

WPITE <4,9 ifeiLlx ,LZv . M/HOur, pvaL | X gIG , TNIC 

w * v our, ^ i( " 91 ^ D,jr '°* 7 ' /> 

•^ 2 , 1 *f&fffi 0 Sy^^ 5 £ gSA?KuS??TSR. p ^ M t ^V 

434 IFIRGR IDU2X ,L2Y 1 .GT.l'lCO TO *92 
X IDrCIDE ( NMl ,11-0X1 

xiD=cIgciiarj 2 ,i)-oxi 
7 IP=CJDC INN2 ,2 l-OT 1 

wasrattunrair 

NEXETXhNl 
60 TO 919 
Ox2=CipriNN2,t 1 

i£lPL* •'«!*» 

[,Uf» = llS*MEMTER*WEXCr 
II=*^R 1 'ATlZ* ,LZV 1/1 


SI0CS3 




440 


0 Y 2 


CD TO SCO 


nil Lrxtll 1=0 
leac ti t =n 
LCxC < 5 » =0 
UkiiDua 

trie j la!: I % 

i? iLCsciiil i!eo.o.ano. n ,ue .riniir.awd. cioc m ,i i.oi.t -.s> . ano.cioe 
iU sri.or.i-.siico ro *tt 
*(■ gofchnwc. 


’ ■ to T# 415 

41* NIVETXII 

oSxscior tn,ii 

OVisCIDTl 11,21 


CKO WVPCRSOLIC CR10 SQUARE SI Cl ION 


NNS >5 IU0-0 xi* XGR IG »FLOA T ( LJfK -1 1 
V 1N0 >012* YQRI 0 *F L 0 A T • L£Y “II 
I r tRCillpI L?X jL 2 V I .CQ.ffliTRCfi 1D< iZX ,LZV > = 1 
XXSQ=X JRD-XBIG 
rYSO-VlND-YBIG 

MMuf® 

20Y2,x*s4,YVSC,HVPE,NENST,LCLiI,LCLY,NCqR0,XX,yV 

397 FCA4ATI/,2X,4HHERE,4I10 t /,2X»llFia.5,//,2X,4I10,5X,2F12.*,//» 

GO TO 500 
294 GO TO 250 

. LIRE SEGMENTED CONTOUR TRAVERSE NOW COMPLETE. NEXT. DIVIDf THIS 
* TRAVERSE INTO NRInC EGUAL SEGMENTS. THIS NUMHER IS FUNCTION NOT 

ONLY OF LENGTH OF TRAVERSE, BUT OF INCOMING ARGUMENT SMOOTH, WHICH 
CONTROLS DEGREE OF SMOOTHING DES1RE0. 

730 IFlNCQRB.LE.MAXCROiGO TO 722 
HAACRO=NCORO 
WHAT=PV*L 

732 NRINC :H VP TOT /M INC * l . 

HARC=HYPTOI/FLOAT(SAlNC> 

IFE.NOT.CLOS . ANO -NCORD.Cr. 1 , AMO .HR IRC .GT. I >60 10 f 34 
IF ICLOS. AMOiiNCORD.GT . 1 .AND .NRINC .01 .2 1 GO TO 734 
GO TO 5Q2 

•iEXT. SET UP ENTRY AND Ex I T SLOPE DATA FOR .FIRST SEGMENT fcEFORE 
. ENTERING MAIN CURVILINEAR INTERPOLATE AflO PLOT LOOP. 

734 XftCC=XX 

riicry 

IF INXINC«HANC. LT.. 75* OASMIX=. FA£SC« 


rticrv 

IFINRIMC 

M:0. 

> CNOrXA 
YfMO=YY 


00 740 I - l ,NCORO 
IF (R.MIPPSm.GE.HANClSO TO Th2 
M-H»HIPPI il I 
XfNO=XEND»CCRD<I, II 
740 TENDS VCHO ♦CORO 1 1 *2 1 

742 X?IHAnZ-HI/HIPPSII1 
XCNO=XEND«X*CORDII,II 
VC MO = VC MO »X*C0R0(',2I 

n«K=s|rt <xxso«xxso«ytSQ»VYS4i 
IF «HVPC.GE..o3plIOO TO 7S0 

HI pigt^n&ifniriuK!.. 


r WP**IP 4 • IF I n 

wnwier 


4F12.6) 


751 A(:XX 


4>co ro 7Si 

SCANS 

-SANG 

7E9 
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IF * W»m I PPS i J » . cc . H4NC )«® 1o 75f 


IRUWlH 

?SI XSIHAM>MI/MIPPS* J) 
»e=xe-4«coKDf j, j i 
Yb '. Yb>iJl «CORO < J 1 2 I 
757 xxS0=X6CG-xe 


YViavr»CG-va 

Nrj»£ ;sM T J XX SO *x «SQ • YY5QWY5Qt 
IF {MVFC.LT .<0001 ICO 16 71*7 

sinaacsvvso/hypc 

COS#*C r«XSe|HYPE 

A=.s*ATAN2 I SANGtCOSBAC -CANC*SI 


SArSINlAl 
CXSC05I* ) 


ANG *C 05BAC -C A WC *S I N | * Ct C*KC*COt »-»C* i Af/D’li IM&ACl 


scwr=siMrAc*CA >cos§ac*sa 

CCMTrC055*C*CA ~S INBAC tSA 
CNTI5=ATAIj2ISCNI,CtNT> 

ccf8|=cc8r 

C fSTCfl HAIM CURVILINEAR INTERPOLATE A NO PLOT LOOP 

?F ((SPc^N^^N&f N?>SP TO 762 
IFICLOSIGO TO 760 
SCUT=SANC 
COUT-CANG ' 

SO TO 200 
760 SuUT-|>t«ii 
COuT SCCCNT 
GO TC 200 
762 XIND=X* 

YINOSYV 

2 AMC =WA AIC *7 LOA T * LOPE « I ) 

HlO. . 

00 IBM Irt.MCORO 

IF CH*HLPPGiri.fiE.2AMC>Co 10 IBB 

M=M*HIPPSUI 

. x i«c=x imd*co«oi i ,1 i 

H* T I MO s Y I NQ' LORE) ( I • 2 I 
SO TO 760 

704 X s (ZAMC-U i /HI P PS (I I 
XTMO-X IMP *k»COPOI 1 » 1 1 
VlNO=y INP>X«CCPO( I ,21 
766 XXS0rxlR5-XfN0 
VYSQrYIMO.rCMD 

Hv7t:jQRT I XKS(i *X *5S*V Y50* T YSG I 
IF IHVPE.LT. .OOQIIGO TO 767 
SIAF ORiYY 50/HYPE 
COSFOR=XX SO/HYPE 
MTPFORSHTPE 

* > t»?AT*M2 IS1NFOR*CAI|G“C05FOP*S AMS t COLFOR»C AoG *SIN70R*SANGI 

j A«SIn (A I 

CA SCOS I A I 

SflUT :S ARC*CA ♦CANG*SA 
COOT=C AHG«CA-SANG*SA 
CXCT-ATANXlSOuTf COUTt 

zoo 

Jf*||M^Lf55csi?&ffcJ5? & T .E.5T .OR.YLMOoLT.SOUTH.OR.YCMO.ST. WORTH) 

IF IlNCS l«0 TO 444 
INCS :. TduC . 

CALL PLOTlMCG-WiSr.vaCS-SOUTH.Jt 

440 1=IPFR „ 

IFifMSHlX .AN0.I .|0.2J IiIDpLUT 

1 S^t=£uOC .YCM 0 - 50 UTM. X A 

TOT=TOT«HfA 
60 TO 790 
C 

C KSIN MAH INTERPOLATION FOR SEGMENT 
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IN 

325 

324 

32’ 

ill 

32« 

Hi 

>33 
3 3* 

m 

33’ 
3 3! 
7 3» 
3 "i.’ 
3 h l 
3 «.2 
3.3 
3***i 
345 
344 

111 

3*9 

353 
3*1 
3«: 
352 
JE* 
3«5 

354 

ill 

K* 

342 

ni 

745 
JU 
34’ 
34A 
St9 
3 73 

m 

373 
3 7* 
375 
7 74 
377 
174 

K! 

ill 

a 3 

2e* 

3 e 5 
lo* 

3 e 1 
129 
191 

15 ] 

3 >3 

!»•. 


- 4(20 COXfMC or TRUE iNTXY.AfJCLC 


eiCAKCl 1 « 

SEMI A NO CENT ARC 
sour AND CGU 7 ARE 


HSS 


and cosine or rituc exir anc 


5 * |N 6 AMD CANS ARC SItiC AMO COSINE OF STPAJCMr LINE EnTRY/TAIT CONNECIOR 

BEGIN SECTION* IMA t INTERPOLATES AMD PLOTS THRU 5 CCMCNT 

<W0 5 :SEWT *CANG-C£NT *5 ANO 
C:CCNT*CANG»SCNT«SANG 

c:nc=cang*c«sahg»s 

S:SOUt«CINC-cau 7 « 4 IMO 

Ei'TIItffSBSWlTM?. 

tx:iBc« 

TV- 76 CG 

?}«: 8 ..c 

H»PU*IX=NYP-CLI I 
H25:.«S*HVP 
•*5u T VP s I YP«QINC 
OlsTVP 

inoi.or.H25)oi-H2& 

ORsT VP-M25 

i|NG-cfe^i^”C?*f9p*t2* 102-Dll 

TX=TX*QTNC*COS ISZNGI 
TV=TV«QINC*SIN<SINGI 

ENO snake INTERPOLATION SECTION— -TRY AND FIGURE IT our AND CO NuTS- 
NCXT STORE POINTS THRU THIS SEGMENT TCP FINAL ADJUST AMO PLOT 

IF ININC.LT. 27S ICO TO *53 

MRIT£lb s *5«lPYAL 

m 


b 4*5“ }PY AL . -SCC , VBES .XCliD .TEND 
I/.2X.12HN1NC SKUiOFF •2X.&FJ2 .21 
.. .^yfoPF MESSAGE RECEIVED MFRC. ARRAY XXPLGT is too small. 
roR INFREQUENT MESSAGES. OONT Uo4ry ABOUT IT. SINCE LACK OF 

CLOSURE IS ADJUSTED OUT. IF MESSAGE PERSISTS, EITHCP INCREASE 
_ PLOT INCREMENT LENGTH OR SIZE OF XXPIOT. 


..... CO TO *55 
453 NINC:NlNC*l 

XXPLOTiNInC.IIsTX 
XXPL0T(NINC,2l=fV 
IF IT YP .LE .HfPMAx ICO TO 450 

adjust- for closure error, then Plot curve among mis secment. 

NCRS tXCND-XXPLOT IN INC, 1 1 1 /FLOAT INI NCI 
YCRr I YEND-XXPLOT IN INC , 2 1 ) /Ft OAT ININCI 
UlQ. 

nwBcso 


403 


<04 


BEGIN SECMENT PLOT LOOP- — -DASHED OR SOLID CURVES 

SUBROUTINE ENOCH IS CALLED TO LABEL LINES 

00 410 I=I,NINC 

u:(!«x£ft 

V=V*VCR 

X=XXPLOTII ,1 )«U 
VSXXP L 0 T(I. 2 I«V 

!?AjrJ:y£ S :A 0 ?A¥'* T * Ei ‘ sr * 0 ** T * L ^ •SOUTH. OR. V.GT.YoRTMICO TO AOS 

ISiiKttSL 10 603 

IF l IpCK.CQ. 2 )C ALL PLOT I XBEC -WES T , YBfC -SOUTH, S> 

?&=i 

AURsRPP -WEST 


*ux=xpp-vest 

WettT » 


Jo ISO 10 4 Q<I 

.EC.OICALL cn 5 crixux,yuy>pval,i> 


IF IDASHIX.ADO.IPCR.LO. 2 iCO 10 415 
<15 call plotixux.yuv.joploti 
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3 tl 

3*4 

4 22 
HZ\ 

Sb! 

oct 

4GS 

4ri 

‘•S' 

4. i 
‘oC9 

on 

4 
4 
4 

4 14 
* 4 l? 

. m 

18 

019 

420 

421 

422 

423 
42* 

4 25 
Hit 
42? 
428 

5 29 
4 3 “* 
431 
4 32 
4 32 

4 3* 
4 J4 
4 2b 
4JT 
4 28 
434 

4 *; 

4 4 I 

442 

443 

444 
4 4 5 
44b 

447 

448 

449 
.i»9T, 5 



M co Tb tif 


TfiNHkWO.Gl.aiCO 10 
I0PL0T=1 
NHAROSLDASHI 
„ GO TO bof 
5 JO NSCF T shsor 1* l 

rr iMsoFT.cr.oico to 409 

10pLOT=2 
NWT=tO*SH2 
409 iflk: 2 


40f 


MUNCrMUMC* 1 

IF I I .Wpr.OOLlOS .OO.YuY.Ut .YMO* *40 


THAX : TUT 
XMAX?XV* 

fW ?! 410 


r* HO 


11 A0 * kit^fcolo u -irS? f »?soo?m , 

I lb I AOl 1 T-14 


no cSMiif 


mi ( n 


EHO AOJUST AMO PLOT SECTION 

AIU6G :NliOG *fAUNC 
TOT :TQT»OINC*FLOAT<HJNCI 
rr IT0T.LT.THAX1G0 TO 790 


_ WAITM4,A|b2|PVAL ,TMAX 

442 ; ,J AH A n / * 2 X f 1 6 H fi £ A C H E 0 THAX ON ,2F12.4| 
HTPMSHVPFOK 


190 


<78 ,1 ?rq.NRINCICO TO 800 


yog 


m\n$\ 

EMTlft-lxC T 

XMLGSXEND 

■florrtwo 

;lHOSAlHD 

CNQrriND 

omTINOC 

CUC, HAIM CURVlLlHEAff INTERPOLATE AnO Ptor LOOP 


i.rHDilL 


*°° $p£ i xiux*!iHt 0R ** RaT * D0L * aSc0R ’ VH *** ll ’ T * *°i ,C8 To 501 
YPP: THAX ♦SOUTH 

„ CALL CnOCR I XP.AX • YllAX|PVALj2) 

HOS I WC =WUGG 

vallim^pval 

502 RE TURN 

two 
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• 1 • ^PSfe^K'f&Ki^kv^U' u* •«, , 

COAMOM /QINDCQAHONUA.WlOl ^OCH.qUS^yUS 
OIaENSION 013) 


THIS SUBROUTINE IS CtUCo TO lABCL C0NI0U«4 
* • 

DurSORTl IX -XL »!(•(* -XUS I ♦ IT-YLASl*IV>TLA 5 *» 

IT (0H.Ll.2.*H0Nur».0»,«C«U/1.LE .0. ICO TO 25 
JJjrC 

inico.ra.2i4o to ik 

lSIT-HOONI 
IS < x-u 00 t| 

(sly i 

H=1 


it libuttt 

Otl»=*Bb I 
K = 1 

ohuisui 

??,!?■ I Ice.onGO to 10 


DM:0 1 1 1 
•Cl *1 


10 CONTINUE . , 

CO t 6( 1< ( 1 k • 14 • Ik i .* 

12 VA0=-MCMgH/2. 

60 *T0 105 
Ik r«D=.o? 

60*TC 130 

u npi.oi 

t AO* -HONUM/2 . 

60 tO 20 

1* j 4 ||: -hOWNH”. 02 

60 TO 100 
2C XASCsXAQ 

r ao=v*vao 

JJJ:JJJ«1 

inJUJ.E0.2IC0 TO 333 
CO TO H4 

333 CAtL *JOWeKUAG,VAQ,M0NUN,OVAi. ,2/0.,2» 

„ jjj=0 

334 CONTINUE 

CALL PL0T(k,y,3l 
ALAS:* 

„ »UAS:y 
JS RETURN 
IDO X AD: - . 7* •NONUH 

ITIPVAl .«f .9.5.04.PVAL .l« . I » .5 I I XADlXAD-HOIfuM 

} nPVAL.CE.99.5.0R.PVAL.LE.|-9.S) IXAD:X AO-HOMUH 
f IPX AL.tE .999 .5 >*R .PVAL.LC.t -9 9.*) IX AO=X AO-HONUH 
Ml .CQ.t I XAD-.S*X AD 
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END 
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SAMPLE PLOTS 
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RUN NO i LOO 7. 

DISCHARGE VELOCITY i C.84CM/SEC 
DISCHARGE TEMPERATURE! 18.4*C 
NIND SPEED (MAX) > 4.SIM/3EC 
CURRENT ( UQCA3EE FLOW ) i 4.8 CM/SEC 
TOTAL SIMULATED TIME i 1-01 HRS 

«10 l 

LENGTH SCALEC METERS) 6 i 



ISOTHERMS RT K= 1. 

LRKE KEOWEE-I RIGID-LID MODEL) 
SIMULRTIONS FOR FEB. 27 1979 
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RUN NO* LOO 7. 

DISCHARGE VELOCITY i 0.04CM/SEC 
DISCHARGE TEHPERATUREi 10 .4*C 
HIND SPEED (MAX) i 4.51M/3EC 
CURRENT ( UQCR3EE FLOW 1 a 4.0 CH/3EC 
TOTAL SIMULATED TIME « 2.01 HRS 

MlQ 1 

LENOTH SCALE! METERS) 0 «-' 00 6 . 1 ,0 ° 



ISOTHERMS HT K= 1. 

LAKE KEOHEE-C RIGID-LID MODEL) 
SIMULATIONS FOR FEB- 27 1979 
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RUN NO i LOO 0. 

DISCHARGE VELOCITY * 7.42CM/3EC 
DISCHARGE TEMPERATURE* 31 .7 # C 
HIND SPEED (MAX) a 3.09M/3EC 
CURRENT ( JQCA39E FLOM)a 1.1 CM/SEC 
TOTAL SIMULATED TIME a 1-01 HRS 


LENGTH SCALE! METERS) 


0.00 


VELOCITY 3CALECCH/3EC) 


0.00 


* * * \ 


VELOCITIES AT K= 1. 

LAKE KE0WEE-C RIGID-LID MODEL) 
SIMULATIONS F0R FEB. 27 1979 





VELOCITIES RT K= 2. 

LAKE KEOWEE-C RIGID-LID MODEL) 
SIMULATIONS FOR FEB- 27 1979 
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VELOCITIES AT K= 3. 

LAKE KE0WEE-( RIGID-LID MODEL) 
SIMULATIONS FOR FEB. 27 1979 
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RUN NO a L0O 0. 

DISCHARGE VELOCITY i 7.42CM/3EC 
DISCHARGE TEMPERATURE i 31.7*C 


HIND SPEED (MAX) a 3.09M/3EC 
CURRENT! J0CAS3E FLOW ) a 1. 1 CM/SEC 
TOTAL SIMULATED TINE a 1.01 HRS 



Velocities rt k= 4. 

LAKE KEOWEE-C RIGID-LID MODEL ) 
SIMULRTIQNS FOR FEB. 27 1979 
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VELOCITIES RT K= 5. 

LAKE KEOWEE-C RIGID-LID MODEL) 
SIMULATIONS FOR FEB. 27 1979 
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